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In this Isswe 

The voice on the telephone one day last year was HP statistician Tim 
Read's, Pointing out that not ail technical contributions are made by en- 
gineers, he made us an offer we couldn't refuse — to invite a group of HP 
statisticians to submit papers on applications of statistics. The results are 
five excellent papers that cover an impressive range of topics: designing 
electronic hardware, calibrating instruments, setting specifications, estimat- 
ing reliability, and managing software projects, In "Statistical Issues In Setting 
Product Specifications," page 6, we learn how the science of statistics helps 
ensure that specifications are useful for describing the products perfor- 
mance, relevant for competitive comparisons, measurable and verifiable, and valid over the 
required environmental ranges and time intervals. This paper takes a tutorial approach to Its 
subject. 'Circuit Design Using Statistical Data Analysis/' page 12. tells how measurements on 
prototype amplifiers have been analyzed statistically to guide the design of adjustment-free 
amplifier boards. "Statistical Calibration of a Vector Demodulator,' 1 page 18 h describes the design 
of the statistical demodulator calibration algorithm that is programmed into the HP 8981 A Vector 
Modulation Analyzer, The algorithm can be used to calibrate not only the built-in demodulator but 
also external demodulators, "An Availability and Reliability Tooi for Computer Systems/' page 
26, describes a software program that computes system availability based on reliability data for 
the system components. The program helps R&D engineers make trade-offs in designing for 
reliability. "Project Management Using Software Reliability Growth Models,' 1 page 30, relates one 
HP Division's experience with a model-based method for predicting the duration of software 
system testing. The model also helps management decide when to release a product. 

Half-inch reel-to-reel tape recorders are widely used for archival storage of computer data, 
backing up data stored on disc memory, and interchanging data between systems, The HP 7980A 
Tape Drive, besides being faster and more reliable than its predecessors, is interesting because 
it's an autoloading drive. Designed to mount in a rack with the reels and tape path in a horizontal 
plane rather than the more conventional vertical plane, ft takes up much less rack space, but its 
reels and tape path aren't readily accessible for the user to mount reeis and thread tape. However, 
no problem. All the user need do is open the door and slip in a reel of tape. The HP 7980A 
automatically threads the tape onto the takeup reel. How does it do that? Its designers tell how 
in the article on page 36. The drive's state-of-the-art control electronics (page 43) provide features 
such as variable-velocity rewind, which chops 25% off the 120 seconds it used to take to rewind 
a 2400-foot reel of tape. 



i 
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In the February 1 988 issue, we presented the design story of a family of digitizing oscilloscopes 
based on waveform recorder technology . The primary objective in these instruments design was 
maximum fidelity within a given bandwidth. Featured this month is a digitizing oscilloscope designed 
for maximum speed, which means maximum bandwidth for capturing singie-shot phenomena. 
The key element in the HP 54711D Oscilloscope is the advanced one-gigasample-per-second 
analog-to-drgital converter (ADC) shown on the cover. Thanks to this ADCs speed, the HP 541 1 1 D 
can capture and display signals containing frequencies up to 250 megahertz with a single look. 
This wide single-shot bandwidth is needed for digital system diagnostics, high-energy physics, 
transient analysis, and other applications After a brief introduction (page 58) by Joe Millard, who 
was project manager for the HP 54111D's development, three authors from HP Laboratories 
describe the fast-ADC technology (page 59). Four digitizers are interleaved in this design, which 
uses gallium arsenide, silicon bipolar, and silicon NMOS integrated circuits and thick-film hybrid 
circuit technology. The wideband hybrid preamplifier and attenuator circuits at the instrument s 
input are described in the article on page 67. Filtering, bandwidth, noise, and waveform reconstruc- 
tion issues are discussed in the paper on page 70, which shows how the six- bit ADC can provide 
seven or eight bits of information for sufficiently low-frequency input signals. 

-R.P. Dolan 



What's Ahead 

In the August issue, eight articles will discuss the design of the HP PaintJet Color Graphics 
Printer, with particular emphasis on the design and manufacturing of the print cartridge Aluminum 
gallium arsenide (AIGaAs) red light-emitting diode lamps are the subject of another article, HP-RL, 
the Hewlett-Packard Representation Language, was used for expert systems research in HP 
Laboratories until recently. Well have a paper summarizing its contributions and another paper 
on MicroScope. an expert system for program analysis, which was originally written using HP-RL. 
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Statistical Issues in Setting Product 
Specifications 

A primer on the use of statistics in specification setting. 
by Sherry L Read and Timothy R.C, Read 



PRODUCT SPECIFICATIONS are frequently the first 
chance that a potential customer has to evaluate a 
new product, The purpose of such specifications is 
to characterize the product's performance and to warrant 
a minimum performance level to the customer. Con- 
sequently, il Is exlremely important that the specifications 
accurately reflect the quality and performance of the prod- 
in :L Conservative specifications that downplay the true per- 
formance are likely to benefit competitors selling similar 
products who provide a more realistic assessment of the 
specifications. On the other hand, specifications that are 
too ambitious for the product are likely to cause high man- 
ufacturing cost and possible customer dissatisfaction. 

The choice of product specifications is driven by cus- 
tomer needs, product design, competition, and manufactur- 
ing capabilities. Final setting of specifications should ac- 
count for all these aspects, as well as for statistical evidence 
of product performance from a production pilot run. In 
particular, the process of setting specifications needs to 
reflect the accuracy and repeatability of a single production 
unit and the variation between different units, 

Examples of typical product performance questions are: 
■ With what level of certainty will the harmonic distortion 
of this signal generator lie within a certain limit? 
With what confidence will the noise power of a single 
100-MHz crystal oscillator vary by less than 10 percent 
over a two-week period? 

What is the probability that a product meeting test 
criteria at room temperature will meet or exceed specifi- 
cations over all warranted environmental conditions? 
In this paper we illustrate how statistical methods of 
data analysis can help product designers, manufacturing 
and marketing engineers provide efficient and accurate an- 
swers to such questions, 

Objectives 

What are the objectives for the process of setting specifi- 
cations? 

Specifications should be: 

Useful in describing product performance to customers 

and field offices 
g Relevant for competitive comparisons 

Measurable and easily verifiable in production at service 

centers, and by customers 

Valid over the appropriate environmental conditions 

and time intervals. 

To meet these objectives and provide maximum inf orma- 
tion to the customer, product specifications are frequently 
presented in a two-tier structure,* The top tier describes 



warranted specifications and the second tier describes sup- 
plemental characteristics. The supplemental characteris- 
tics provide information useful in product applications by 
giving typical, but nonwarranted t performance parameters 
for the laboratory environment. In this paper, we shall con- 
centrate on the setting of warranted specifications. How- 
ever, the same general methods can be applied to supple- 
menial characteristics. 

Statistics, the science of drawing conclusions and mak- 
ing decisions in the face of uncertainty, aids in achieving 
many of the objectives for specifications. Based on produc- 
tion pilot units, statistical methods are used to estimate 
the expected performance of future production units. En- 
vironmental effects can be modeled so that not all units 
have to be tested over all environments, Specific graphical 
tools provide a clear picture of performance to facilitate 
the decision process by marketing and manufacturing when 
setting the final specifications. Statistics is the key In under- 
standing and efficiently quantifying all sources of variation. 

Statistical Model 

The underlying model for the specifications of many 
Hewlett-Packard products is [Must rated in Fig. 1. In I his 
diagram, production margin is the difference between the 
average product performance and the test line limit [TLL], 
which is the pass/fail limit used by the production line at 
final test under standard environmental conditions, Delta 
environmental (DE) represents the possible change in per- 
formance of the product over the environmental extremes 
specified on the data sheet (e.g., temperature and humid- 
ity]* and drift represents the change in performance over 
the calibration period or other specified time interval. Mea- 
surement uncertainty (MU) accounts for possible measure- 
ment errors in the equipment used to characterize the prod- 
uct. The customer guardband represents any additional 

"Ttiis structure was originally proposed 1g the Nation a< Conference tfl Standards 
Laboratories {U S.A ) by John Wmck In 1977 
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Fig. 1 . Representation of the components that make up a 
warranted product specification. TLL is the test fine limit used 
by the production line at final test 
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guardband considered necessary to ensure thai in a worst- 
case scenario, any customers unit will perform better than 
the published specification [this additional guardband is 
sometimes set to zero). The right end of Fig. 1 (SPEC) rep- 
resents the warranted customer specification as described 
by the product data sheet. 

When finalizing product specifications, it is necessary 
to look at Fig. 1 from two directions: from the SPEC inward 
and from the observed pilot data outward. The former di- 
rection involves setting internal criteria (the TLLs) to 
guarantee that the product performance meets specifica- 
tions over all warranted conditions. The latter direction 
involves setting specifications that will allow a high yield 
in manufacturing and simultaneously minimize the chance 
that good units will be rejected and maximize our confi- 
dence that the units passed will perform as specified. We 
shall focus on the statistical aspects of moving from the 
production pilot units outwards to the preliminary specif i~ 
cations. However, it is important to realize that the final 
specifications are set by iterating between both approaches 
to ensure that the final product has a high yield and satisfies 
the market needs. 

Using our experience of specific product introductions, 
we will describe the relevant statistical methods needed 
to complete the estimation of each of the quantities in Fig. 
1 based on production pilot data, The same approach is 
also used to adjust specifications on more mature products. 

Estimating the Test Line Limit 

Based on the assumption that the production pilot units 
are representative of the product characteristics that can 
be expected in production, the goal is to use these pilot 
units to calculate the TLL so that the vast majority of pro- 
duction units will be better than this limit. In particular, 
we need to ensure that the products will be readily man- 
ufacturable over the long term a! a low cost. 

If we use the convenient assumption that the production 
pilots come from a Gaussian [normal] distribution wilh a 
population mean of ^ and standard deviation of u, then the 
interval (/* - 2a + ^t + 2(j) contains 95% of the distribution 
and [jjl - 3rj,jjt + 3cj] contains more than 99% of the distri- 
bution. This normality assumption allows us to model the 
expected production pilot distribution and the correspond- 
ing percentiles using only the two parameters }i and it. 



Robust Estimators 



The sample mean and sampie standard deviation are ihe best 
estimators of the focaiion and scale of a Gaussian or normal 
distribution However, a few farge errors may occur m the data, 
«g fn a distribution with heavier tai*s than the Gaussian In 
this case, estimators that ptace Jess emphasis on extreme values 
may provide more appropriate estimates for the mean and stan- 
dard deviation 

One such set of estimators consists of the trimmed mean and 

the Winsonzed standard deviation To calculate these estimates 

trim off g observations from either end of the set of n ordered 

1 observations x 1 < x 2 < < x n , where g Js the ciosest integer 

to (0 1 )n Then the 1 Q% inmrned mean is given by: 

and the Winsorized standard deviation is calculated as 

*h - {[(g + IHx^, - x,)* + (x g , 2 - xf + 
+ <x n _ g _--x r ) ? + (g + 1)(x n _ 9 - xf]f(n-2Q-^r 2 . 

The Winsonzed standard deviation formula treats the large (or 
small) observations that were trimmed as if they were equal to 
the largest ( or smallest) observations in the trimmed sampJe. 

When the trimmed mean and Winsonzed standard deviation 
are used, the effective sample size is the size of the trimmed 
sample. Hence, when finding the K factors for use in x T + Ks,. 
enter the tables with n = n - 2g. 

Other robust estimators, such as the brweight, are described 
in References 1. 2. and 3 listed on page 1 1 



While there are good reasons to assume normality, it is 
important to check the valiilily of this assumption before 
proceeding to estimate the percentiles of the production 
pilot distribution. Histograms are frequently used to verify 
the symmetry of a distribution. However, the investigation 
nf normality is much clearer when the data is plotted on 
a Gaussian quantile-quantile piot as illustrated in Figs. 2 
and 3, 

Here the data values are plotted against the values that 
would be expected from a Gaussian distribution. If the data 
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Fig. 2* Quantite-quantiie plot of spurious distortion data on 
a voltage scale 
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is Gaussian, the values should lie approximately on a 
straight line. Such plots as these help to check assumptions, 
to determine the appropriate scale in which to set specifi- 
cations, and to detect measurements or units that appear 
not to conform to the bulk of the data, To illustrate the 
appropriate choice of scale. Figs. 2 and 3 represent data 
on spurious distortion for an HP Stanford Park Division 
product on both Ihe linear (voltage) and the log (dBc) scale, 
respectively^ Fig, 3 shows that the measurements on the 
log scale are more nearly Gaussian than those on the linear 
scale in Fig. 2, Consequently, the calculations for determin- 
ing the TLL of this product should be done in dBc, 

A different situation occurs in Fig. 4 r which represents 
a quantiie-quantile plot of intermodulation distortion from 
some pilot-run models of the same product. While the bulk 
of the data appears to be nearly Gaussian, one instrument 
(circled) appears to be different, exhibiting extremely good 
performance, This outlying instrument should be ac- 
counted for in the estimation of the mean and standard 
deviation used to calculate the TLL. The following section 
outlines the methodology for dealing with such points. 

Estimating /i and n- 

Consider now replacing the mean fj. and the standard 
deviation rr with their standard estimates 



x = 



+ X. + 



and 



s = {u%-x) 2 + (x 2 -x] 2 



+ x n )/n 



+ [x n -x) 2 ]/(n-l)} 1/2 , 



respectively, where Xj represents the measurement value 

from pilot unit i (i = 1,2 nj and n is the total sample 

size, We are faced with two potential problems. First, con- 
sider the situation where one or two extreme points are 
found in the sample, such as in Fig, 4 for the intermodula- 
tion data. If these points occur only at the opposite extreme 
from where the TLL is to be set, then under certain cir- 
cumstances we may wish to use estimates of u. and a that 
are less sensitive to (i.e,, more robust against) individual 
extreme departures from a Gaussian distribution. The stan- 
dard estimates x and s given above are notoriously non- 
robust. An example of how one extreme point can throw T 
off x and s is illustrated in Fig. 5 below, Here both sets of 
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Fig . 4, Quanvte-Quantite plot of two-tone intermodulation dis- 
tortion data, showing an extreme point. 



Observation Set 1 
x = 5.0 s = 3.3 r 



Observation Sei 2 
x = 50 s = 3.3 n = 11 




(a) 



10 



15 




(b) 



10 



15 



Fig, 5, Both sets of observations have the same values for 
x and s, but the values for (b) are inflated because of the 
extreme point at 75 

observations have the same values for x and s. However, 
clearly, x and s are inflated for the set of observations in 
Fig. 5(b) because of the one extreme point at 15. 

There are a number of procedures for robust estima- 
tion, 12,3 For example, the standard estimates x and s for 
the intermodulation data in Fig. 4 are —63.53 dBc and 4.59 
dBc, respectively, whereas the more robust 10% trimmed 
mean and Winsorized standard deviation (see "Robust Es- 
timators," page 7} for this data are -62.36 dBc and 1.64 
dBc, respectively- When such methods are used, the en- 
gineer must be certain that the extreme points are well 
understood, and could not equally well have occurred on 
the other extreme of the distribution, Unless the mecha- 
nisms causing such extreme observations are understood, 
robust estimators should be used with extreme caution in 
this context. 

The second potential problem in replacing ju and o by 
their estimates is that probability statements such as "the 
interval [p. - 2rr, /x + 2cr) contains 95% of the distribution" 
are no longer true in general. The amount of the underlying 
distribution covered by the interval (x — 2s T x + 2s) or 
(x — 3s t x 4- 3s] depends on how accurately x and s 
estimate ju, and cr. which in turn depends on the sample 
size n. For large samples, the statement "the interval (x - 2s, 
x + 2s} contains 95% of the distribution" will be approx- 
imately true (as a result of the Central Limit Theorem). But 
how large a sample size is needed to achieve reasonable 
accuracy? 

This question is answered by reference to suitable toler- 
ance limit tables. 4 5 Given the sample size, these tables 
provide a level of confidence that the interval does indeed 
cover the requested percentage of the distribution. For 
example, based on a sample of 10 units, if we want an 
interval that covers 95% of the distribution with 95% con- 
fidence, we wmdd have to use (x — 3.4s t x + 3 + 4s) rather 
than (x - 2s, x + 2s), This tolerance interval procedure 
accounts for the random variations in the estimates of both 
ft and cr, and provides an interval that will include at least 
a specified percentage of the population w r ith a high degree 
of confidence. 

In conclusion, our procedure for setting the test line 
limits involves three steps. First, verify the distributional 
assumption of a sample from a Gaussian distribution. Sec- 
ond, estimate the mean and standard deviation of the per- 
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Propagation of Error with 
Multiple Sources of Variability 



For a single sc Mton cr can 

be estimated by 

s = «<x, 

where the* = 2 , n) represent measurements on n units 

However, it frequently happens that the test characteristic ss 
measured as a known function of other maependem measured 
variables For example, suppose that y^ , y 2 , and y 3 are the heights 
of three types ol discs. If one stacks the discs on top ot each 
other then the height ot the stack is h ■ = y T + y 2 4 y 3 . Given 
tolerance limits on the individual discs, what are the tolerance 
limits on the overall stack height? Define s r , s> and s 3 to be the 
standard deviations of y t , y 2 , and y 3 . based on n,. n 2l and r\ 3 
data points, respectively Then the law ot propagation of errors 
states that the standard deviation of h is given by 

% = Vsf + si + s^. 

Note that this formula requires only that the measured variables 
y 1 , y ?1 and y 3 are independent of one another Further distribu- 
tional assumptions are not necessary. 

Tolerance intervals can now be constructed for h according 
to the methods discussed in ihis paper, using an effective sample 
size of n = n, + n & + n 3 - {3- 1). 



foxmance using a standard or robust procedure. Finally, 
using the K factors for tolerance intervals (based on the 

production pilot run sample size), set TLLs at x + Ks or 
x — Ks as appropriate. These limits will provide the first 
step in setting customer specifications and provide knowl- 
edge, with high confidence* of what can be expected from 
future manufacturing, assuming a stable manufacturing 
process, 

Estimating Delta Environmental 

To guarantee that the product will meet or exceed its 
warmn ted performa nee spec Lfical ions over al 1 e n vironmen- 
tal conditions, il is necessary to estimate the expected 
change in product performance from the standard test con- 
ditions to the worst-case environment. In most cases, this 
can be estimated as the average measurement in the worst- 
case environment minus the average measurement in the 
standard (or test) environment, that Is* 



DE 



*Aid< 



Some products will experience not only a mean shit I 
over an environmental range, hut also an increase in the 
variability of their performance. This increase in variahi I i l \ 
can be accounted for by the following equation: 



DE = x„ 



x sui + ECfccpv " B *dJ* 



Fig. 6 illustrates both situations by depicting the distribu- 
tion of a population of units in the standard operating 



environment and in the worst-case environment If s tmv = 
s^, as in Fig, 6a, then DE will be equal to x«. nv - x 5 < d , If 
Scnv is greater than s fitd (Fig. 6b). then DE will be equal to 
the difference in some percentile. By accounting for DE, 
customers can be assured that the units tested to the TLL 
in the standard environment will meet or exceed specifica- 
tions over the warranted > ondil 

Estimating Measurement Uncertainty 

Measurement uncertainty accounts for possible changes 
in measurements caused by the inherent variation of the 
test or measurement equipment. The most common method 
of estimating this variation is to combine the specifications 
provided on the test equipment data sheets for each piece 
of equipment integrated into the test system. For example* 
if there are four pieces of equipment making up the mea- 
surement system, with individual specifications • r^. ±r 2 , 
-t 3 , ±r 4 . then typically we can estimate the respective 
standard deviations s P s 2 , s :|t and s 4 as one half or one 
third of these values (it they are given iis 9ri"o fir 'W \, limits, 
respectively). Using the root-sum-square [RSSJ principle 
(see "Propagation of Error with Multiple Sources of Vari- 
ability," this page), the measurement uncertainty for the 
system is calculated as 



MU = V Is, 2 + s 2 3 + s ;i 2 + s 4 % 

Setting Preliminary Specifications 

To account for all the sources of variation, the prelimi- 
nary specifications relative to manufacturing capabilities 
would be defined as 



SPEC = x + V (Ks) 2 + (3MU) 2 + DE. 



or 



SPEC = x . - V(Ksf + (3MU) a + DE. 





When s- 



(a) 



Std 



Env 



A 



X\ 



When s_ nv >s* 



Std 



(b) 



— i — 
Env 

-DE- 



Flg, 6, Estimation of defta environmental (a) when the mean 
ts a function of the environment but the standard deviation is 
constant, and (b) when both the mean and the standard de- 
viation are functions of the environment 
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Fig. 7. Boxptots of common- 

mode rejection data for different 
instruments and channel numbers 
over aii environmental conditions. 



as appropriate. Again we use the RSS principle to combine 
the unit variation (s 2 ) and the measurement uncertainty 
(MU 2 ). Delta environmental is treated as a bias. This is the 
preliminary SPEC based on evidence of manufacturing ca- 
pabilities, but these manufacturing capabilities do not cap- 
ture all the can si derations in choosing the final specifica- 
tions. Much of the final decision must be based on cus- 
tomer, competition, and business considerations. 

Graphical Tools 

To facilitate discussion and understanding of the pre- 
liminary specifications, graphical displays of the pilot data 
with TLLs and preliminary specifications are useful. In 
particular, boxpiots {see +k Boxpiots , ,+ page 11 j are impor- 
tant tools for displaying the magnitude and dispersion of a 
number of sets of data. For the HP 5183 A Waveform Record- 
er, the dispersion of measurements for each pilot run instru- 
ment over the various environmental conditions is shown 
in Fig. 7, and the aggregate pilot run data for each of the 
tested environments is shown in Fig. 8. Fig. 7 allowed the 
group responsible for setting specifications to identify the 
particular channels that performed exceptionally well [e,g., 
#250.1) or showed little variation over temperature and 
line voltage (e>g. t #247.2). Fig. 8 indicates that there is little 
difference between environments. These graphical tools 
help identify cases that are different from the rest, a neces- 
sary step in improving the consistency of products and the 
manufacturing process. 

Boxpiots can be used to evaluate specifications over a 
range of measurements. For example. Fig. 9 is a plot of the 




Environment 

Fig. 8. Boxpiots of common-mode rejection data for different 

environments for all instruments. 



pilot run results of the same measurement on the various 
HP 51 83 A Waveform Recorder ranges. This kind of graphi- 
cal display helps determine the ranges where the specifica- 
tions should change. If there are any ranges that are particu- 
larly variable, they can be pinpointed for more intensive 
testing, providing a basis upon which to optimize the effi- 
ciency of the internal testing process. 

Once preliminary specifications are established, many 
of [he specifications will he as expected and will require 
little if any discussion, while other specifications may re- 
quire some specific action. A systematic approach to cal- 
culating the preliminary specifications is important. This 
includes validation of assumptions, a clear understanding 
of the sources of variation, and graphical displays of all 
the pertinent information. This approach will do much to 
ensure that the decisions about final specifications result 
in a valid representation of a product's performance upon 
which customers can depend, 
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Boxplots 



Boxplots are a graphical tool displaying (tie salient features of 
on of a set of data— name Sion, tfse spread. 

any a :tnd outlying data ponts Because of their com- 

pactness, boxplots are espec*& - ng several 

sets of data on one graph 

t shows a h; slog ram and a boxpior of a single set of data. 
In this boxplot, the maximum value is not joined to the t: 

ni S ker (i e the line extending from ihe box), and is iden- 
tified separately as an unusual . define the d fference 
between the 75ih percentile and the 25th percentile as the inter- 
quartile range (IQR) then unusual values are those values that 
are more than 1 5 IQR away from the appropriate percentile. In 
The case of an unusuaf value, Ihe whisker will end at the largest 
(smaJlest) data value that «s within 1.5 IQR of the 75th (25th) 
percentile. 

Reference 

1 PF VcllefnanarvdDC Hoagl-n Applications t &$ r 'pyfrng adaptors for,- 

Data Analysts Dujtbury Pre- : tgr 3 



* 6 
o 

m 5 

c 
o 

2 

t 3 
o 



Si 

E 



(a) 

70 - 
60 - 

i 

gso 

1 

o 40 - 

| 

■ 

S 20 - 

JD 

o 

10 -I 



-h 

10 



H 1 j 1 h 

20 30 40 50 60 

Characleristic of Interest 



Maximum Value 
(Extreme PointJ 



70 




75th Perce nttte 
50th Percentile 

25th Percentile 
Minimum Value 



m 



Fig. 1. (a) Histogram and (b) boxplot tor the same data 
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Circuit Design Using Statistical 
Data Analysis 

A methodology for setting limits of uncertainty on amplifier 
output power based on measured prototype data facilitates 
the design of amplifiers that minimize output power 
variations. The result is a printed circuit board requiring no 
adjustments. 

by Karen Kafadar and Lynn M. Plouse 



A COMMON SITUATION facing the circuit designer 
is the selection of design parameters for certain cir- 
cuits so that specified targets and variations in out- 
put power levels are met. This is the case for several circuit 
boards that are part of the HP 86792A Agile Upconverter, 
a fast-switching synthesizer that can upconvert complex 
modulation with up to 32-MHz bandwidth to anywhere in 
the 10-to-300n-MHz range. These boards contain several 
circuits such as uiupliHers, pads, frequency dividers, power 
splitters, switches, and filters. While some circuits, such 
as pow r er splitters and dividers, involve few design choices, 
other circuits, such as amplifiers and pads* involve param- 
eters that can be chosen in a variety of ways to nif^l target 
output levels, The aim is to select values for those param- 
eters so that variation uf the output power is minimized, 
[| the overall variation in I he output power is small enough, 
no adjustments are necessary, thereby reducing cost and 
labor requirements in production. 

Statistical methodology has been applied to this project 
in several ways. Measurements on the different types of 
circuits provide information on output levels and expected 
variation. Such data, collected in the prototype phase, is 
used to select design parameters for amplifiers so variation 
of the output power can be minimized. Variations caused 
by other sources, such as temperature and Input power 
variations, are assessed. Finally, the results from the lab 
pilot phase are compared with those from the prototype 
phase for prediction of performance in production. 

This paper discusses the application of statistical 
methodology to seven printed circuit boards used in the 
prototype design phase. This project illustrates the power 
of statistics and the value of a statistical analysis of data 
collected early in the design stage. The success of ihis 
approach is verified by the close prediction of lab pilot 
units from the prototype, and, more important, by the labor 
savings from designing a printed circuit board requiring 
no adjustments. 

Circuit Layout 

The new r synthesizer upcon verts a modulated input sig- 
nal. Several fast switches allow fast hopping of the carrier, 
Since each switch drives a mixer, its output level is critical. 
If the output power is too high, the mixer's spurious outputs 
(spurs) increase, degrading the overall spur performance 



of the instrument. On the other hand, if the switch output 
power is too low, the signal-to-noise ratio gets worse, Op- 
timally, the output levels of the switches should be within 
a 3-dB window, regardless of output frequency, tempera- 
ture, or components used to assemble the circuit, 

Often the target output levels on circuit boards of thr 
type are achieved by adjustments for each pole of the 
switch, Such a procedure requires at least 16 adjustments 
for these particular boards, depending upon how I hey were 
implemented. At a minimum rate of $12 per adjustment, 
this is highly costly in production. A better way uses statis- 
tical analysis to determine the loss and variation of the 
circuits in the switch modules, [n particular, some of these 
circuits (e,g,, the amplifiers that drive the circuits) can be 
designed to minimize output variation caused by variation 
of the incoming power. 

Sources of Variation 

Having defined the goals of the circuitry, data on various 
circuits must be collected and analyzed to help realize 
these goals. Depending upon availability of components 
and required test time, different circuits of each type are 
measured. This data provides estimates ol output power 
levels and variations as a function of input power. Several 
sources of variation in these measurements can be identi- 
fied: 

Variations in incoming power levels 
* Errors in input power level transmitted by the circuit 
■ Temperature biases on the output 
9 Circuit component variation 

to Variation in the manufacturing of the circuit board ma- 
terial 

Measurement error. 

In this study, the incoming power is actually the output 
power from an earlier portion of the instrument, which 
itself has a prespecified variation by design, The transmit 



Switch arid 
Low-Pass Filter 



Fig. 1. A typical power path on one board at a single fre- 
quency (470 MHz) 
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sion of the error in this incoming power depends upon the 
circuit (see next section). 

AU circuits are required to perform over a temperature 
range to pass environmental testing; hence the offsets in 
power levels because of extreme heat (70*C) and cold (5 a C) 
are estimated (see "Temperature Testing/" page 14). 

Having assembled all circuits on printed circuit boards, 
variations in lots of hoard material or components may 
affect output power; again, such variation is best estimated 
by taking measurements on different boards (see "Compo- 
nent Variation," page 15). 

Finally, this analysis involves two components at vari- 
ation that come from measurement error. The first arises 
directly from the test system and is assessed most typically 
from customer specifications. The second component in- 
volves an error imposed by the analysis, because input 
power levels are not identical for measurements of circuits 
taken on different occasions [different temperatures, com- 
ponents, etc.). To make fair comparisons of the output 
levels, most data had ta be linearly interpolated to a sum- 
mon input power scale. This common stale was an interval 
within the span of the data confined to the integers (e.g., 
- 10 dBm, -9 dBm, .... 10 dBnx). The maximum error in 
linear interpolation is given by one eighth of the second 
derivative of the function and estimated by one fourth of 
the second divided difference- Only rarely did this maxi- 
mum error exceed 0.02 dB. Generally, this error was a 
negligible component of the overall uncertainty in the 
power level. However, in other less precise measuring situ- 
ations, this error might not be negligible. 

The variations from these different factors are estimated 
and combined Into an overall standard deviation by a root 
sum of squares [RSS) method. This method is almost always 
justified by the fact that the variance (square of the standard 
deviation) of a sum of independent random quantities 
equals (he sum of thr variances, regardless of how (be 
individual random quantities are distributed. In some I 
where two factors were measured simultaneously (e.g,, 
temperature and boards), the interaction between them was 
negligible and was therefore ignored. 

For simplicity , and in accordance with recommendations 
by the Bureau International des Poids et Mesures, 1 all vari- 
ations are given in terms of plus or minus one standard 
deviation fSIJ or o], except where noted, When a Gaussian 
(Utopian) model for the data applies. 2-576 multiples of 
the population SD should include 99% of the population. 
When only ten units have been measured and the Gaussian 
model is assumed, :U)18 multiples of the estimated SD 
includes ^5% of the population with 90% confidence, 2 A 
rough translation of a customer Specification of a test sys- 
tem into a lo uncertainty is mostly guesswork; often one 
half to one third of the customer specification is treated as 
one standard deviation. Customer specif icatiou of the 
power meters in this study is 0.2 dB; hence tr Ml -- 0.07 
dB. (Throughout this paper, u will refer to an estimate oj 
the variation from a specific source, which will be iden- 
tified by the subscript on o\) 

Transmission of Error 

In this section, we illustrate the statistical methodology 
by considering one path from one board corresponding to 



one frequency input and output; this path is shown Ln Fig- 
1 . Power is transferred along this path through a series of 
its that includes a power splitter, an amplifier, a 
►me pads la absorb excess power. The input 
power level from the preceding circuitry for this path is 
estimated at 3.91 dBm ^0.94 dB. The power splitter incurs 
a loss of 3.8 dBm ±0.16 dB. Hence, using RSS. the power 
preceding the amplifier is at most 0*11 dBm i:0>95 dB, 
Increasing the loss of the amplifier input pad will lower 
the output power but not change the variation of the trans- 
mitted power. The goal is to design the amplifier so that 
i output power after the switch is -6*5 dBm. with 
a 99% uncertainty interval of ±1.5 dB, 

Knowledge of the transmission of error variation through 
a circuit provides the greatest leverage in reducing overall 
variation in the output power of the switch. Consider a 
typical response curve from a switch in Fig. 2, It is clear 
that the output variation transmitted around P { , is L 
than that around P v To a first approximation, this output 
variation depends upon the slope of the curve at the input. 
Taylor's theorem provides the basis fur this iaw oi prop- 
agation of error" : 3 

output = i[P] = f(P ) + (P P f] )f(PJ- 

= constant 4- [P - constant) x constant. 

Since 

SD (constant + X) - SD(X) 

SD (constant x X] - constant X SD(X) 

w r e have that 

SD(output) - |f'(P u )| x SDfinpin | ( 1 ) 




P ±066 dB 
Power into Switch fdBm) at 470 MHz 



P f *0.66 10 
dB 



Fig. 2* Typical switch response curve at 470 MHz, showing 
tnpul power atP a ±0 66 dB to achieve target output of -6 5 
± 0.5 dBm 
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Thus, output variation is related to the input variation 
through the slope oi the I unction relating the output lo the 
input near the point til input power. One may question the 
accuracy of this variation when we use only the leading 
term in the Taylor series expansion. Examples indicate thai 
such approximations are in fact quite satisfactory. 4 Further 
uses of this law of propagation of error will he demonstrated 
later in connection with combining several sources of vari- 
ation into a single uncertainty statement, 

To meet other specifications, the switch must operate 
where its response is linear. Using equation 1 8 and given 
the maximum allowable standard deviation of the output 
and the slope of the curve of the designed switch, one can 
calculate the maximum allowable variation in input power 
to the switch. To illustrate, the switch in Fig. 2 has a target 
output of -6.5 dBm and an allowable variation of ±1.5 
dB around this target. The allowable overall (RSS] standard 
deviation of the output power is therefore 0.6 dB, since 2.5 
multiples of this standard deviation equals the desired 
± 1.5-dB window with 99% confidence. If the output power 
variation from components is 0.3 dB, then the allowable 
output variation from transmission of error by the switch 
is Qi£ dB, since D.3 dB and 0,5 dB yield an RSS result of 
approximately 0*6 dB. The slope of the switch function at 
the input level (denoted P n in Fig, 2) corresponding to an 
output of — 6.5 dBm is 0.76. Thus, the maximum allowable 
variation in input power to the switch is 0.5/0.76 = 0,66 
dB around a target input power of P = 4.69 dBm. 

Fig, 3 shows a family of amplifier curves designed with 
three values of collector current. Notice that the effect of 
the collector current is to vary the point at which the slope 
of the curve starts to decrease, that is, the point of amplifier 
compression. The price lor this smaller slope (and hence 
smaller output variation] is less output power. Recall that 
sufficient power out of the amplifier is necessary to drive 
the switch and obtain the target — 6.5 dBm out of the switch. 
(It the amplifier output exceeds 4-69 dBm. the loss of the 
pad before the switch can be increased to absorb the excess 
power.) In addition, we seek those curves whose variation 
at the output of 4.69 dBm or above does not exceed 0.66 
dB. Fig, 3 shows that the slope, and hence the output vari- 




Allowable fnterval 
(-0.76 to 0.11 dBm) 



H 

-6-4-2 2 4 6 

Power into Amplifier (dBm) 

Fig, 3- Amplifier response curves for three values of collector 
current at 470 MHz. 



atiorc is smaller al higher inpul power levels [amplifier 
compression), but recall that, on average, only 0.11 dBm 
of power is available at the input to the amplifier. 

With a maximum limit of 0,11 dBm of input power, we 
calculate the slopes of the curves corresponding to the 
three amplifiers in Fig. 3. These slopes are estimated by 
successive first differences and are shown in Fig. 4. Notice 
that the slopes do not drop below the maximum allowable 
slope of 0.66 dB until the input power reaches at least 
-0.76 dBm, 1.7 dBm, or 3.7 dBm, respectively. Only the 
15-mA curve satisfies both the incoming power constraint 
and the maximum slope constraint. 

The allowable interval (-0,76 to 0.11 dBm] is indicated 
in Fig. 3, Since the slope is smallest at 0.11 dBm, one might 
reduce the output pow r er of 7.31 dBm after the amplifier 
with a 2.62-dB pad. For other engineering reasons, a 1-dB 
pad is placed before the amplifier. This increases the slope 
only slightly to 0.70. The output of the amplifier at P in — 
- 0,89 dBm is 6.62 dBm, so a 2-dB loss on the pad preceding 
the switch ensures the target - ft, 5 dBm a I the switch out- 
put. 

This same analysis is applied to all circuits in the pro- 
to type stage of the design. An additional consideration is 
the effect of temperature on the circuits. This effect is dis- 
cussed in the next section. 

Temperature Testing 

Circuits must perform under a variety of environmental 
conditions, One such set of conditions is extreme temper- 
atures. Measuring at least two circuits of each type at 5°C 
(cold), 25°C [ruom) t and 70*C (hot) allows an estimate of 
I he effect of tempera! ore on the circuits. This effect on 
some circuits is minimal; for example, for the power splitter 
in Fig, 1, the change in power is only ±0.05 dB T depending 
on whether the circuit is hot (-) or cold [* ). Likewise, 
amplifiers are largely unaffected by temperature. Fig. 5 
shows a set of six curves: two amplifiers, each measured 
a I room. hot. and cold temperatures. The curves are nearly 
identical. Heat decreases the output power by an average 
of only 0,2 B dB, and cold increases the power by an average 
of only 0,33 dB. (The temperature biases are nearly zero 
when the amplifiers are in compression*] Regardless of 
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Fig. 5, Six amplifier response curves measured at 5*C. 25"C, 
and 7(TC for two amplifiers 

output frequency, an overall 1*7 uncertainty of 0.27 dR 
accounts for variation in both temperature and amplifier 
components. 

The temperature effect is not negligible for switches. Fig. 
6 shows the difference in power between room and hot 
temperatures and between room and cold temperatures for 
two switch assemblies, Notice that this effect is not constant 
for oil switches. In the region of interest (P in — to 5 dBm), 
the power at 25 Q C is roughly 0.8 dB higher than at 70°C for 
niie switch and 0.96 dB higher for another switch* Simi- 
larly, one switch is 0,4 dB lower at 25 C than at 5 a C, while 
a second switch is 0,3 dB lower at 25°Q The temperature 
biases for the two switch assemblies are averaged: AHot = 
-0,88 dB and ACoid = 0,35 dB, The variation in these 
averages between different switch assemblies constitutes 
part of the component variation, which is discussed in the 
following section. 



Component Variation 

Component variation can be assessed in two ways: man- 
ufacturer specification or actual measurement data. The 
first approach requires a translation of the manufacturer 
specification into multiples of cr Usually the best guess is 
to treat the specification as a 3cr limit. In this study, this 
approach was limited mostly to characterizing the variation 
of purchased low-pass and bandpass filters. More often, at 
least two circuits of each type (power splitters, printed 
circuit board filters, amp libers, dividers, switches , XOR 
gates] were measured to provide more real is tit; determina- 
tions of the uncertainties in the components. Such uncer- 
tainty necessarily includes measurement error. 

We have seen the effects of component variation both in 
the amplifiers and in assessing the temperature biases on 
the switches. The component variation, like the tempera- 
ture effect, is small for the amplifiers: lrr estimate of the 
combined variation is 0.27 dB. The component variation 
in the switch circuits is higher, as Fig. 7 illustrates. Here, 
only a portion of the switch curve is plotted to exaggerate 
the differences in the curves resulting from measuring two 
different circuits on two different boards. The lower curves, 
representing measurements of two circuits on board A, are 
clearly closer to each other than they are to either of the 
upper curves, representing measurements of two circuits 
on board B« The board variation is substantial; here, the 
difference is approximately 0.7 dB. This variation may be 
caused by various other components on the boards , such 
as the low-pass filter that follows the switch. On a given 
board t measured differences between circuits are smaller. 

A difference of two measurements can lie translated into 
a Icr estimate by dividing the difference by V 2 . (In general, 
for Gaussian data, an estimate of u from the difference 
between the largest and smallest of n observations is ob- 
tained by dividing this difference by Vn, when n<10. See 
page 317 of reference 4.] This additional vinutjon is in- 
cluded in the overall uncertainty statement, as we discuss 
in the following section. 
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Fig. 6. The upper curves show 
room-temperature (25 li C) output 
power relative to hot-temperature 
(7tTC) output power as a function 
of input power for two switches, a 
and b The tower curves show 
room-temperature output power 
relative to cold-temperature (5*0) 
output power for the same two 
switches. 
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Power into Switch (dBm) 
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Fig. 7. Response curves for two switches on board A (tower 
curves} and two switches on board B (upper curves) 

Combining Uncertainties 

The root sum of squares [RSS) rule and the Taylor series 
rule are special cases of the general law of the propagation 
of error, which can be used in a general way to combine 
estimates of variation into a single overall uncertainly state- 
ment. The general law r of propagation of error says thai, for 
a function of k independent variables 



y = f{x t , x^ t 



x k \ 



the variance (square of the standard deviation) of y is ap- 
proximated by 

Var[y) = 2 Idf/dxfVail^) 

w T here the partial derivatives are evaluated at the sample 
means of the k variables [see page 363 of reference 5), 
Notice that when k - 1, taking I he square root of both 
sides leads to equation 1 above. In this paper, the response, 
output power* involves variations in input power, compo- 
nents, board version, and measurement error. Temperature 
effects are treated as additive biases, The simplest model 



for output power is an additive one: 

output power = Afinput power) + component effects 
4- board effects 4 ME 

where A( input power) is an amplifier function like those 
show T n in Figs. 3 and 5 and ME stands for measurement 
error. According to the law of propagation of error, 

Var(output) = A' 2 (P tl J x Var(input) + Var( components) + 
Var( boards] + Var(MK) 

so it follows that 

SD(P am ] = |A' z (P in )o-f n + <x; omp + <xL,rd t <v.\ 11 - 

which is the familiar root sum of squares formula. Notice 
thai nowhere have we assumed a specific: probability dis- 
tribution for any of these variables, only that the variables 
are independent. 

We are now ready to combine the previously identified 
sources of variation into a single uncertainty statement, 
applying this law after every circuit. We treat the power 
levels at the three different temperatures separately, since 
5°C and 70°C are extreme conditions, and in a given operat- 
ing environment, the temperature is not likely lo vary over 
such a broad range. Also, the measurement error <r M ,. : has 
been included In I he measurement of component and board 
variation. 

Consider our previous illustration at room temperature, 
for which P in = 3.91 dBm ±0.94 dB. The power splitter 
and the pad incur pow r er losses of - 3.8 dB ±0.16 dB and 
1 dB ±0.00 dB respectively. Thus P in to the amplifier is 
-0.89 dBm ±0.95 dB. The amplifier increases this power 
to 6.62 dBm, where the slope is Q.7 t so the output variation 
is reduced to ±0.70 dB, Amplifier component and temper- 
ature variation of 0.27 dB raises this uncertainty to ±0.75 
dB. A 2-dB loss results in P ir] = 4.62 dBm ±0.75 dB to the 
switch* which itself has a slope of 0,76 at P in . Thus the 
output power is -6.48 dBm ±0,57 dB. Combining this 
with the board variation (±0.40 dB) by the root sum of 
squares method, we have a final uncertainty of ±0.70 dB. 
A similar path of analysis follows from the power levels 
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Fig. 8. Triangles indicate the ex- 
pected values and the 2 Su itmits 
of the power output from each 
stage for three temperatures: 
room (right, color), hot (left, color), 
and cold (left black). 
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Fig, 9. Room. hot. and cold response curves for an average 
of four prototype synthesizer boards , and a room-temperature 
response curve for a lab pilot board 

at hot and cold temperatures, 

A graphical display of all these numbers is shown in Fig, 
8. The expected power level at each stage is indicated by 
the center of the triangle on the vertical axis corresponding 
to each circuit. There is one triangle for each temperature: 
room on the right h hot and cold on the left. The length of 
the base of Hie triangle extends to 2,5 multiples of the 
standard deviation of the power output at each stage, Thus, 
considerable information is conveyed in this display: 
changes in power levels across circuits, variations in power 
levels at different stages, and differences in power levels 
at extreme temperatures. The reduction in variation from 
the beginning to the end of the path can be seen, and tin- 
target P out ea - G.5 dBm is attained within the desired 
limits at room temperature. 

Comparing Prototype Analysis with Lab Pilot Results 

The path of power shown in Fig, 1 can be formulated 
into a single functional relationship between input power 
to the power splitter i-inci output power from the switch. 
This functional relationship is computed using linear inter- 
polation in the data collected on amplifiers and switches 



and is shown in Fig. 9. A lab pilot board was measured for 
its power transfer between the initial and final points of 
the circuitry; it is shown in Fig, 9 for room temperature 
only. The agreement in the general shape of the curves 
between observed lab pilot measurements and exp»- 
values based on the prototype is excellent. Any absolute 
difference between the curves can be explained by imped- 
ance mismatches; pads can be adjusted accordingly- Once 
these pad values are ascertained for the lab pilot, target 
power levels are achieved within the designed limits, with- 
out the need For hitter adjustments of pads in production. 

More Complicated Analyses 

Not all power paths in the instrument are as straightfor- 
ward as the one illustrated in this paper. Fig, 10 shows a 
display similar to Fig. 8 of the uncertainties in a power 
transfer involving incoming power at one of two different 
frequencies. While the approach to amplifier design and 
estimation of component variation is similar, the estima- 
tion of variation resulting from input power is complicated 
by the presence of incoming power at another frequency- 
Thus, measurements must be taken by varying power levels 
at both frequencies. Other factors, such as board type and 
circuit components, can also be varied, and a statistically 
designed plan to take measurements on several different 
factors simulaneouslv is encouraged, 5 
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Statistical Calibration of a Vector 
Demodulator 

Circuit performance characteristics must he determined 
and accounted for in instrument calibration. This paper 
describes an algorithm for statistical demodulator 
calibration to guarantee high precision in the demodulated 

signal. 

by Karen Kafadar 



VECTOR MODULATION INSTRUMENTS such as 
the HP 8780A Vector Signal Generator 1 and the HP 
fly BOA Vector Analyzer" provide a test system for 
measuremenl applications in such fields as digital micro- 
wave radio, communications, and radar. Some of the 
technologically new developments in this series of prod- 
ucts are the modulation capabilities (complex, wideband 
signals), the high-resolution CRT display, and the 350-MHz 
bandwidth. 

The precision measurement capabilities offered in these 
products raise the issue of the internal calibration of (he 
instruments, All instruments need to be calibrated to ac- 
counl lor differences in components and test conditions. 
The accuracy and precision of these internal self-adjust- 
ments are crucial, since they have direct effects on the 
characterization of the device under test. The hardware 
measurements made for purposes of internal calibration 
must be combined in a logical and intelligent fashion to 
yield the correct calibration factors, 

Statistical methods play an important role In the calibra- 
tion routine. This paper describes a statistical approach to 
the calibration of the HP 8981 A Vector Modulation Ana- 
lyzer. The next section briefly introduces the use of statis- 
tical methods in such situations. The measurement process 
fur a vector demodulator is Hum described, using the HP 
H981A as an illustration. Two sections provide the techni- 
cal details of the calibration algorithm both with and with- 
out an accurate vector signal generator (e.g.. the HP 8780 A) 
as a reference. Finally, examples of real data illustrate the 
computations in the statistical algorithm and formulas for 
the uncertainties in the calibration factors. It will be shown 
I ha I the algorithm can be used to calibrate not only die de- 
modulator in the HP 8981 A but also external demodulators. 

Statistics in Calibration Routines 

Calibration consists of comparing a set of measurements 
from an uncharacterized instrument (e.g., an HP S081A] 
with a defined reference standard (e.g,» an HP 8780A] ac- 
cording to a measurement algorithm. Thus the calibration 
model relates the observed measurement readings to the 
reference standard. A model is never perfect, and the dif- 
ference between the model and reality can be character! zee i 
according to two types of errors that may be present: 

Part ot thrs paper appeared in MSN & Communtcalforts Technology. March 198S Repf inled 
by permission ol the publisher Copyright EW CownunicaTions Inc 



Systematic errors resulting from incomplete specifica- 
tion ol the calibration model 

Random errors, or small, unpredictable fluctuations that 
affect every measurement but are themselves unmeasure- 
ahle. 

Statistical methods can be used to address both types of 
errors, First, the calibration should he designed to allow 
identification of possible systematic departures between 
the model for measurement and the observed data, Second, 
the influence of random errors can he assessed when es- 
timating the parameters in the model that relate tu the 
measurement process and their uncertainties. In this paper, 
we shall assume that the random errors are independent 
and identically distributed according to a symmetric distri- 
bution; that is, that the individual errors cannot be pre- 
dicted in either size or direction, and that the chances of 
an erroneous measuremenl being loo large or too small art 1 
roughly equal. (Diagnostic tools for checking the validity 
of these assumptions and the consequences of their viola- 
tion are discussed in connection with the examples later 
in this paper,) 

Rarely is a single estimate of a target quantity sufficient. 
For example, reporting a sample mean without its standard 
error provides no information on the reliability of the data 
that went into that sample mean. The same is true for 
estimates of the parameters describing the measurement 
process, These estimates are the calibration factors, and 
certain limits of fluctuation may be desirable (e«g M gain 
adjustment accurate within 1%), Therefore, an important 
part of the statistics in a calibration algorithm is the deri- 
vation of associated measures of uncertainty for the calibra- 
tion factors. 

Vector Analyzer Measurement Process 

\n the HP 8981 A, the phase- modulated signal can be 
expressed mathematically in terms of its frequency compo- 
nents as 



s(t] = £A k (t]cos[w k t + </> k ) 

= |I k (t)cos[* k t) + Q k (t)sinKt) 



(1) 



Equation 1 expresses the signal in terms of its demod- 
ulated components I k {t] (in-phase) and QtJt) (quadrature). 
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Components that vary ovbi time can be observed by examin- 
ing I k {t) and QJt) over time at a particular frequency (u k . 
Ideally, at a fixed time, a stogie-frequency signal with con- 
stant amplitude is demodulated at 

I = cos(a) 

Q - t:os(a - (90° - <£)) 

where a is the angle of the vector signal and 6, the error 
in the quadrature phase of the demodulation, is zero. 

Since the HP 8981 A has the ability to display signals in 
either lime or vector form, amplitude and phase informa- 
tion must be recovered accurately. This requires special 
software functions to estimate and subsequently account 
for any amplitude or phase offsets. Proper calibration of 
the HP 8981 A thus allows simple measurements of 
amplitude and phase modulated signals. 

Fig. 1 is .i Mm k diagram for a simplified demodulator, 
such as the HP 8981 A, which decomposes the modulation 
on a signal into its in -phase and quadrature components. 
Accuracy and precision in this decomposition are essential 
to assess degradation of the signal's amplitude, phase; or 
frequency, or to assess cross-talk between the two channels. 
Since I and Q signals are used often in measurement sys- 
tems, we will discuss the calibration scheme for the Q-vs-I 
vector representation. 

Six sources for adjustments in a vector demodulator cao 
be identified, based on the block diagram of the instrument. 
These six sources correspond to stages in i he transforma- 
tion of the idea/ I/Q signal Into ihe m ensured, or actual ( 
I/Q signal. The calibration process uses statistics to estimate 
the sizes and directions of the adjustments attributable to 
these sources and then corrects the output so that the de- 
modnlaled I Q signal agrees with the ideal I Q signal. 

Let (I,QJ denote the pair of points representing the 1 and 
(Jchiiiinelsal B particular point in time at a given frequency, 
The HP H9W1 A hardware transforms the Coordinates [I,Q) 
itltO measured coordinates (x.y) through the following 

series of operations: 



[[,Q) -» Rotation - 
dc Offsets —* (x.y) 



Axis Sheer -* Compression 



Rotation can be viewed as a phase change in the ideal 
ignal, which can be introduced by any of the components 
before the demodulation, The axis sheer creates two addi- 
tional adjustments: quadrature adjustment, resulting from 
inaccuracy in tln> 90 split, and gain imbalance, resulting 
from different power levels in amplifiers or unequal power 
losses in other components after the split. The LQ measure- 
ments are rescaled by amplifier coin press ion. rind tic offsets 
may be added to the measurements < r any stage. 
Notationally, let us represent these adjustments as: 

7 = gain imbalance [1-channel/Q-channelJ 

(f> - angular difference between ideal quadrature (90°) 
and actual quadrature 

8 - angle of rotation (sometimes called lock error] 
p - scaling resulting from compression 

l r( sa dc offset io I-channel 
t^ = dc: offset in Q-ohunnel 



The four steps in transforming (KQJ to |x,y) can be ex- 
pressed mathematically as follows: 

Rotation through angle fc 



R: 



cos0 -sinB 

sinfl 



This transforms (hQI into (icostf - Qsinft Isinl •¥ Qeos0), 
Sheer transformation (quadrature adjustment #. gain im- 
balance 7): 



y 

sin<f» cos4> 



This causes two additional adjustments; the ratio in the 
two channels (I/Q) is now y, not I. and the tangent of 
the angle between the two axes i.LQ) is now [90" - &), 
not 90 p , 

The measurements are uniformly compressed according 
i i ,i factor p. 

The measurements are offset by a fixed amount 1 1,, in 
the I-channel. Q^ in the Q-channel). 
Putting these four effects together, the resultant point 
[x T y] can be expressed as a function of the ideal input (IjQ) 



l.:l 



SlUth oos<£ 



COS 

sinfl 



B -sinfl IT I " 
6 cosfl [ Q J 






- 



[i;J 



[2] 



Tho unobserved variables t y .e v are assumed to be uncor- 
rected random errors in the measurement of the x and y 
coordinates, That is, 

x = I + yiitcuzO}] -y[psmQ]Q + e v 
= ot () -i- citjl + ot^Q + £ x 

V - Qii + (pcosftsiiu/) • />sinflcos<£)I 

f [ -psin6fsin<£ + pcostfcosc£|Q +■ e v 
= ft, + 0,1 4 &Q t £ v 
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fig. 1 . Simplified vector demodulator. 
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where: 



a 2 - -*ypsinfl 

A - Qo 

/3 T — pcasfkintf) + psin#cos</j 

= [a]5in<£ — a 2 cQSff>)?y 
Pz = - ps i n fls i n *£ + pco s BdQ a $ 

The errors of interest can be recovered via the equations; 

Iq = ct f I- channel offset) 
On = A j (Q-channel offset) 
g**lf£+ /S) 1 ' 2 (Compression) 

6 = tan'^-ajs/aj [Rotation] (3) 

7 = M + 9S#{jjS t ^l)] 1 "^ ffitj gain imbalance) 

=* tan 1 [[ot 1 ^ 1 + aaPaJ/fa^^-oty/S,)] [Quadrature error) 

Notice that, apart from rotation, the ideal I and Qcan be 
recovered via the equations 

I = (X . - loJ/^p 

Q = [7(y - Qu) - sin f [x . - I E] J|/7p cos <fi t 

Estimation with Accurate Reference 

At this point, we need to distinguish between the case 
where ideal (I,Q) values are provided and the situation 
where only magnitude information about the (I T Q) points 
is known, Ideal (LQ) measurements can he obtained using 
a precise signal generator such as the HP 8 780 A, so that 
signals having several accurately defined phase modes can 
be guaranteed. Since this situation lends itself to the easier 
algorithm, it will be discussed first. 

It is evident that an input signal of constant frequency 
and magnitude and variable phase, namely 

s(t) = AcosM + <h). t = 0\l t . . , 

results in an I/Q display for which 

I 3 + Q 2 = (Magnitude) 2 , 

w T hich is a circle. It can be verified that the transformed 
points [x t y) lie on an ellipse. The parameters of this ellipse 
are related to the adjustments in the instrument as outlined 
above. 

Since we have six adjustments [parameters), we need a 
minimum of six phase states to obtain an exact fit for the 
parameters In equation 2, 

The 8PSK mode of the HP 8780A allows rapid generation 
of an I/Q signal at eight evenly distributed phase states, 
allowing two extra phase states for model confirmation. 
The coordinates (I,Q) of the points in any one of these 
states will be considered ideal, transformed (by the six 
adjustments) into measured coordinates [x.y]> For simplic- 
ity, we label the eight ideal states sequentially at those 



points having phases k(45*) t k = 0, 1, ..., 7 [Fig. 2). Thus, 
State I Q 






1 





1 


f 


f 


2 





1 


3 


-f 


f 


4 


-1 





5 


-f 


-f 


6 


II 


-1 


7 


1 


-f 



where f = 1/VY = 0.7071068, [Numerical association 
with each state is arbitrary.) 

Suppose n/8 points are measured in each state, so n mea- 
surements in all are taken. (This requirement is not essen- 
tial, but the equations are greatly simplified by encouraging 
orthogonality into the problem.] Denote the coordinates ol 
these measurements by (x lT y 1 ), . , . , (x n ,y n ). Associate the 
point with its ideal (I,Q) pair as described above: (Ii r Qih 
. , . , (I^QjJ. Then we need to fit « = [a^a^a^]' and fi = 
[pQifii, fi 2 Y in the two relationships; 

Xj = a + o^Ij + t* 2 Qi + e X[ 

y 3 = ft + |il, - #A + ^ 

The estimates of [ct^a^a^], namely [d n ,a T ,6^), can be 
obtained by a least sqoares algorithm. If we let M be the 
matrix of n rows and 3 columns corresponding to the n 
observations and ihree "carriers" (constant term, I, Q). that is, 



t.o -r 
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Q» 
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1.0 



Fig. 2, in-phase (I) and quadrature (Qj modulation states for 
8PSK vector modulation 
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M 



i> Qi 



Q„ J 



and calculate 



-VI 



</(n - 3) 



then it is well known 3 that 



=VT 



- 3). 



q = (M'M)''M'x 



MM M'j 



w 



are the least squares estimates of ot and /?, where x = 

(x lP x 2f . . . , x n Y and y — (y lt y 2 y n )' are column 

vectors and the prime sign (') indicates transpose. Since 
there are n/S measurements in each state, and the states 
are evenly spaced around the circle, the matrix M is ortho- 
gonal, and so 





n 







M'M = 





n/2 













n/2 



M'x =[1^. jj*& | t xA 

and thus 

<i - (M'M)~ l M'x = [ave(x), 2ave (x -I), 2ave(xQ)] r 

= (M'M) _1 M'y = [avefy], 2ave [yl), 2ave(yQj]'. 

The estimates of the target paraoieters (I , Qq> p, Q t y> <j>) 
can then be unscrambled using equations in [3] above. 



The divisor is fn — 3) instead of n because three parameters 
have been estimated in each channel, If the departures from 
the model have a Gaussian distribution, we would exped 
most of the e Si to fall within 2s x ur 3s x . and likewise for 
e y ♦ This criterion may serve as a simple check for unusual 
measurements (but see "A Caution/' page 24], 

Using the root mean squares of the model departures 
(residuals), we can calculate the standard errors of the pa- 
rameters in the model. Again assuming these residuals are 
Gaussian distributed, the standard errors of (& n , S t , &J are 
given by 



Thus, 



s x V diagonal of [M F M] 

SE[a ) = SE[ij = s x /V n 
SE(/3J = SBf&] - VV^ 



The standard errors for the other target parameters are a 
little more complicated, because they are functions of ot v 
a 2 , /3 lt and /3 2 * We rely on propagation of error formulas* 1 
for approximations to the standard errors of functions of 
variables. The approximations are simplified if 

2s y 2 /(np 2 } « 1, 

a condition that will often be satisfied if the fit is good [s y 
is small) and the number of points is moderately large. 
Then: 



Standard Errors 

We need some measure of uncertainty in our estimates 
of the parameters. For example, after calculating an esti- 
mate of y, say 7, how far off might that estimate be from 
the true value? One such measure is the standard error 
[SE) of the estimate. An approximate 95% confidence inter- 
val for the estimate is given by 

estimate ±2(SE), 

To calculate the standard errors for the six parameters, 
we first need a measure of goodness of fit in our model. 
This can be calculated most simply as a "root mean square* 1 
of the differences between the observed measurement x s 
or y, and its fitted value predicted by the model. Denote 
these differences, or residuals, by: 

e*> = *i - «o - <Mi " «*2Qi 
e yfc = y, - $0 - 0iU ~ 4 2 Qi 
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Fig. 3. Data for the first example (see text). Stx measurements 
were taken in each of eight states. 
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SE(/5) ~ s y V2/n 

SE(-y) ^ s x V / 27n 

SE(tan0J - s^V^Tn/c*; 

SE(tan0) ■* [(sj| + &?)[! + 1^)^/^^ 

where tantf and tanc£ are the estimates of tan// ami Iuik/j. 
So approximate 95% confidence intervals for the estimates 
of the errors are: 

7 ± mm) 

to ± 2SE(IJ 

0o±2SH[Q n l 

p_± 2SE(p) 

tan '{tanfl ± 2SE[tan0]| 

tan _1 [tan> ± 2SE(tan>]|. 

This calibration algorithm thus proceeds as follows: 
Measure (x it yj) for i = 1, <>., n 

Associate stale (I-jQJ for i - 1 n 

Calculate: ct — [cit Qp d 1( 6t 2 ) r 

Calculate target parameters: y,(&,pJhlu,Qn 
Calculate: goodness of fit (s x ,s v ) and residuals [e x .,e v ) 
Identify unusual residuals 
Calculate standard errors. 

Inaccurately Known Phase 

We now consider the situation where phase cannot be 
associated precisely with each measurement, Generally, 
amplitude information is more reliable, I hat is, 

I 2 + Q 2 = [I,Q]'(I,Q) = 1. 

Two possible signals I hat satisfy this requirement axe: 
A signal with continuously varying phase (e.g., two 
tones, offset in frequency] 

A signal with eight discrete phases [e.g., an fiPSK mod- 
ulated signal and its coherent carrier). 
Inverting equation 2, 



I 
Q 



Id lc i 



p X R J S 



x - l " e * 

y - On - e y 



and thus 

{I.qj'(I.Q) = (5) 

because (R'R) = (R'YfR" 1 } is the identity matrix [R is 
orthogonal). This can be simplified by noting that 



er¥firV lycoB^y 



* (7cos<£]~ 



cos$ -sinob 



cos<£ 
- sin<£ y J 



Let us assume that the errors in measuring x and y have 
zero means, are uncorrelated with one another, and have 
common standard deviation <t. Then, on average, 

ave(€ x ) = ave{€ y ) = 
ave(ejj - ave(€^] = a 1 

and so equation 5 becomes: 

[x-I„P - 2 7 (sin<W(x-I n )(y-Q JJ ] + 7% "Co) 2 + 



[y 2 + 1 J<r - y 2 p 2 cm 2 d> - 0. 



(6) 



If measurement error is negligible compared to the mag- 
nitude ol the adjust merits, the fourth term in equation 6 
can be ignored. 

Equation 6 is nonlinear in the parameters, of which there 
are now only five (clearly I lie rotation, or lock error, cannot 
lie estimated without relative phase information). Al- 
gorithms for nonlinear least squares often rely on gradient 
methods or on linearization of the problem via a Taylor 
series expansion. In addition, several methods for approx- 
imate confidence intervals far the parameters have been 
proposed. 3 These two issues are described in the Appendix 
on page 24. 
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Fig, 4. Histograms of the x and y residuals (e x , e y ) for the 
ftrst example The histograms show the differences between 
the transformed measurements in each channel and their 
ideal values as shown graphically in Fig. 2 
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Examples 

This first example serves as a "control." Six [x T y) mea- 
surements were taken in each of eight states (Fig* 3). All 
calibration factors are essentially at their nominal values. 

(a ,d lf 4 3 J= (0,000054. 0.16564, -0.068486) 
( O . l . 2 ) = { - 0.002694. 0.068058. 0. 163904) 
L B = 0.00005 volts (I-chaiinel of i 
Oo = - 0.00269 volts (Q-channel offset) 
p= 0,17747 volts (Compression) 
0= tan" 1 (0.39206) = 21,41* (Lock en 
y= 1.0100(1 Q gain imbalance | 
<£ = tan" 1 (0.00152) = 0.087 1 D (Quadrature 
error), 

Root mean squares; s* - 0.0012162;s y = 0,0008686 volts. 

Approximate 95% confidence intervals (estimate ±2SE): 

I :( - 0,00044, 0.00055) volts 
Qo:( - 0.00295 T - 0.00243) volts 

7;[1.0053* 1,0147] 

p:[0,l 7711. 0.17783) volts 
Um<f>:[ - 0.001 44, 0.00448) -* <j> : [ - 0.08°, 0.26°) 
ianft(0.39153. 0.39253)^> 8: (21.38°, 2l.43 c ). 

Fig. 4 shows the histograms of the x and y residuals. 
Clearly there is one aberrant measurement. This measure- 
ment occurred in state 4(1 = -l t Q = 0), A comparison 
of these estimates with those obtained with a more robust 
procedure ( dow ewe igh ting this point) verified I ha I (he 
point is not so extreme that it affected the quality of the 
overall fit. 

Using the algorithm for unknown phase (see above] , 

1^:0.00019 (-0.00044, 0.00055) volts 
Cfe; - 0,00243 ( - 0,00295, - 0.00243) volts 
-y:1.0100 (1.0053,1 0147) 
fr.i)A 774^ [0. 17711,0,1 7783] volts 
tan(^]:f].00230 [ - 0.00318, 0.00622) — <fr 0.13° 
(-0.18°. 0.36°). 



Notice that these solutions fall well within die confidence 
intervals of the linear (known-phase) method. 

The second example involved a large quadrature adjust- 
ment. Hi^ht !xa 1 measurements were taken in each of eight 
states (Fig. 5). 

[&DAi*ft 2 ) ,B (0-089621.0.143384. -0.292603) 
| n . ft, ft) = (0.01 1 144. 0,227279, 0.244292) 
tu = 0,089621 Volts (I-channel offset) 

Q,, - 0.01 1 144 Volts (Q-ehannel offset | 
p= 0.33367 Volts (Compression) 
0- tan" 1 (1.11516) = 48.12= (Lock error) 
y= 0.97656 (L'Q gain ratio] 
<$ = tan" x ( - 0.36852] = - 20.09* [Quadrature 
error), 

Root mean squares: s x = 0,0019726; s y = 0.0015534 volts. 

Approximate 95% confidence Intervals (estimate ±2SE): 



I tl : (0.08913, 0.09011) volts 
Q Cf :(0.01076. 0.01153) volts 
7:[0.97586, 0.97726} 
p;[0.33312, 0.33422) volts 
tan<£;(- 0.36263. -0.36891)-*: 



20.25°, -19,94*) 



LaiiW.fl. 10963.1.12069)^ :[47.97 c , 48. 26°) 
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Fig. 5. Data for the second example. Eight measurements 
were taken in each of eight states. 
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Fig, 6. Histograms of the x and y residuals (e^ r e y ) for the 
second example 
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Appendix 

Solving Nonlinear Least Squares Problems 



The accompanying paper discusses Ihe 50=11 [ron To estimating the parameters in a 
modeE ihal Is linear in Ihose parameters, e.g., 



a, J, 4 (i 2 Q, + ' 



J = 1 n. 



from which it follows that 



= % - ftM tl y l i%)/f , (jt it y l !fl). 



-f - Wftf - Go) 



tiW-(W'W) 'W1, 



(3) 



Equation 3 >s iterated until satisfactory convergence is obtained, (For the examples 
in the accompanying paper, three to four iterations assured convergence [0 Ihe fifth 
significant figure,) This method of finding (he least squares solution to a nonlinear 
problem is known as a Gauss -Seidei iteration and has nearly quadratic convergence 
when the data is known to fit the model well. The algorithm and other methods are 
described in reference 6 listed on page 25. 

The coverage of approximate confidence intervals for parameters sofved by nonlinear 
least squares algorithms such as this one has been investigated (see reference 5 listed 
on page 25) The simplest and easiest approximation to the standard error of *K fc 
gtven by 

<i|klh diagonal element of (WW) -1 ] 1 '* 

where it is the standard devialion of the residuals. 

However, an interval of ' 2 of these "standard errors' around the estimated param- 
eters may have coverage as low as 75% instead of the nominal 95% Hence, tor 

conservative purposes., we recommend using the multiplier 3 for an approximate 95 D k 
confidence interval Further details can he found in references 5 and 6 listed on page 25. 



The estimates are given expltcitfy m equation 4 of the paper. 

When the model is nonlinear in the parameters, such as in equation 5 of the paper, 
one approach is to approximate the lunotion by a first-order Taylor series expansion 
and solve as for the linear problem. For a function, say f. of the data values (x,,yj 0! 
only one parameter say a, 

= lfa.y>;#) - t, 

- F{x,y.:tf - \tt rt^d-ffx.. y,.w).d^ + e, 



(1) 



Let Bq be a prior guess for the parameter if, Since there are n values of f(*,,y,;<y and 
' (K„y -"rj. we can caf curate a new n by replacing ihe quantities in equation 1 by their 
averages. 

Equation S of ihe paper actually has five parameters; U |-y,p.l<].Qij.<W. The multi- 
variate extension of Taylor's theorem becomes 

fftfflrt) - tai,Vi:ft) <■ [tfU^.^K'W %). 

where [tff(X|.y, ;^,)^^J Is a matrix of partial derivatives of the function f with respect to 
each of the five parameters IfyJ. Let t, denote the column vector of n values of the 
function f(x |k y,:et), i = 1, .... n and let W be the n.'-5 matnn ol vaJues of ihe function 

^\x, r %f^)^t\ r 1 = 1 n. k - 1 , , .., 5, all evaluated at the current guess IV Then, ideally, 

f{x,.yi;0) B 0. 30 



f2j 



This problem is in the same format as that for the linear least squares problem . where 
W was the t\x3 matrix M of rows consisting of 1s, I, . and Q,. and was the vector of 
parameters denoted by (k. 

Like equation 4 of the accompanying paper, the solution to equation 2 is given by 

{« - AM = -(WW) 'W'f d 



Fig. 6 shows the histograms of the x and y residuals. 
Using the algorithm for unknown phase (see above]. 



I o :0. 088 71 volts 

(i,:U.01 187 volts 

7:0.97723 

p;0. 3 39 7 8 volts 

tan M\ -0.36460-*: 



20 031 M.i?<, J].:i6°). 



Notice that a large quadrature adjustment [ - 19. 25*) and a 
slight adjustment in the gain imbalance are indicated. The 
other errors are considered negligible. 

A Caution 

If ihe model departures [residuals) are not Gaussian dis- 
tri billed for any reason (e.g, f outlying points as in the first 
example, or misiden titled stales), or if there are unequal 
numbers ol measurements in each stale in ihe known-phase 
algorithm, the estimates of the target parameters and their 
standard errors may be seriously biased. Least squares fit- 
ling procedures have a tendency to make all residuals about 
I lie same magnitude, even if. in fact, the residuals in all 
measurements were tiny, except for one. One indication 
of this case is if the median of the residuals is not zero, 
even though their average is, Bui I his criterion would be 
ineffective for detecting two residuals in opposite direc- 
tions. To reduce ihe effect of such aberrant measurements, 
the algorithm is made more robust by assigning a weight 
to each measurement based on the magnitude of its depar- 
ture from the model. This feature is essential to avoid seri- 
ous bias in the calibration factors. 67 

Summary 

A vector demodulator can be calibrated using a statisti- 
cally-based algorithm to determine the size and direction 
of the adjustments necessary to guarantee accuracy of the 
output. The uncertainties in these adjustments can be de- 
termined based on some assumptions about the measure- 
ment process: 

The ideal (I.QJ measurements of the signal in the two 
channels have been transformed to (x,y) measurements 
according to the following series of operations (listed in 
order of occurrence): rotation (lock error), sheer (quadra- 
ture, gain imbalance), compression* dc offsets. 
Errors in the measurements in the two channels (x,y) are 
uncorrected. 

When the relative phases among ihe points are known, 
data is collected in phase-stepping mode for eight slates, 
and the state (0, ...j 7] can be associated with each mea- 
surement. It is convenient if the number of measure- 
ments is the same in each of the eight states. 
When phase information is known only imprecisely, 
measurements are made that ideally lie on a circle in 
ihe J-Q plane. Input signals must have constant 
amplitude and varying phase (e.g, T continuously, as in 
a two-tone setup, or discretely, as for 8PSK). 
Based on these assumptions, a straightforward linear 
least squares algorithm can be used to determine the sizes 
and directions of the adjustments when relative phases of 
the incoming signal (such as 8PSK from the HP 8 780 A) are 
known. If this is not the case, then a linearization of the 
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nonlinear least squares problem can be solved, with some 
Loss in precision of the adjustments and increased compu- 
tational time. Since this approach can be used more gener- 
ics calibrate external demodulators as well, it is the 
algorithm currently programmed in the HP 8981 A, In both 
[mate 95% confidence intervals for the adjust- 
ments can be derived. 

tuples with data collected on the HP 8981 A using 
the HP 8780A as a reference signal illustrate the implemen- 
tation of the algorithms. The examples illustrate the impor- 
tance of incorporating an efficient and robust calibration 
algorithm using A principles to control the 

hardware. Proper modeling of the system combined with 
the appropriate methodology results in accurate calibration 
factors as well as an assessment of the uncertainty in these 
factors. 
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An Availability and Reliability Calculation 
Tool for Computer Systems 

This software program helps R&D engineers make trade- 
offs in designing for reliability. It also provides marketing 
representatives with a tool to evaluate the reliability and 
availability of customer-defined systems. 

by Wutf D, Rehder 



AS OFFICE NETWORKS, computer clusters, multi- 
processors, and other ensembles of hardware de- 
i vices are integrated into large systems, the tradi- 
tional problems of reliability, availability- and serviceable 
ity/supportability (the RAS features] gain a new and in- 
creased importance. The advent of complex systems that 
provide a variety of solutions rather than just computations 
has shifted the computer dependability issue from compo- 
nent reliability requirements to the customer-oriented con- 
cept of system availability^ 

A new algorithm for modeling the RAS features of large 
systems has been implemented in a prototype software tool 
for internal use within Hewlett-Packard. Before illustrating 
the benefits of the RAS tool, it is important to clarify the 
distinction between device reliability and. system availabil- 
ity and serviceability. 

Fig. 1 shows how a system might degrade I h rough several 
stages during which useful work can still be done. For 
instance, a printer or workstation maybe down temporari- 
ly^ or may not be functioning optimally. In a situation like 
this the system is to some degree imperfect , but still usable, 
hi most cases, a system with such operational faults is si ill 
available for proper, albeit reduced, service. 



Designing for this possibility of controlled degradation 
and high availability in the presence of faults requires the 
inclusion of means to handle fault diagnosis, system recon- 
figuration, and operational recovery. Preferably, this is all 
done intrinsically by a local control or support processor 
or remotely from a central support facility. 

Not every fault results in a failure, even though the sys- 
tem as a whole may contain a transient or permanent error. 
There exist many examples of such fan il -tolerant system 
behavior. A well-known example of a transient error is an 
alp ha- particle strike of a memory cell causing a one-bit 
error, [f the system is capable of one-bit error correction, 
then system operation is not affected at all. As another 
example 1 deconfiguring certain defective parts, from mem- 
ory and processor boards to peripherals, does result in a 
degraded system, but not necessarily in system failure. 
Similarly, a defect (bug) in a soil ware program need not 
corrupt the system operation. 

On the other hand, if the effect of an error escalates 
sufficiently to prevent the continuation of acceptable ser- 
vice for the user, then the system state falls below the 
threshold of availability, and a penmment failure occurs 
that requires external interference and repair. 



High 
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E 

1 



1 
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Time 



Fig. 1 « States of a degrading sys- 
tem. 
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RAS Calculations and Sensitivity Analysis 

For a specific system configuration, the RAS tool com- 
putes the availability A. mean time between failures MTBF. 
annualized failure rate AFR. mean time to repair M'lTK, 
and expected total up and down time per year, given the 
failure and repair rates of the individual devices. Further- 
more, the uik ^ added and the ratio AFR/KS 1 
is computed, where KSVAS = thousands of dollars of \ ■ . 
added shipments, 

To make data entry easy and the output format clear and 
comprehensive, segmented windows on the screen prompt 
;ser for input values and the soft keys are activated to 
browse back through the data and edit it if needed. Results 
of (he calculation can he inspected on the screen or can be 
printed in the form shown in Fig. 2, Here an abbreviation 
such as 2(1) means that out of n = 2 devices present only 
k = 1 is needed for the system to be operational. 

Besides providing the kind of snapshot of the major fig- 
ures of system merit shown in Fig. 2 P the RAS tool's flexible 
capabilities can be employed for the purpose of a "What 
fl£ M * J sensitivity analysis. It will answer questions about 
the impact of additional redundancy not only on system 
availability, but also on its complexity, its cost, and the 
ratio AFR/KSVAS- Likewise, some simple trade-off deci- 
sions about where to employ redundancy and how to in- 
crease service can be made. For instance, will a 2-out-of-3 
[triple modular) redundant power system design more ef- 
fectively increase RAS than, say T selective dual porting? 



There are many complex technical issues and questions 
of resources, market requirements, and competitive pres- 
sures that determine the outcome of many trade-offs, prob- 
ably even more so than the purely numerical suggestions 
the RAS tool is able to give, Bui for an initial assessment 
of the relative impact of different configurations on the 
RAS features, the tool is valuable, fast, and easy to c 

Redundancy Modeling 

The basic model assumption requires that, the system 
configuration be broken up into a series of device classes 
made up of devices or units or components, which may be 
field-replaceable units (FRUs) or more complex devices 
like peripherals. It is in these device classes or modules 
that the RAS tool allows any kind of k-out-of-n;G redun- 
dancy This means that at least k out of the a units present 
have to be working (to be G = good) for the specific device 
class to be considered up, or operational. For instance, a 
3-out-of-3:G printer constellation means all 3 printers are 
necessary, l-out-of-2:G processors is an assembly of two 
parallel (masking) processors, and a 92-oul-of-llS:G setup 
of workstations means that the device class of 1 IB worksta- 
tions is considered available and acceptable to the user 
even if as many as 26 of them have failed or are off-line, 

Another useful property that we believe is new in the 
RAS tool is the possibility of including device classes com- 
posed of several different but functionally equivalent units. 
The device class of 112 workstations may, for instance, be 



ACME Computer System 
(All numbers are fictitious) 

Contribution to System Failure Rate tor Each Device Class 



ACME Computer System 
(AH numbers are fictitious} 



AFR(%/year) 


15SB.67 


MTBF (hours) 


562,02 


MTTR (hours) 


4,74 


Cost(KS) 


92.70 


AFR'KS 


16.81 


Availability (%) 


99.16 


UPTIME (hours) 


6666.79 


DOWNTIME (hours) 


73.21 



Assumes 1 year = 8760 hours of operation 
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Printers 
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1(1) 


Tape Drives 
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Tape Drives II 
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" Total number of devices (Minimum that must be functioning) 
Fig. 2. Results of RAS tool calculations 
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The Language of Dependability 



The reliability of a system as a function of ftme R(t), is the 
probability thai the system will survive throughout the entire rnter- 
val [G\t], Reliability as a dependability measure is mosl useful 
for systems in which reparr cannot take pEace (space missions), 
or in which the whole system's operation is so critical that even 
short-lime repair or temporary degradations are impossible, or 
when repair is loo expensive. 

The availability of a system as a function of lime, A(t), \s ihe 
probability that the system is operational at time I If the steady- 
state availability exists, that is, if the system availability stabilizes 
over time at a constant value A, then A expresses the average 
time during which the system is capable of performing useful 
work Such useful work is done in a machine or system state [hat 
is acceptable to and defined by the user or customer Availability 
is useful as a figure of merit primarily in systems in which service 
can be delayed or denied for short periods while the system 
functions in a degraded mode 

Fault-tolerant computing means the continued correct execu- 
tion of a specific algorithm in the presence of operational faults, 
For this to work, the usually detrimental effects of defects must 
be overcome by employing hardware and/or software tech- 
niques Typical hardware methods include redundancy, but usu- 
ally software ana hardware have to work together for maximal 
effectiveness. As a consequence, fault-tolerant computing be- 
comes a topic of design and architecture rather than being an 
issue of hardware components. 

Trade-Offs 

White high reliability is a no-compromise business, high avail- 
ability and fault-tolerant computing require delicate trade-offs A 
trade-off needs to be made whenever a response variable V (for 
instance, customer satisfaction) depends on at least two conflict- 
ing variables, one of them pulling V up. the other pulling V down 
(for instance, the benefits versus ihe cosis of high availability) 
Trade-offs are unavoidable design decisions for complex sys- 
tems, especially at the high end, Ideally, trade-offs would op- 
timize every important technical and market-driven response vari- 
able under the constraints and boundary conditions of resources, 
expected return on investment, technology, and competitive 
pressure. 



composed of 80 HP 1 Pit) Computers and 32 HP Vectras, or 
any other combination. While the 15(1 and the Veclra have 
different failure rales, they are, as workstations, function- 
ally equivalent (even if, for example, one is faster than the 
other), 

Markov Chain Algorithm 

McGrady 1 describes an algorithm to calculate the avail- 
ability of a k-out-of-n:G subsystem with equal or different 
components. Instead of reproducing this procedure in its 
generality here, we give a simple example of a parallel 
subsystem, that is t a l-uut-of-2:G device class. Say thai this 
is the jth subsystem (of a total of m subsystems in the 
system) and has component availabilities A ]r i = 1 P 2 T 
which are computed from the respective failure rates a n 
and repair rates 4* } \' 



A„ = rV( ^ + A„J. 

The repair rate 4*\\ ' s related In ihe corresponding mean 
repair time MTTK,, by 

tfa - ^ftuVMTTRji. 

if a full duty cycle of ti7fi0 operating hours per year is 
assumed. The subsystem availability is 

A, = 1 - (1 - AjiJfl - %|: 

McGrady 1 S routine to compute such device class avail- 
abilities is based on the idea that A, can be written as 

Aj = A n A i2 [l f (1 - %1% + (1 - A#)/Afrl. (T) 

Generalizing to more complex k-out-of-n;G subsystems 
leads to similar expressions. For instance, the correspond- 
ing expression for a 2-out-of-3:G device class (sometimes 
encountered in power systems} is: 

Aj = A J: A i2 A J3 |l + (1 - A n ]/A n + (1 - A, 2 ]/A l2 + (1 - A| 3 )/A j3 |. 

(2] 

See reference 1 for more examples. 
Finally, the product of all the subsystem availabilities 
gives the total system availability A H ; 



II *j 

j-i 



The question is how to compute the system failure rate 
and system MTBF from the knowledge of the subsystem 
availabilities A, and the system availability A s . As it turns 
out T knowing the subsystem availabilities and their math- 
ematical forms [e.g., equations 1 and 2 above) is sufficient. 
The method comes from the theory of Markov chains, 

Markov chains are capable of modeling the different pos- 
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Markov Model 



Fi g . 3 . R AS tool block dta gram .The tool mo dels de gra dauons 
and redundancies at the CPU and peripheral levels, and 
computes reliability measures for general customer-defined 
configurations. 
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sible states (e.g.. of degradation) of a system. In particular. 
Markov chain tools allow the calculation of steady-state 
(matured) failure rates A 1 for all subsystems. The steady- 
state failure rate \, for the jth subsystem can be expressed 
as a weighted sum involving the d evict? availabilities A, 
and the repair and failure rates of the devices in the subsys- 
tem. The details are given in Ross. 2 For example, the failure 
rate of a parallel subsystem turns out to be: 

X, = ^jtAjiCl - A t2 ) + X i2 A i2 (l - An). 

where Ajj and A }2 can be expressed in terms of the 4> }] and 
Aji of the devices in the subsystem. 

The failure rate X ft of a system that is a series configuration 
of subsystems is the sum of the subsystem values A r 

Finally, the system MTBF T MTBF^ is given by the recip- 
rocal of K ■ 



MTBh\ 



i/jL, 



whlcfa is very familiar to engineers. 
The average repair time MTTR S for the system is equal to 

MTTR, = MTBF S {1 - AJ/A*. 

Fig. 3 is a block diagram showing the major building 
blocks of the RAS tool* 

Lotus Compatibility 

Currently the RAS tool is installed in an HP 3000 Series 
58 Computer. Future versions are being prepared to run 
on faster machines, A simplified algorithm for small to 
medium systems has been implemented in a Lotus 1 * 1-2-3" 
spreadsheet format to make its use mure interactive. 

As an example, Fig. 4 shows how the device class failure 
rates change when the configuration from the example in 
Fig. 2 is altered as follows: 



2 parallel SPUs instead of 1 SPU 

2 out of 3 4M-byte memory boards instead of 3 boards 

3 out of 4 tape drives I instead of 2 out of 3 tape drives I 

The Lotus version of the RAS tool can handle the relia- 
bility impact of such configuration changes as well as 
changes in component failure rates, repair times, or prices 
in a few minutes, giving the designer a quick evaluation 
of design trade-offs and providing the sales representative 
in the field with a hands-on tool to discuss RAS options 
with a customei 

A more ambitious enhancement to include stress factors 
(e.g, P temperature) is in its investigation phase. 
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Project Management Using Software 
Reliability Growth Models 

At HP's Lake Stevens instrument Division, the Goei- 
Okumoto software reliability growth model has provided a 
means for more accurately predicting the duration of system 
testing and determining when to release a new product. 

by Gregory A. Kruger 



SIGNIFICANT IMPROVEMENT in software reliabil- 
ity calls for innovative methods for developing soft- 
ware, determining its readiness for release, and 
predicting field performance. This paper focuses on three 
supporting strategies fur improving software quality. First, 
there is a need for a metric or a set of metrics to help make the 
decision of when to release the product for customer ship- 
ments. Second, accurately estimating the duration of sys- 
tem testing, while not directly contributing to reliability^ 
makes for a smoother introduction of the product to the 
marketplace. 

Finally, achieving significant improvement is easier 
given the ability to predict field failure rates, or perhaps 
more realistically, to compare successive software products 
upon release. Although this third strategy will be discussed 
in this paper, the emphasis will be on choosing the right 
software reliability metric and confidently managing the 
testing effort with the aid of software reliability growth 
models. For a more thorough discussion of estimating soft- 
ware field failure rates, see Drake and Wotting. 1 

Defects per Hour as the Quality Metric 
Two principal reliability measures are: 
Defect density expressed as defects per thousand lines 
of noncomment source statements [KNC3SJ 
Number of new defects found per hour of system testing. 
The first measure will show progress on successive soft- 
ware projects, while the second measure is more helpful 
for real-time project management. 

The user is unlikely to perceive a software product in 
terms of lines of code, since this is an invisible dimension 
of the product, The customer perceives quality in terms of 
how often a problem is encountered. This is analogous to 
a hardware product. We speak of the mean time between 
failures (MTBF) for repairable instruments, and in fact, 
failure rate and its reciprocal, MTBF, are appropriate mea- 
sures of software quality. It is likely that, however different 
the sizes in lines of code, if two software products have 
the same failure rate during operation t they will be per- 
ceived as having the same level of reliability by the customer. 
Although ascertaining the customer-observed failure rate 
is difficult at best, measuring the failure rate during system 
testing is straightforward. For each successive time period, 
a week for example, record both the engineering hours 
invested in testing and the number of new defects found. 



Dividing the number of defects by the hours spent testing 
gives an estimate of the "instantaneous" new-defect-f hid- 
ing rate. As defects are found and removed from the soft- 
ware, this rate declines. The decreasing new-defect-finding 
rate and the corresponding increasing mean time between 
finding new defects provide the criteria by which progress 
and the ultimate conclusion of system testing can be 
judged. 

Reliability Growth Modeling 

There is considerable statistical literature on modeling 
the reliability growth process of finding and fixing defects 
in a software product. Available models can be classified 
into four broad categories: 2 

Input-domain-based models 

Fault-seeding models 

Time-between-failures models 

Fault-count models. 

Inpul-dnmain-based models can be characterized as a 
basic sampling theory approach. One defines the input 
domain of a software system to be the set of all possible 
input data the software will encounter during operational 
use. The reliability can then be estimated by taking a rep- 
resentative sample from the input domain and looking at 
the resultant failure rate when the sample data is input to 
the system for execution. Although the input domain sam- 
pling procedure will not be strictly followed here, the con- 
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cepts embodied in this technique turn out to be helpful in 
designing the testing process and understanding the limi- 
tations of making estimates of field performance based on 
defect-finding rates observed during system testing. 

The approach in fauft-seedmg models is to seed a known 
number of defects into a software product. The software is 
then tested using a process that presumably has equal prob- 
ability of finding a seeded or an indigenous defect, The 
numbers of indigenous and seeded faults found are used 
to estimate the system reliability. This concept has several 
implementation problems and draws considerable skepti- 
cism from project managers and designers alike. 

The preferred data in time-between-fai lures models is 
the observed running times between successive failures. 
Typically it is assumed that the time between two consecu- 
tive failures follows a distribution whose parameters are 
dependent upon the number of defects remaining in the 
software during that time interval. These distribution pa- 
rameters can be estimated from the observed time-between- 
fai lures data. However, the actual running time between 
successive failures may be a difficult measure to extra* I 
from the software testing process, 

In fault-count models, the variable of interest is the 
number of defects observed per specified time interval. 
The basic idea is that the number of failures observed per 
time interval can be modeled according to a Poisson pro- 
cess. The Poisson distribution is widely used to model the 
number of occurrences of some event in a time or space 
interval, As previously noted, the defect-finding rate statis- 
tic is relatively easy to capture. For this reason, the focus 
of this paper will be on the category of fault-count models. 

John Musa's execution-time model, discussed by Drake 
and Wolting, 1 is an example of a fault-count model. Al- 
though their derivation differs from Musa's. Goel and 
Okumoto 3 propose essentially the same model. Because of 
its simplicity urn I inluitiveness and thisauthor's preference 
i'i »r ihe derivation of the model from the underlying assump- 
tion of a Poisson process, the Goel-Okumoto model will 
be used here, 

The Goel-Okomoto Model 

The model assumes that between software changes, the 
number of defects observed during each hour of test will 
follow the Poisson distribution with a constant average, A, 
Note that this defect rate includes the occurrence of both 
new and repeat defects, If the observed defects are fixed 
immediately, the Poisson parameter A will he decreasing 
every hour of testing (or at least decreasing with every code 
correction). In this case, we have a nonhomogeneous Pois- 
son process and can expect the defect-finding rate to de- 
crease exponentially until software changes cease. Con- 
versely, as the defect rate declines, the cumulative number 
of defects found will increase and asymptotically approach 
a constant. Since there is a lag in the correction of defects, 
the ideal statt! of instantaneously fixing defects when found 
can be approximated by counting only the new defects 
found each hour 

The functional form of the model is as follows. The 
cumulative number of defects found by time i is giver by 



and the instantaneous new-defect-finding rate at time i is 
given by 

l(t) = ra'(t] = ape"* 1 . 

In this model, a. is the expected total number of de' 
in die software. f$ is the initial defect -finding rate divided 
by the expected total number of defects in the software, t 
is the cumulative hours of system testing, and ni'[t) denotes 
the derivative of m(t). 

Figure 1 shows the typical shape of l(t) and m(t). 
It should be noted that the parameter a is actually an 
estimate of the total number of defects the system testing 
process is capable of detecting. The input domain concept 
is useful for understanding this point. If there is some 
stratum of the input domain that the system testing process 
is not covering* the parameter a is actually an estimate of 
the total number of defects minus those generated by data 
from that specific stratum. 

There are three key assumptions underlying this model: 
* All new coding is completed before the start of system 
testing. 

Defects are removed with certainty and without intro- 
ducing new defects. 

Testing methods and effort are homogeneous, that is. the 
last hour of testing is as intense as the first, 
The first assumption ensures that the defect-finding rate 
follows a monotonia- pattern. While the second assumption 
is difficult to guarantee, instituting strict security measures 
will minimize the rate of defect introduction. The phasing 
of distinct test activities, the introduction of new testing 
tools, the addition of new team members, and software 
testing burnout make meeting the third assumption dif- 
ficult, as well. However, good planning of the testing pro- 
Cess and averaging the defect-finding rale over a calendar 
week will tend to smooth out differences in testing effort 



Fitting the Model on a Project 

Consider the application of the Goel-OkumotO model to 
a 90-KNCSS firmware product This [imiluct had com- 
pleted system testing eighteen months earlier, but records 
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were available from which a data set could be constructed 
to test the feasibility of the model. For each of 43 weeks 
of system testing, the cumulative test hours, the instantane- 
ous new- defect -fin ding rate, and the cumulative number 
of defects found were calculated- However, because the 
number of hours invested in system test during the first 
several weeks of testing was highly variable, the first 19 
weeks of data was condensed into five data points, each 
representing approximately 200 hours of testing. This left 
29 sets of data values, each representing at least 200 con- 
secutive hours of testing. The nonlinear least squares fit of 
the the l(t) model to this data appears inFIg< 2. The obvious 
departure of the actual defect- finding rate data from the 
model at about 9000 cumulative hours appears to be 
explained by the addition of several experienced designers 
to the testing team. This is a clear violation of the assump- 
tion that the testing effort remains homogeneous. Using 
the parameter estimates obtained from fitting l[t), the m|t) 
model and cumulative data are plotted in Fig. 3* Parameter 
estimates and standard errors from fitting l(t) are; 



Estimates: $ 



2307 
0*000384 



Standard Errors: ov, = 362 



0.000102 



There are some statistical advantages to fitting l[t) to the 
defect -finding rate data and using the resulting parameter 
estimates to plot m[t), By the nature of cumulative data, 
the residuals from fitting m(t) are serially correlated. How- 
ever, the residuals from fitting l[t) to the instantaneous 
defect-finding rates are well-behaved, more closely con- 
forming to the standard assumption of normally distributed 
errors with mean zero, On the other hand, l[t) can only be 
fit if there are enough hours of testing and defects found 
to generate reasonable estimates of the instantaneous de- 
fect-finding rate. If a project will be undergoing a small 
number of test hours so that blocks of time (such as a week) 
cannot be aggregated to give an estimate of the defect-find- 
ing rate* then it is difficult to obtain parameter estimates 
by the method of least squares. In this case, maximum 
likelihood estimation should be used to fit the m(t) fund ion 
to the cumulative data, This should not be interpreted to 
mean that maximum likelihood estimation is an inferior 
technique to be employed when the least squares method 
is not feasible. On the contrary, it can be demonstrated that 
maximum likelihood estimates are among the best, The 
intent here is to show what is possible with simple least 
squares procedures. 

Using the Model in Real Time 

The project discussed above, project A P demonstrates 
that real software testing data does follow the non homoge- 
neous Foisson process model. Another project, project 0, 
an application software product consisting of 100 KNCSS, 
demonstrates the potential for using the defect-finding rate 
metric: and the Goel-Okumoto model in managing the sys- 
tem testing process. 

Setting the Release Criteria. First, a defect -finding rate ob- 
jective must be incorporated into the existing set of software 
release requirements. The total release checklist would 
then include such things as all functionality complete, all 
performance specifications met, no known serious defects, 



and the defect- finding rate criterion. The question remain- 
ing is, "What would be an appropriate value for this defect 
rate?" Project A concluded system testing wilh a rate of 
0.03 defects per hour. As it turned out T the firmware proved 
to have extremely high reliability in the field. However, it 
is not clear that a 0.03 defect-finding rate is appropriate 
for all projects — particularly on software rather than 
firmware projects. In addition, arguments can be made that 
the relative severity of each defect found should be used 
to express the release criteria in terms of a weighted defect- 
finding rate. Ultimately, based upon work at other HP di- 
visions and the belief that the defect weighting would be 
quite severe, a release requirement of 0.04 weighted defects 
per test hour was established. A team of four project en- 
gineers reviewed all defects found and scored each on a 
1-9 scale with 1 being a cosmetic Haw and 9 being most 
serious. These scores were then normalized by dividing 
by nine. The team recognized that the software reliability 
modeling process would be dependent upon their ability 
to weight the defects consistently throughout the duration 
of the project. 

Applying the Model. Once the release requirement is set, 
the Goel-Okumoto model becomes useful as a real-time 
management toot. Having fit the model to the most recent 
data* an estimate for the total test hours required to release 
is obtained by calculating where the curve crosses the re- 
lease requirement. The management team can then review 
weekly plots of the defect-finding rate and estimates of the 
total test hours to release. These pints enable the team to 
predict the project's conclusion more accurately and to 
judge progress along the way. in addition, the engineering 
staff conducting system testing may find the data to be 
morale building* Without some measure of progress, system 
testing is nflen perceived to be a never-ending process. 

Project B ultimately reached the release goal after 5760 
hours of testing over 22 calendar weeks. Figs. 4 and 5 show 
the final plots of the I(t) and m{t) functions fit to Project B 
system testing data using nonlinear least squares. The pa- 
rameter estimates and standard errors from fitting l(t) are 
as follows: 
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Fig. 3. Cumulative defects each week for project A firmware 

and fitted curve. 
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One can simplify the analysis by first taking the natural 

log of the defect-finding rate data. A quick check of the I[t) 
function % T erifies that taking the log transformation makes 
the model linear. Fig. 6 shows the regression line fit to the 
logged data. The benefit of operating with a linear model 
is that standard techniques can be used to set confidence 
bounds on the estimates of total system test hours needed 
to reach the release requirement. 4 As can be seen from Fig. 
7, the estimates of total test hours proved to be both consis- 
tent and accurate. The project concluded system testing 
within 300 hours of the estimate made eleven weeks previ- 
ously* The accuracy of these estimates far exceeded our 
expectations. 

Estimating Field Failures 

Two different but related estimates can be made concern- 
ing the number of defects the customer will observe in the 
first year after release. The first is the expected number of 
times the user wi 1 1 experience a defect in the system inclu d - 
ing the recurrence of a previously observed defect. The sec- 
ond is the expected number of unique defects encountered . 
New And Repeat Defect Rate. Including repeat defects in 
estimates of a software system's failure rate can be consid- 
ered as taking the perspective of the customer. The cus- 
tomer perceives quality in terms of how often a defect is 
encountered. The user may not be able to tell if two or 
more defects are actually caused by the same fault in the 
code. Estimating the combined new and repeat defect rate 
may be possible using the l(t] model fit during system test- 
ing, Recall that the defect -fin ding rate declines because 
defects are removed from the software as they are found. 
Once QA ends* defects are no longer being fixed. Therefore, 
defects should now be observed according to a Poisson 
process with a fixed average number of defects per hour 
of use. The best estimate available for this is the defect-find- 
ing rate observed at the conclusion of system testing — U.Q4 
weighted defects per hour for Project B. The question is. 
how good an estimate is this for what the customer will 
experience; f 

The answer is dependent upon the difference between 
a customer's and an HP engineer's defect-finding effi- 
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Fig. 4. Weighted weekly defect -finding rates for project B 
software. 
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Fig > 5. Cumuiatt ve defects each week for project B software 

ciency, It is likely that the test engineer, rather than the 
customer is the more efficient at finding defects. If an 
acceleration factor, K, could be identified expressing the 
ratio of customer hours to one hour of system testing, the 
customer's defect-finding rate could be estimated as IfK 
times the defect rate at the conclusion of system testing.* 
The simplifying assumption here is that the customer's 
defect-finding efficiency is constant. Evidence to date 
suggests that a value of K = 50 might be a reasonable esti- 
mate of the acceleration factor. The estimated combined 
new and repeat defect -find ing rate would then be 0,0008 
weighted defects per hour. 

New-Defect-Finding Rate. Considering the number of 
unique defects in the software is taking the perspective of 
I he company; we want to know T if code revisions are likely 
to be necessary. When considering repeat failures, we 
stopped on Ihe exponential Failure rate curve and extended 
a fixed horizontal line at the defect-finding rate observed 
at the end of system testing. To estimate unique failures, 
we continue down the l(t) curve for the first year after 
release to shipment, Conversely, the cumulative number 
of defects found should continue up the m[t) curve. 

Ii is a trivial matter to estimate the number of unique 
defects the system testing team would find in an additional 

*5ee referernju I lur another perspective? cm accploralic-n factors. 



&£ 



90 D p Confidence Bounds 
on the Regression Line 




1000 



2000 3000 4000 

Cumulative Test Hours 



5000 



6000 



Fig, 6. Linear model tor project B after log transformation 



JUNE 19aa HEWLETT-PACKARD JOURNAL 33 



)Copr. 1949-1998 Hewlett-Packard Co. 



t hours of testing. This is simply: 

m(T + t) - (total found so far) 

where T is the total number of test hours accumulated IP 
date. In fact, the total number of defects remaining in the 
software can be estimated by 

a — (total found so far). 

For the Project B software, 1239 - 1072 = 167 weighted 
defects. 

An estimate of how many of these remaining defects will 
be found in the field can be accomplished by defining a 
n^w m(t) function for the customer. The a parameter can 
be set equal to the estimated number of defects remaining, 
The ft parameter can be obtained if we are once again will- 
ing to use an acceleration factor. For the Project B software, 
the cumulative number of unique defects found after t 
hours of customer use is given by: 



m(t] = 167(1 - e 



■<i imiMH-ji-iiJi 



For example, in 2000 hours of customer use we could 
expect 3.2 unique weighted defects to be found. 
Field Failure Rates in Perspective. Estimates of customer 
defect rates and defect totals should be taken with a healthy 
degree of skepticism. It is a risky proposition to propose 
acceleration rates and attempt to estimate customer ob- 
served defect rates. Experience La date clearly shows that 
our customers are not finding as many defects as estimated 
by the model. However, these two projects have not pro- 
vided sufficient insight for us to estimate confidently the 
relationship between model projections and customer ex- 
perience. After more products are released according to 
this process, comparisons will be drawn between the defect 
rate observed in-house and the subsequent rate in the field. 
An empirical estimate of the acceleration factor may be 
obtained in this way. To date, the results from these two 
projects are as follows. 

Recall that on Project A the analysis was conducted on 
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unweighted data — all defects were counted equally regard- 
less of severity. From the model, the number of unweighted 
defects estimated to be remaining in Project A is 2307 - 
2167 = 140, After 18 months in the field, there have been 
three customer reported defects. All of ihese have been 
resolved to the satisfaction of our customers. Three defects 
found out of an estimated 140 remaining suggests a ratio 
of estimated to actual defects of 47 to 1. 

As has been stated, the estimated number of weighted 
defects remaining on Project h 1 is 167. After 12 months in 
the field, there have heen four defects identified, all of 
which have been resolved in a subsequent release. After 
applying the same weighting criteria that was used during 
system testing, these four defects represent 2,4 weighted 
defects. Since this product has not been in the field as long 
as Project A. we could multiply by 1.5 to give 3.6 weighted 
defects after 18 months. The ratio of estimated to actual 
defects would then be 46 to 1 — amazingly close to the 
results from Project A but still conjecture at this point. 

Conclusions 

The defect-finding rate metric and the Goel-Okumoto 
model proved, on one project, to be a real contribution to 
managing the system testing phase of the software develop- 
ment process. Judging from the fit to historical data, the 
model would have been helpful on an earlier project as 
well. Through Weekly plots of the data and the model, the 
management team received feedback on testing progress 
and information to aid resource planning. Ultimately, there 
was a clearer decision as to whether to release the product 
for customer shipments. 

In cases with limited test hours, the straightforward ap- 
plication of the Goel-Okumoto model would be to fit the 
m(t) model to the unweighted cumulative defect data using 
the method of maximum likelihood. If there are enough 
test hours to allow r reasonable estimates of the instantane- 
ous defect-finding rates, fitting Ift] usinf" least squares tech- 
niques is possible, although not superior to the maximum 
likelihood procedure. Based upon the results from these 
two projects, it appears that the modeling process works 
for both unweighted and weighted defect data, Extrapola- 
tions from the model to field conditions have proved to 
overestimate the number of defects customers will en- 
counter. Further study is required to estimate confidently 
the relationship between model estimates and customer 
observed failure rates. However, we are continuing to 
tighten the failure rate release criterion as new products 
enter final system test. 

An alternative approach to making the release-to-produc- 
tion decision would be to plot the mean time between 
finding new defects, l/l(t). This metric could then be used 
to make an economic decision on when continued testing 
w*ould be too costly. 
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A Reliable, Autoloading, Streaming 
Half-Inch Tape Drive 

Designed for rack mounting, this compact tape drive 
cleverly channels air flow to load different-sized half -inch 
tape reels automatically, It also features higher performance 
and reliability than its predecessor. 

by John W. Dong, Kraig A. Proehl, Ronald L. Abramson, Leslie G. Christie, Jr., and Douglas R. Domel 



HP's NEWEST HALF-INCH TAPE DRIVE, ihe HP 
7980A (Pig. 1) t is an autoloading, reel-to-reel, hori- 
zontally mounted, streaming drive that reads and 
writes in two standard nine-track formals: B250 GCR and 
1600 PE. ia This maintains compatibility with previous 
drives and tapes, an important feature because half-inch 
tape is still a significant standard in the computer industry 
for backing up, archiving, and interchanging computer 
data.^ 

The HP 7980A provides computer system users with a 
reliable, low-cost, backup device for large amounts of on- 
line disc memory. It does this with higher performance 
and reliability and lower cost than its predecessor, the HP 
7978B, which it replaces. It reads and writes at 125 inches 
per second, roughly 60 percent faster than the HP 7978B, 
The HP 7980A can rewind a 2400-foot tape in less than 90 
seconds, reducing overall data transfer times significantly. 



The HP 7980A is 40 percent more reliable, a result of the 
increased use of VLSI components to reduce parts counts 
even further than was achieved on the HP 7978B, The 
monthly maintenance cost of the HP 7980A is half that of 
the HP 7978B, Horizontally mounted in a standard-width 
rack cabinet, the HP 7980A is 8,75 inches high, a third the 
size of the HP 7978KJ. Phis saves valuable floor space and 
allows better use of rack cabinets. 

The horizontal mounting means that the user normally 
cannot access the tape path to load the tape manually. 
Hence, an autoloading feature was designed into the HP 
7980A, The operator simply places the tape reel in the door 
opening and closes the door. The autoloading sequence 
starts automatically once the door is closed, leaving the 
operator free to do other tasks while waiting for the tape 
to load in a nominal time of half a minute. There is no 
need for an EZ-LOAD cartridge around the tape reel, which 




Fig. 1. The HP 7Q80A Tape Drive 
is a half-inch, reel-to-reel, 125-tps 
streaming tape drive designed for 
systems with disc backup require- 
ments of greater than 400 mega- 
bytes. Uke its predecessor, the 
HP 7Q78B, this high-performance 
drive operates with both 6250 
GCR and WOO PE standard for- 
mats . The HP 7B&0A can automati- 
cally load and thread any stze reel 
ranging tn diameter from six to 
10.5 inches. This drive can be 
horizontal fy mounted m a 19-mch- 
wide rack enclosure for better 
floor-space utilization, and it sup- 
ports IBM/ ANSI-compatible for- 
mats for software distribution and 
data interchange between HP and 
non-HP systems. 
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is required for autoloading on the HP 7976A Tape Drive. 
The HP 7976A autoloads only 10.5-inch tape reels: the HP 
7980 A autoloads all standard half-inch tape reels from six 
to 10.5 inches in diameter. The earlier HP 7978 A 'B Tape 
Drive must be loaded manually and cannot autoload tapes. 

In addition, the HP 7980A supports the use of 3600-foot 
half-inch tape, under certain guidelines, which the HP 
B doesn't. This increases data capacity by 50 percent 
over standard-length tape reels. 

The HP 7980A was developed and released in 40 percent 
less time than the previous tape drive. This was mainly a 
result of two factors. The first was keeping together an 
experienced core group of engineers from the HP 7978A/B 
development team to architect and design the HP 7980A 
drive. The second was concentrating on a core system de- 
velopment of the HP 7980A, that is, having a well-defined 
product and not adding additional features and configura- 
tions along the way. Such added features tend to prolong 
development cycles. This does not mean that these addi- 
tional features are not eventually added, but they are 
worked on according to need after the core or base system 
is released. 

Tape Path 

The HP 7980A is a totally integrated tape drive, simul- 
taneously incorporating a small form factor, an autoload 
feature, design-for-assembly concepts, low cost, and high 
reliability. The major design objective established to ac- 
complish these goals was design simplification. 

The HP 7980A has a very simple tape path as shown in 
Fig. 2. There are only two rolling elements: the speed sensor 
and the buffer arm roller. There is only one additional fixed 
tape guide. The oxide side of the tape contacts only the 
tape cleaner, the magnetic tape head, and the tape displace- 
ment unit. The tape displacement unit, located between 
the tape head and the tape cleaner, contacts the oxide side 
of the tape only during repositioning and while the tape 
is stopped. The tape displacement unit pushes the tape off 
the very smooth surfaces of the head and the tape cleaner 
to prevent the tape's sticking to these smooth surfaces dur- 
ing high temperature and humidity conditions. 

The buffer arm assembly [buffer arm, spring, and roller) 
helps take up slack in the tape during servo starts and 
stops. It also establishes the tension on the tape, The buffer 
arm roller and fixed guide, along with the speed sensor, 
guide the tape in a precise manner over the head. The 
speed sensor measures the velocity of the tape and feeds 
it back to the servo system. 

The half-inch tape reel is centered, seated, and locked 
by the supply hub. The autoload blower forces air through 
the door. The louvers in the door direct the air onto the 
tape reel to lift the end of the tape. The tape end is then 
carried by the air flow around the buffer arm roller and 
the fixed guide, It then goes over the tape cleaner and the 
magnetic tape head. The tape is finally sucked onto the 
take-up reel after passing around the speed sensor. 

The drive motors are located directly heneath and are 
attached to the supply and take-up hubs. The blower is 
located between the two drive motors. 



Integrated Autoload and Tape Path Design 

The HP 7980 A tape path is very simple. This greatly 
reduces costs and assists the design of the autoload mech- 
anism. It Simplifies the (ask of blowing the tape end off 
the supply reel and threading it through the tape path and 
onto the take-up reel. 

The casting is an integral part of the tape path and au- 
toload mechanism. The walls of the casting help determine 
the way the tape is blown around by the air flow created 
by the blower. The surface roughness of the casting is very 
important, The surface of the casting around the supply 
reel area must provide enough friction to permit the supply 
hub to autocenter the smaller tape reels. Conversely, the 
casting surface near the buffer arm assembly must be 
smooth enough so that the tape doesn't catch on any surface 
features while autoloading. The buffer arm shape is critical 
to autoload and servo success. The design of the air deflec- 
tor on the speed sensor is crucial to how well the tape end 
moves around the speed sensor roller and attaches to the 
take-up reek 

The autoload success rate is greatly affected by the design 
of the door ramp and the door louvers. These are designed 
so the HP 7980A can autoload all standard half-inch reels 
with varying amounts of tape on each reel. 

Holes in the top cover eliminate the need for gaskets or 
tight tolerances on the door-to-bezel fit and the top cover-to- 
casting fit, They do this by eliminating air flow reversal, 
which can occur when autoload air flows out the door-to- 
be^el and top-cover-to-casting cracks, rather than down the 
tape path. The ribs bb the lop cover fit inside the casting 
to reduce air leakage and allow loosening of the top-cover- 
to-casting fit. 
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Fig. 2. HP 7980 A tape path components 
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The speed sensor is placed in the tape path for maximum 
tape wrap on the roller. This prevents the tape from slipping 
on the roller during tape acceleration and deceleration. 

The buffer arm and its spring are designed and placed 
SO I hat tension variations caused by the tape coming off 
the supply reel al different angles relative to the aim are 
minimized- An additional roller could have been used to 
feed tape at a constant angle to the buffer arm, but was not 
put in for several reasons: cost, simplicity, and ease of 
autoloading the tape, 

Air System Design, The autoload air system was designed 
into the deck casting and the tape pain from the beginning. 
Reduction of air ductwork was a major cons i Herat ion. The 
less the air has to pass through winding passages, the lower 
the air pressure and volume requirements. For this reason, 
the tape path needs to be simple, Reducing the air pressure 
and volume reduces the autoload air blower's size, noise, 
and cost. 

The blower is placed between the two drive motors since 
that provides the shortest air passage. The air is sucked 
into the take-up reel, into the blower, and then blown out 
of the blower into the door area to repeat the path [see Fig, 

The door is used as an air passage for several reasons. 
No space is taken up by air ducting to blow air into the 
tape path, thus allowing the maximum amount of width 
for the card cage. In addition, using the door as an air duct 
allrjws the air to blow across the entire front of the tape 
reel for maximum design flexibility. The door louvers were 
then designed to blow air to autoload all standard reel si^es, 

The air velocity required to move the tape down the tape 
path was calculated from momentum and Bernouli equa- 
tions. The pressure drop was calculated in a similar man- 
nen The blower was then selected based on these calcula- 
tions and the tape path, with the blower, was measured 
for pressure drops and velocities. The measurements were 
fairly close to the theoretically calculated pressure drops 
and air volumes. 

Autoload Algorithm. The HP 7980 A autoload algorithms 
are designed to load all sizes of half-inch tape reels in a 
minimum amount of time, while making sure that the user's 
tape is treated with as much care as possible, The autoload 
process is a combined effort of firmware, mechanical de- 
sign, and electronic sensors. Four sensors are used by the 
autoload process to monitor its progress and to indicate 
when an error condition is present. The first of these sen- 
sors is the door sensor, which consists of two micro- 
switches connected in series that detect whether the front 
door or top cover is open. When this sensor detects a door 
closure, the autoload process begins. Opening the top cover 
at any point halts the autoload. Under the supply hub is 
the reel encoder sensor. This optical sensor detects the 
three reel encoder flags attached to the supply hub. These 
flags interrupt the sensor beam when a reel is properly 
seated on the supply hub, The lape-in-path sensor is an 
optica] sensor which is positioned across the tape path 
preceding the head. The tape leader will interrupt this 
sensor beam when it enters the tape path. The final sensor 
is the optically encoded speed sensor. This sensor provides 
two quadrature pulse trains which are decoded to give 
position and velocity information. 



Once a door closure is sensed, the autoloading operation 
begins. The first task is to detect whether a reel is already 
threaded through the tape path, a condition that is most 
likely to be present after a power failure. This condition 
is detected by the tape-in-path sensor and by turning both 
the supply and take-up motors in opposing directions at a 
low T rate. If a tape is already threaded, activity will be seen 
on the speed sensor. In this case, the tape does not need 
to be automatically threaded and the servo loops can be 
closed. 

If a tape is not already threaded, the load fan (blower) 
is turned on and the supply hub is slowly rotated in a 
I nunterclockwise direction to center and seat the reek At 
this {joint, the reel encoder sensor is checked continuously 
for evidence that the reel is seating itself on the hub prop- 
erly. When the reel is seated correctly, all three reel encoder 
flags will interrupt the sensor beam once per revolution, 
Time interval measurements are made ttj ensure I hat all 
three flags [not just one or two) are present, and that they 
have the correct relationship. If all three flags are not de- 
tected, the supply hub is shaken back and forth quickly in 
an attempt to get the reel to locate itself properly tm the 
hub. A reel that will not properly seat after shaking a second 
time will be rejected and a MISLOAD will be reported. 

The reel encoder sensor and flags are also used to regulate 
the speed of the supply hub during these open-loop oper- 
ations. The dc motors must turn at relatively low speeds 
when centering a reel and feeding the tape. These speeds 
require relatively low voltages at the motors. Because of 
voltage offsets, motor constants, and temperature changes, 
the rotational speed could not be adequately controlled by 
a single voltage command. Hence, the voltage command is 
constantly adjusted based on the time measured between 
each reel encoder flag pulse. This provides rather gross, 
but sufficient control over the rotational speed of the motor 




Fig. 3. HP 7980 A tape path showing air flow for automatically 
loading and threading a tape. 
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in its open-loop mode. 

After the reel is found to be seated properly on the hub, 
it is ready to be locked into place. First, however, the tape- 
in~path sensor is polled to determine if the tape is 1 
seen continuously in the tape path. If t: case, the 

reel has been loaded upside down. When the reel is spin- 
ning counterclockwise, the tape should only be momentar- 
ily flopping into the tape path. If the reel is found to be 
upside-down, the supply hub is turned in the clockwise 
direction until the tape clears the path, The door is then 
opened and the INVERT message is reported. 

If the reef is not inverted, the hub locking routine can 
proceed. The supply hub motion is stopped and the hub 
lock solenoid is engaged. The hub is then rotated clockwise 
by applying a voltage ramp. As the hub rotates, three lock- 
ing feet come up out of the hub and grab the tape reel, 
holding it securely to the hub. During this locking opera* 
tion. the reel encoder flags and sensor serve as a check on 
whether the lock is a success. If the reel is properly locked 
on the hub, one of the encoder flags will interrupt fife 
sensor beam as the voltage ramp is applied. After a success- 
ful lock, the hub lock solenoid is disengaged and the supply 
hub is free to spin again. This locking cycle is described 
in more detail later, 

The supply hub is now rotated again is the coun- 
terclockwise direction in preparation for the inertia check. 
The inertia check is used to determine the size of the reel 
being loaded. This information is used to set up some au- 
toload and servo parameters that can be optimized accord- 
ing to the reel size, The inertia check is performed by ap- 
plying a step voltage to the supply motor at the instant a 
reel encoder flag breaks the sensor beam, The time it takes 
for the next reel encoder flag to come arunod and break 
the sensor beam is proportional to the inertia of the reek 

As the hob continues spinning counterclockwise, some 
additional error conditions are checked. The check for an 
upside-down reel is repeated. A check for the tape leader 
being stuck to the reel is done. In the same manner thai 
tape continuously in the tape path signals an inverted reel, 
a tape that never breaks the tape-irvpalh seosor beam indi- 
cates that the leader is stuck to the reel This can often 
occur because of static electricity or because the leader is 
jammed under the reel flange. If the tape leader is indeed 



stuck, the supply reel is spun counterclockwise at high 

speed in an attempt to free the tape end. Failure to free the 
tape leader at this point will abort the load process. 

Next, the tape is ready to be threaded. The supply hub 
continues to spin in the counterclockwise direction, and 
the tape-Ln-path sensor is monitored for the tape leader. 
Once the tape leader is sensed in the tape path, the supply 
reel continues to spin for another half second (to pull the 
end of the tape back to the beginning of the tape path]!, at 
which point it reverses direction and starts feeding the tape 
down the tape path assisted by the air flow. The t ape- in- 
pat h sensor is checked continuously to make sure the tape 
stays in the path during the feeding operation. If the sensor 
fails to detect the tape in the path, an error condili' 
flagged and the load is aborted. 

As the threading proceeds, the speed sensor is monitored 
for activity, if the tape has been correctly fed down the 
tape path and has caught onto the take-up reel, the speed 
sensor begins to spin. The speed sensor is then used to 
calculate the amount of tape that is wrapping around the 
take-up reel. If activity is not seen at the speed sensor 
within a certain period, or the required number of wraps 
are ncil completed within another certain period, the tape 
will be pulled back out of the tape path by the supply reel 
and the autoload will be retried, The HP 7980 A will attempt 
to load the tape in this manner five times before reporting 
failure to do so. 

After it has been determined that the tape is properly 
threaded, the load fan can be turned off and the servo loops 
closed, The take-up and supply motors are driven slowly 
m oppositioo to tension the tape. The tension arm position 
is monitored by an analog-to-digital converter (ADC) and, 
once it has reached the center (UV] position, the tension 
loop is closed. The tension integrator is turned on and the 
loop is given about a half sruiml to stabilize. Next, the 
microprocessor closes the velocity loop and digitally con- 
trols its operation. The tension shutdown circuitry is en- 
,ihh i] In prevent the drive Irom damaging itself or the tape 
during some sort of catastrophic failure, Any failure to 
establish tension during this process will cause the load 
operation to be aborted and a MISLOAD reported* 

The autoload process takes *AU seconds if all goes welt. 
Then the drive can be put on-line and all normal reading 
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Fig. 4* Cross section of supply 
hub. The foot is rotated upward 
and outward from the supply hub 
by holding the hub fixed and rotat- 
ing the sleeve underneath. This 
causes the dumbbell -shaped tie 
rod to move from a sf anted posi- 
tion io a position more upright, 
thus pushing the foot up and out, 
Reel encoder flags pushed down 
by the presence of a reel are 
sensed by the optical sensor to 
verify proper reel loading. 
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rind writing operations can be performed. 

Hub Lock Mechanism 

The first major objective in the autoload process is to 
determine the presence of. seat, and secure reels of tape. 
This process is complicated by the availability of four dif- 
ferent si sees of half-inch-lape reels ranging in diameter from 
six inches ttj 1 Qt$ inches, the latter being the most common 
size. The tape reel is secured onto a supply hub which 
suspends the reel at the proper height, holds it secure, and 
rotates the reeh 

When the door of the tape drive is closed, the autoload 
algorithm starts immediately* The supply hub rotates in a 
clockwise direction and an off-center protrusion atop the 
hub contacts the inner race of the reel of tape. The rotational 
motion drags the reel in an inward spiral pattern until the 
reel eventually drops onto the hub's mounting plane. Three 
small plastic springs (the reel encoder flags mentioned ear- 
lier) mounted around the mounting plane are pushed down 
by the weight of the reel and trip the reel encoder optical 
sensor. Thi ush fit three springe assures planar contact qJ 
the reel on the hub, which helps prevent tape edge damage 
by keeping the reel as parallel as possible to the tape path- 
A fourth spring may be deflected by the write enable ring 
on the reel of tape to relay that status as welh 

At this point, the reel is secured to the hub by means of 
three feet that rotate up and out of the hub to clamp the 
inner race of the reel, The feet are driven by a four- bar 
linkage that translates the rotational motion of the hub into 
a vertical rotation of the feet. At the end of their travel, the 
feet are locked into position in much the same manner as 
a toggle switch with an overcenter design. 

The mechanism that drives the feet up and into place 
has foor basic parts [Fig, 4j. The first is the foot itself, 
which pivots about an axle swaged into the hub base. The 
foot is driven by a barbell-shaped tie rod which translates 
the horizontal rotation of the hub into the vertical rotation 
of the foot. The barbell, in turn, is driven by a sleeve-shaped 
part that fits around the base of the hub and is held in 
place by a spring and retaining ring, in normal operation, 
the sleeve rotates with the hub. During the reel locking 
stage, a solenoid engages an arm that stops the motion of 
the sleeve, and while the hub is rotated, causes the sleeve 
to rotate relative to the hub. The foot is pushed up and 
into place, depressing the wave spring when the foot con- 
tacts the tape reel. As the hub rotation continues, the barbell 
pops over center and starts to retreat. The relative motion 
at this point is ended by stops in the hub that limit the 
motion of the sleeve. The lock cycle is double-checked by 
the reel encoder optical sensors. Sensor placement is such 
that one of the plastic springs always interrupts the sensor 
at the moment the sleeve hits its stops, A neoprene pad 
embedded in the surface of the foot ensures that the reel 
does not slip relative to the hub. The hub is unlocked by 
simply reversing the process. In the case of a power failure, 
the entire lock or unlock cycle can be performed by man- 
ually turning the reel and engaging the arm, 

Many obstacles had to be overcome to ensure a successful 
design. The first of these was materials selection. The ma- 
terial for the plastic springs must have a high endurance 
limit, low creep at moderately high temperatures, and a 



consistent and high spring constant. The material chosen 
was ULTEM 1000 by General Electric. Material for the slid- 
ing parts on the hub was tailored to minimize wear and 
friction, and to resist deformation under load. The second 
obstacle w T as cost, We made extensive use of design for 
manufacturability to decrease part count and minimize as- 
sembly time. The result is a hub with no screws that is 
completely assembled from one side. 

Integrated Tape Path 

All the tape path components, except for the reel motors, 
are mounted onto a precision head plate which is mounted 
on the deck casting. The magnetic tape head and the tape 
cleaner are permanently mounted onto the head plate in a 
manner similar to that used for the earlier HP 7978A Tape 
Drive. 4 This greatly reduces the number of tight tolerances 
required on the deck casting because of tape path require- 
ments. The buffer arm and speed sensor assemblies are 
removable. This allows for easy replacement of the head 
plate, speed sensor, and buffer arm assemblies. Each of 
these assemblies is designed to be modular and inter- 
changeable. There are no service adjustments on the HP 
7980A because any adjustments required are done at the 
factory. No routine maintenance is required* except for 
regular cleaning of the tape path. 

The tape path components perform several tasks in guid- 
ing Ihe tape over the magnetic tape head. They make sure 
the tape is wrapped around the head properly and consis- 
tently. They also guide the tape past the head at a precise 
angle (skew] and at a proper height (tracking) so that the 
written tape is interchangeable with other tape drives. Nor- 
mally, there are two precision stationary or fixed guides 
that perform these functions. The HP 7980 A tape path in- 
corporates one precision guide Into the buffer arm assembly 
and the other guide is integrated with the speed sensor 
roller. 

Speed Sensor 

The speed sensor integrates the functions of an optical 
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encoder and a precision tape guide. The particular chal- 
lenge was to create a rolling element, needed for transmit- 
ting tape speed to the encoder, that would guide the tape 
with the accuracy of a fixed guide. 

Fundamentally, a roiling guide is simply a vertical cylin- 
der with flanges on the top and bottom. The upper and 
lower flanges restrict the vertical limits of the tape's travel. 
According to the ANSI standard, the width of a half-inch 
tape must he between 0.496 and 0.500 inch. Clearly- the 
distance between the roller flanges must be at least 0.500 
inch not to damage maximum-width tapes, but this allows 
a 0.496~inch>wide tape to wander 0,004 inch between 
flanges. However, in the simplified tape path the speed 
sensor must also guide the tape for proper skew and track- 
ing across the head. This makes 0.004 inch of vertical wan- 
der unacceptable. 

A traditional solution of this problem is to put a fixed 
guide between the roller and the head. A fixed guide is a 
nonrotating cylinder with a top flange set to a precise 
height, A spring-loaded washer pushes the bottom edge of 
the tape so that the top reference edge of the tape stays in 
contact with the top flange. This eliminates tape wander 
since t ht; distance between the washer and the top flange 
i in vary with tape width. 

The HP 7980 A speed sensor combines features of both 
rolling and fixed guides. Fig. 5 shows a cross section of 
the speed sensor. The roller, which Includes a top flange, 
spins with the tape and transmits tape speed (through the 
shaft) to the optical encoder. The roller has no bottom 
flange, but there is a spring-loaded ceramic washer below 
the roller. The washer is restrained from rotating, but can 
move vertically to push the tape against I he fixed height 
of the roller's upper flange. Thus, the roller is fixed verti- 
cally, but moves rotationally, while the washer is fixed 
rotational I y and moves vertically. 

To keep the bottom of the spinning roller from rubbing 
on the non-spinning washer, the roller width is slightly 
smaller than the minimum tape width. Thus, the semicir- 
l ufir arc of tape wrapped around the roller pushes the 
ceramic washer down and away from the roller. The spring 
is placed under the center of gravity of this semicircular 
an of tape. This balances the forces on the washer so that 
it remaios horizontal and does not damage the bottom edge 
CtJ die tape. 

Buffer Arm Assembly 

The design of the buffer arm assembly presented several 



special challenges. The buffer arm assembly must provide 
tape buffering, proper tape tension, a servo position signal, 
proper tape guidance* and an overt ens ion shutdown signal. 
The buffer must be lightweight enough to maintain servo 
bandwidth, but strong enough to hold the tape height 
within a tenth of a millimeter. Because the HP 7980 A is 
an autoloading drive, the buffer arm must also act as an 
im to load the tape correctly. To achieve these require- 
ts» a thin- wall aluminum die-cast part is used. Its flex- 
ibility allows an air baffle, a spring post, a stop, and a slot 
for overtension shutdown to be incorporated into a single 
part. 

During testing, ii was discovered that the tape resonant 
frequency dropped when the drive repositioned frequently. 
It was determined that during such reposition cycles more 
air became entrapped in the tape stack. As a result, the 
effective tape length was increased, decreasing its spring 
constant. To overcome this problem, a Coulomb damper 
is used. A cantilever spring is placed between the buffer 
arm and fixed guide so that it produces approximately four 
inch-ounces of frictional torque. This dissipates enough 
energy from the system to allow servo stability without 
affecting other parameters. 

In a typical tape drive, the tape, after leaving the buffer 
assembly, enters a fixed guide for proper skew alignment. 
To minimize space requirements, simplify service, and re- 
duce parts count, the fixed guide is combined with the 
buffer assembly. The buffer arm pivots about the center of 
the fixed guide as can be seen in Fig. 6. The fixed guide, 
made of a stainless- steel ring and two ceramic; washers, is 
bonded on the buffer base to a u\015-inm tolerance. This 
allows the buffer assembly to provide all the tape guidance 
in the front half of the tape? path, but still be removed by 
loosening only ihree screws. 

Another feature is the way the position of the buffer arm 
is sensed. The sensing assembly had to be small enough 
lo pass through a 40mm-diameter hole in the head plate, 
hut require no adjustment if the buffer is replaced, To 
satisfy these requirements, a small ceramic magnet is 
mounted on the buffer shaft. Changes in the rotating mag- 
netic field are sensed by a linear Hall-elf etM IC. The plastic 
magnet holder allows each magnet to tie bonded in its 
correct calibrated position. With this arrangement a signal 
that is linear within ±8% over the buffer arm range is 
achieved, requiring no further calibration for interchange. 
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Design for Manufacturability 

The HP 7980A tape drife Is designed to be built easily. 
Tin- major improvement over the previous HP 7978B drive 
is its size, The HP 7980 A Is about a third the size and b -ill 
the weight of ihe HP 7978B, This makes il much easier to 
build and handle. 

All of the major components are attached to the main 
tiirk casting. The main deck casting is moved down the 
manufacturing line as major components and assemblies 
are attached to it, 11 is then flipped over to add components 
onto the top and tested, Most of Ihe components and as- 
semblies are attached to the casting In a downward move- 
ment and fastened with a common sell-lapping screw. 

Reduction of parts is another major Improvement. The 
number of mechanical parts was reduced from approxi- 
mately f510 on the HP 7978B to 340 on the HP 79ftuA. The 
number o! unique mechanical parts was reduced (rum ap- 
proximately 260 to 140, All these factors help substantially 
reduce the time to build the HP 7980 A. It is, therefore, a 
much more manufacturable machine than its predecessor, 
the HP 797HB. 
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Streaming Tape Drive Control Electronics 

by Bradfred W. Gulp, Douglas R. Dome!, Wayne T. Gregory, Jeffery J. Kato, Gerod C. Melton, Kraig 
A, Proehl, David W. Ruska. Virgil K. Russon, and Peter Way 



IN DEVELOPING THE HP 7980 A Tape Drive, the de 
team leveraged its expertise gained from the develop- 
ment of previous streaming tape drives. Design goals 
included a low factory cost, a small form factor, an auto- 
loading capability, and a selection of interfaces. To achieve 
these goals, advancements were made in many key areas 
of tape drive design. These areas included the drive archi- 
tecture, controller software, read/ write electronics, tape 
servo system, and front panel 

Three-Box Control Architecture 

The HP 7980 A Tape Drive incorporates a control architec- 
ture based on the three major functional blocks of a stream- 
ing tape drive: interface, buffer, and drive, Each functional 
block contains a microprocessor-based controller and data 
path electronics to hand] e unit il.it,*. V\g. \ shows (he three- 
box architecture ot the HP 798UA. 
Interface. The interlace block handles all communications 
with the host computer. The HP 7980A contains an HP -IB 
{IEEE 48ti) interface. Currently. tw T o other interfaces (FER- 
TEC and SCSI) have been developed. Drives with one or 
the other of these interfaces are available as the HP 88 781) A. 
Each interface accepts tape commands from the host. Inn in- 
fers the user data, and returns the appropriate Hiatus. To 
perforin a task, the interface fi rat relays the tape command 
to the buffer controller, User data is then sent to or obtained 
from the buffer memory. Upon completion of the task, 
status is obtained from the buffer controller. 
Buffer. The buffer block is responsible for streaming perfor- 
mance, The buffer controller oversees the use of I hi- Tj12K- 
buffer memory, The buffer accepts lape commands 
from the interface and forwards (hern to the drive control 
ler. The controller transfers user data between the interface, 
the buffer memory, and the drives formatter. Status is ob- 
tained from the drive controller and forwarded back to the 
interface. Streaming performance is maintained through 
the use of immediate response on write operations am] 
read-ahead on read operations, Immediate response on 
writes is implemented by the return of status as soun as a 



full record of daia is accepted into the buffer, but beh > 
tie tape. This allows the host to fn 
memory and begin transferring the next record. Wb 
write error occurs on the tape, the but ler will perform the 
necessary write retries independently of the Interface and 
host computer, Read-ahead cm reads is implemented 
issuing additional read commands to the drive so that 
streaming is maintained until either the buffer is full, a 
command other than a read is received from the luted 
or a read error Occurs, This trees the host and interface 
irom the constraint of performing all command and status 
processing overhead within the interblock gap time. Over- 
head processing time can be averaged out over several rec- 
ords. The buffer also performs all necessary read retries 
independently of the interface and host computer. Incorrect 
data is oot sent to t he host. 

Drive, The drive block executes all tape commands. The 
drive controller oversees the servo system, motor drive, 
tape path, formatter, read write, and front panel, The drive 
controller accepts tape commands from the buffer con t oil- 
ier and executes them. User read data is read from the tape, 
corrected and checked by the formatter, and then sent to 
the buffer memory* User write data Is accepted from the 
buffet memory, encoded by the formatter, and then written 
to the tape. Read -alter- write occurs on all write commands 
lo ensure error-free writes. Status is sent to the buffer con- 
l ml ler. In executing a tape command, the drive controller 
must control tape speed, tension, and position, The drive 
controller relies upon the servo system electronics for 
closed-loop control of tape speed and tension. Tape posi- 
tioning is performed in conjunction with the read elec- 
tronics and formatter and is based on block boundaries. 
The drive controller also manages user requests and infor- 
mation through the front panel. User requests requi triiiH 
buffer and/or interface involvement are passed to these 
functional blocks. 

Each functional block is connected by a common < gfn- 
munication link (COL) to an adjacent functional block. The 
CCL encompasses common hardware and common soft- 
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Fig. 1. HP 7980A architecture 
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ware definitions, All oi' the optional interfaces communi- 
cate with the buffer using the same CCL. The buffer com- 
municates with the drive using a similar CCL. The CCL 
software definition will be discussed later, 

The CCL hardware uses a message bus for the transfer 
of commands and status, and a data transfer bus for the 
transfer of user data to and from the tape. The message bus 
interfaces the controller of a master functional block (closer 
to the host) to a dual-port RAM located on the controller 
of a subordinate functional block [closer to the drive). Fig. 
1 shows the location of the dual -port RAM for each CCL- 
This dual-port RAM is directly and fully addressable by 
the microprocessors of both controllers. Commands and 
slat us are passed between the controllers by writing and 
reading specific memory locations. The data transfer bus 
interfaces the data path of a master functional block to 
FIFO ffirst-in t first-out) buffers located on the data path 
of a subordinate functional block. Separate inbound and 
outbound FIFOs are used to facilitate loop back testing. User 
data is passed along the data transfer bus in a sequential 
manner. Odd parity is maintained throughout with parity 
checking occurring as data is received across a functional 
boundary. An end-of-dala bit is asserted to flag the last 
byte of each data record. 

Interfaces 

Important to the tape drive's acceptance by OEM cus- 
tomers (as the HP 88780A) is a selection of interfaces. The 
drive currently supports HF-IB, SCSI, and PERTEC. In the 
development of these interfaces, common hardware and 
software designs were shared. All three interfaces contain 
a 6809 microprocessor- based controller, About 30 percent 
of the circuitry is common to all three interfaces, As men- 
tioned above, a common communication link definition is 
used to connect an interface to the buffer. Each interface 
completely isolates the rest of the drive from host, interface, 
and protocol specific requirements. By developing these 
interfaces in parallel it became clear what functionality 
belonged in the interface as opposed to the buffer or drive 
functional blocks. 

The development of three sets of interface code in the 
time allotted for one was a major undertaking. There ;ire 
24K bytes of HP-IB code, 40K bytes of SCSI code, and 17K 
bytes of PERTEC code. Design time was reduced through 
the use of a common firmware architecture and the leverag- 
ing of common code. These interfaces share the same top- 
level firmware structure. About 6K bytes of code is common 
among the interfaces. The common code contains dual-port 
RAM communication routines and command processing 
of front-panel requests including diagnostics- With the ex- 
ception of a power-on routine and a table-access module, 
all of I he code is written in high-level languages. The use 
of high-level languages increased design efficiency and 
code readability, Overhead processing time became an 
issue with the PERTEC interface. Instead of rewriting the 
code in assembly language, overhead time was removed 
by rearranging the order of tasks during reads and writes. 
An increase in parallelism between interface and buffer 
was obtained, saving about two milliseconds per command. 
The three-box architecture gives the HP 7980 A and HP 
88 780 A designers flexibility. Not only can new interfaces 



be easily designed and supported, but also options and 
enhancements can be added easily. For example, a new 
buffer design can be integrated without affecting ;my of the 
interfaces or the drive. Parallel development debug, and 
test of the functional boxes increased design efficiency by 
minimizing the impact t h h t problems in any one area had 
on the entire project. Because these boxes are functional 
entities, they could be turned on and debugged separately. 
Separate turn-on and debug included functional testing as 
well as the testing of communiciiHons with neighboring 
blocks, The integration of the first complete HP 79B0A 
produced a working drive within a month. 

Common Communication Link Software 

The common communication link. [CCL] is a single inter- 
face definition for two separate interface links. The HP 
7380 A contains three processors, a host interface control- 
ler , a buffer controller, and a device controller, performing 
specific tasks. The main function of the buffer controller 
is to intercept read and write commands and buffer the 
data within its data buffer to provide fast data access. For 
operations not involving data buffering or command queue- 
ing h the buffer controller typically passes the operation 
(command) directly through. By having common interface 
links on both the interface and the drive, the complexity 
of executing unbuffered commands is greatly reduced. 

The CCL is designed to have a common command set 
for the half-inch tape format. Although the HP-IB interface 
is the primary host interface for the HP 7980 A, the CCL 
command set does not preclude support of other interfaces 
such as SCSI, PERTEC, or a future high-performance inter- 
face. The buffer controller and drive controller implement 
a superset of commands that represents the operations 
needed to support most half-inch tape interface protocols. 

CCL protocol is implemented by a IK-byte dual-port 
RAM, This dual-port RAM is accessible by both the master 
(closer to the host CPU) and the subordinate (closer to the 
drive). Access is controlled by message control locations 
which synchronize the passage of messages across the CCL. 
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Fig, 2, Dual-port RAM use, 
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The dual-port RAM is divided Into two command ports, 
a diagnostic message port, a status area, two interrupt con- 
trols, and two interrupt locations (see Fig. 2). The command 
ports support bidirectional command transfers. A master- 
initiated command port passes commands from the master 
to the subordinate, while the subordinate-initiated com- 
mand port passes commands from the subordinate to the 
master. The master-initialed command port typically exe- 
cutes commands generated by the host, while the subordi- 
nate-initiated command port typically executes commands 
generated from the front panel. 

A command port contains one of two message types: a 
command sent from the initiator or the report returned by 
the target, The contents of the command port are deter- 
mined by a message control location. 

The diagnostic message area is used during certain diag- 
nostic operations when data messages larger than the com- 
mand ports [64 bytes) need to be passed. A diagnostic mes- 
sage can be up to 512 bytes long, followed by a two-byte 
checksum. 

The status area is used to maintain information on the 
drive that changes infrequently. This area Contains infor- 
mation such as tape loaded , write protected tape, and the 
density of the current tape. 

THh interrupt control locations are used to specify areas 
of the dual -port RAM that require that an interrupt be sent 
notifying the other controller of any change in their con- 
tents. Interrupts are most often requested on changes to 
the command port message controls (indicating that a com- 
mand is ready or a report has been returned) and changes 
to the status area. The use of programmable interrupts gives 
more freedom for each controller's implementation. The 
interrupt locations allow the master and subordinate to 
interrupt each other independently. The interrupt location 
is written with a value that indicates the area of dual -port 
RAM that has been changed , 

The CCL commands are divided intn six major opcode 
ranges: write, read, space, general h configuration, and diag- 
nostic commands. Each command range is further sub- 
divided into the individual commands supported within 
that group. For example, the write command group consists 
of write data record, wriln file mark, write gap, and write 
density IB. The opcode ranges simplify command decod- 
ing. Parameters on CCL commands are dependent on the 
opcode and vary in size from none to 57 bytes. For most 
opcodes, the number of parameter bytes is fixed. For ex- 
ample, the write record request opcode contains one three- 
byte parameter specifying the maximum length of the data 
to be sent, the write tape mark opcode contains no param- 
eters, and the write density ID contains a one- byte param- 
eter specifying the density to be written. 

The report type returned on a CCL command depends 
on the command sent. All tape commands return a normal 
report containing status flags, tape position flags, a byte 
count field, a retry count, and an error code. The stains 
flags quickly tell whether the command was successful, 
successful with error recovery, unsuccessful, or rejected. 
The tape position flags indicate the current position with 
regard to beginning of tape (BOT] t end of tape (EOT], and 
end of file (EOF). The byte count contains the length of a 
write record or read record. The retry count indicates 



whether physical retries were needed during the operation, 
and the error code contains detail on a command failure 
or command rejection. 

Non-tape commands return report types specific to their 
operation. For example, the diagnostic lest commands re- 
turn a diagnostic report that can contain the test number 
that failed, the failure code, and up to two field- replaceable 
units (FRUs) that could have caused the failure, 

High-Performance Buffer Software 

The HP 7980A is a high-performance streaming tape 
drive. It is also very flexible in its ability to be adapted to 
different environments and to use performance features in 
a variety of ways. High performance and flexibility of use 
in the HP 7980A are achieved by the features provided 
within the data buffer. 

A streaming tape drive mechanism such as the drive 
controller requires a constant stream of data if it is to main- 
tain high performance. In a typical host-drive interaction, 
irregular or slow data rates can interrupt the constant 
stream of data and cause the tape to stop and perform a 
costly repositioning cycle. The data buffer addresses the 
data rate problems by providing an asynchronous link be- 
tween the interface and drive, allowing each to operate 
independently at its own data and command rates for lim- 
ited time intervals. The HP 7980A has a large buffer in 
which to hold data as data rates are matched, making the 
asynchronous link possible. The buffer holds 512K bytes 
of data, enough not only to match variations in data rates, 
but also to allow the data buffer subsystem to accept data 
throughout a repositioning cycle at an average rate of 400 
kilo byte s/s e co n d . 

The data buffer packs variable-sized records into its buf- 
fer and places an entry in a queue indicating I he location 
and size of the corresponding record within the data buffer. 
The queue is large enough to allow the buffer to be com- 
pletely filled without reaching queue limitations except on 
small records. 

The buffer maintains two independent processes* one 
that interacts with the drive and another that interacts with 
the interface. Any transactions the host computer makes 
with the HP 7980 A Tape Drive are made with the interface 
side of the data buffer, The host transaction is complete 
when data has passed between the host and the buffer and 
necessary reporting and status information has passed. The 
operations necessary to write or read the data from the tape 
occur independently of the host interactions. 

During write operations data passes into the buffer from 
the host, As soon as data is in the buffer, the command is 
reported on with an immediate response report indicating 
1 hat the operation is complete. The queue contains the si^e 
of each record. Data held in the buffer is then written to 
the tape automatically when a sufficient amount of data 
has been accumulated or when an appropriate time-out 
has occurred. The drive is streamed as long as data is avail- 
able in the buffer. 

When a read command is received from the host, the 
buffer goes into a read-ahead mode during which the drive 
side of the buffer controller issues continuous reads to the 
drive controller, Data passes into the buffer and reports are 
placed into the queue. This time they contain not only 
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record size, hut also information regarding the success of 
the read attempt. Data is taken from the buffer and reports 
are taken from the queue to satisfy all read commands from 
the host. As data is removed by the interface, the drive tries 
to keep the buffer fid 1. managing tape motion as necessary. 

The queuing/buffering algorithm implemented in the HP 
7980 A allows a great deal of flexibility with the contents 
of the buffer. When data is read from the buffer, pointers 
are moved forward in both the data buffer and queue. Until 
the information is physically overwritten, it is possible to 
recover the information by moving the pointers back again. 
This capability is built into the read-ahead functions of the 
buffer controller- If a backs pace command is received while 
the drive has been performing read-aheads> the buffer first 
attempts to recover the data by moving pointers back. This 
electronic backspace, as it is termed, gives a great perfor- 
mance advantage in that it eliminates the need In reverse 
the tape physically to position back one record, 
increased Buffer Flexibility, The data buffer has other fea- 
tures that provide greater flexibility in the use of its perfor- 
mance capabilities. Two of these features are greater access 
to buffer functions and [he ability to change key performance 
parameters. Flexibility is also obtained through I lie use of 
manual buffer operations and through configurability; 

Buffer operations described to this point involve two 
sides ol the buffer The drive side performs either delayed 
or anticipatory actions automatically to provide quick re- 
sponse* to the host. Data is disposed of according to conven- 
tions in streaming conditions. 

Manual commands offered by the buffer provide the same 
capabilities that occur automatically during a read or write 
operation, w r ith each operation broken down into suhopera- 
tions* During manual operations, only the specified sub- 
operations are performed, with no automatic actions taken 
aild willi data disposition specified, not assumed. For 
example, a single write command is broken down into 
vvri te- to-buff er. write-to-tape, and remove-record -f rum- 
buffer suboperations, Manual buffer operations are much 
more cumbersome to use and are not optimized for high- 
performance streaming, but offer flexibility in certain key 
areas such as diagnostics and in recovering buffered data, 

Diagnostics make heavy use of manual commands, Data 
loopback operations are performed by writing to the butler 
and then reading data back to the interface. The drive pro- 
vides wellness and error rate tests using manual commands 
to w r rite and read tapes locally. Diagnostic access to manual 
commands also allows access to buffer and tape operations 
from the front panel. 

Recovering buffered data also uses manual commands. 
When a hard error occurs, or if the tape is stopped at EOT 
with unwritten data, the data in the buffer can be recovered 
by issuing read-from-buffer commands. Data recovery is 
not a typical operation, but is possible using manual com- 
mands. The flexibility of manual commands also allows 
current and future interfaces to implement interface spe- 
cific operations. 

The performance of the drive is dependent not only on 
the large buffer, but also on how the drive is set up to use 
the buffer. For example, parameters such as the trip point 
at which the drive starts writing data out of the buffer to 
the tape, the maximum record size to be written , or the 



lenglh of time writes are held can ail affect performance. 

The HP 7980 A is a highly configurable drive. Most vari- 
ables that can affect performance and many that determine 
(he personality ai tin; drive can be configured and main- 
tained Within nonvolatile RAM on the buffer board, Config- 
urations are used by the drive controller, the buffer control- 
ler, and the interface controller. A section of configurations 
is maintained and defined differently for each interface. 

The configuration values are distributed to the appro- 
priate subsystems when the drive is powered up. The inter- 
face also maintains default ROM configuration values. In 
the event that battery backup power is removed from the 
nonvolatile memory, default values are obtained from the 
interface. 

Through configurations, a drive can be set up with a 
distinct personality and set of performance features for a 
particular user. The configuration system also maintains a 
set of locks which can be used to lock the configuration to 
the current value and prevent it from being changed. 
Concurrent State Machines. The data buffer holds a special 
position within the drive in that it maintains CCL com- 
munications with four ports, three of them entirely inde- 
pendent. Each port has its own characteristics and needs. 

The downstream interface port is the source of all host 
commands. Commands received here must be processed 
quickly since all streaming operations into and out of the 
buffer occur here. Streaming commands also have the 
characteristic that when one is received, the next command 
will probably be the same. 

All streaming commands pass to the drive through the 
downstream device port. Quick response at this port is 
critical. If a command is not sent to the drive when a report 
is received, the drive controller may be forced to stop 
streaming and reposition the tape. 

The upstream drive and interface ports are not indepen- 
dent since there is no queuing of upstream commands. 
Commands received from the drive are either fully or par- 
tially processed by the buffer controller. If necessary, they 
are then passed up to the interface controller* Commands 
received from the drive are irregular and are typically ini- 
tiated as a result of human interactions. Hence, processing 
speed only needs to be fast enough for human response. 

Resets can also be received as asynchronous events from 
either the drive or the interface. Resets are infrequent and, 
like upstream commands, do not require fast processing 
response like downstream commands. 

The processing needs of the different ports led to an 
architecture within the buffer controller that allows 
maximum independence and specialization for the pro- 
cessing of each of the ports, Each major input source is 
handled by a separate program which controls the source 
and can be tailored to meet the streaming or background 
needs of each port. An operating system provides a set of 
concurrent processes for the programs to run within. Each 
program runs independently within its own time slice, 
with time apportioned for critical needs as they arise. 

The programs are implemented as state machines and 
the operating system as a concurrent state machine driver. 
State machines have the characteristic that they perform 
an action when a particular event occurs, then proceed to 
the next state and wait for another event to occur. State 
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machines can be controlled effectively by an operating sys- 
tem. The operating system performs state transitions and 
calls to state execution modules. It also allocates time 
appropriately between the different processes during the 
waiting periods. Within the time slice given to each pro- 
cess, each of the potential events for the current state can 
be checked, if an event has occurred, the associated action 
for that event is also taken. Actions are coded as straight- 
line operations without loops, At the conclusion of an action, 
processing returns to the operating system for state transi- 
tions and for other processes, 

The interface state machine handles all communication 
between the downstream interface port and the queue and 
buffer hardware. The drive state machine handles down- 
stream communication between queue and drive port, in- 
cluding all automatic retries of tape operations. The subor- 
dinate state machine handles upstream commands in both 
the interface and drive ports. The reset state machine han- 
dles all resets regardless of location. 

The operating system has no watchdog timer nor the 
associated tuning and communication race conditions as- 
sociated with unexpectedly losing the processor to another 
lime slice. Instead, the operating system allows each pro- 
cess to complete its action before passing control to the 
next process. Time for the heavy needs of the interface 
state machine and the drive state machine is automatically 
allocated by allowing the processes to take the time that 
they need. Processes such as the subordinate state machine 
and the reset state machine are effectively turned off by 
having only a single event to check on. They have no com- 



plexity until the event occurs. 

Each of the processes is independent, but each can have 
communication with other processes and affect their oper- 
ation by passing messages or commands. All communica- 
tion between processes occurs as events in the receiving 
process caused by actions in the sending process, Com- 
munication locations are tightly controlled d- jres. 
the queue being a prime example. The receiving process 
looks at the communication locations for conditions thai 
trigger events in themselves and that require action to be 
taken. 

The buffer controller functions are divided into indi- 
vidual programs, which are independent, have their own 
specialization and definition of tasks* and have localized 
access to hardware, ports , and internal communication lo- 
cations. This provides the structure needed to tackle the 
complex tasks of the buffer controller. Each of the pieces 
has defined functions and boundaries for good process defi- 
nition and implementation. 

Integrated Read Write System 

The design objective for the HP 7980 A read/write system 
was to lower cost while improving performance over pre- 
vious tape products. This objective w r as met through the 
reduction of printed circuit board space, the real- inn* "ii 
t imitation of read gains, the calibration of each read chan- 
nel, and the simplification of the write electronics. 

The HP 7 980 A read/ write board replaces the equivalent 
functions of earlier read T write t and write formatter boards, 
This reduction in printed circuit board space was achieved 
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Fig. 3, Analog read channel. 
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by using surface-mount components, combining the pre- 
amplification and differentiation functions into one stage, 
using en integrated disc pulse detector chip, and develop- 
ing a semicustum phase-locked loop chip. 

The read channel shown in Fig- 3 consists of a differen- 
tiating preamplifier, a gain-controlled amplifier, an auto- 
matic gain control circuit, a gap clamp circuit, a zero-cross 
pulse circuit, a zero-cross qualifier circuit, a calibration 
circuit, and a fourth-order Bessel filter. The input to each 
channel is the output of a single-track read coil. Two sig- 
nals P a zero-cross pulse ZCP and a zero-cross qualifier ZCQ, 
make up the output of each read channel. The zero-cross 
pulse is generated each time the differentiated input signal 
passes through zero. The zero-cross qualifier goes to a log- 
ical zero when the differentiated input signal becomes more 
negative than the lower hysteresis level and goes to a logical 
one when the differentiated input signal becomes more 
positive than the upper hysteresis level. This information 
is used by the phase-locked loop to recreate the digital data 
as it was recorded. 

Differentiating Preamplifier. The preamplifier is an NE592 
differential video amplifier configured as a differentiator, 
The 2 00- kHz [PE format) gain is set for 32 dB with a phasR 
shift of 88 degrees, and the 555-kHz (GCR format) gain is 
set for 40 dB with a phase shift of 86 degrees. 

A typical read channel has a preamplifier stage followed 
by an AGC amplifier and differentiating amplifier. To re- 
duce board space and parts count, the preamplifier and 
differentiating amplifier are combined. The layout of the 
differentiating preamplifier is very critical to achieving 
high common-mode rejection and retaining a high signal- 
to-noise ratio. Fig, 4 shows the difference between the con- 
stant gain and the differentiating configurations for the 
NE592. In the differentiating configuration, any leakage 
capacitance from the junction of capacitor Cl and resistor 
Rl results in the two emitters seeing different impedances 
to ground. This unbalances the differential front end of the 
NE592, resulting in a decrease in common-mode rejection. 
The board layout techniques shown in Fig. 5 are used to 
reduce the noise level and retain the common-mode rejec- 
tion of the NE592. The read and write sections have board- 
level isolation of 5V, ±12V t and ground. Each of the read- 
channel NE592 front ends has individual paths to the bot- 
tom ground plane, individual lines for power, and mini- 
mum trace runs near the differentiating preamplifiers. The 
head signals from the J2 connector to the NE592 inputs 
alternate between layers 2 and 3 which are located between 
ground plane layers 1 and 4. 

Disc Pulse Detector Chip. The preamplifier is ac coupled 
to the AGC amplifier, which drives a differential input, 
fourth-order, passive Bessel filter. The output of this filter 
drives the AGC, zero-cross pulse, and zero-cross qualifier 
circuits, The AGC circuit measures the input signal and 
compares it against the AGC level (an external dc voltage}. 
This voltage difference controls the gain of the AGC 
amplifier to make the peak-to-peak differential voltage out- 
put of the Bessel filler equal to four times that of the AGC 
level. The differentiating preamplifier turns peaks into 
zero -cross data so that a comparator can be used to locate 
the read signal peaks. If the read signal exhibits a tendency 
to return to the baseline between peaks, the comparator 



could respond to noise near the baseline. To avoid this 
problem, a zero-cross qualifier circuit is used. This circuil 
is a comparator with externally controlled hysteresis. The 
differential signal must go above or below the dc hysteresis 
value before the comparator will switch states* The com- 
parator output is fed to the phase- locked loop as zero-cross 
qualified data. The output of a bidirectional one-shot is 
fed to the phase-locked loop as zero-cross data that indi- 
cates the start of a data window. The phase-locked loop 
uses the zero-cross data to generate a pulse that is half a 
bit window long, This pulse is used to clock the zero-cross 
qualified data into a latch, The output of this latch repre- 
sents the data as it was put on the tape. 
Automatic Calibration. During track activity (data records), 
each track independently adjusts its gain via the AGC volt- 
age to an optimum amplitude, The AGC voltage required 
by each channel while reading data is determined using 
comparator Ul and DAC U2 as shown in Fig, 3, During 
gap conditions (no track activity] > the AGC voltage is over- 
ridden by individual gain control DACs, This allows each 
channel to be preset to a gain that produces an output 
signal of optimum amplitude if a nominal input signal is 
present. Detection of data by the phase-locked loop signals 
the channel to enter the AGC mode until another gap is 
detected or the gain increases above a predetermined value. 
DAC U2 and comparator Ul in Fig. 3 axe used to determine 
when the gain is above the predetermined value. During 
writes the gain is also held constant so that low-amplitude 
read signals cause a write retry. 

Track-to-track variation in the AGC voltage is dependent 
on the head and read channel combination and is indepen- 
dent of tape-to-tape variation. Track-to-track variation [gain 
profile) from a high-out put tape to a low-output tape can 
be closely modeled using a single gain offset value across 
all tracks. Gain profiles of GCR and PE densities differ, and 
therefore require separate gain profiles for each density. 
Automatic calibration is performed at the beginning of each 
tape when the tape density is being identified or being 
written. This calibration in GCR involves all nine tracks 
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Fig. 4. Constant gain (left) and dvidt (right) configurations 
for NE592 preamplifier. 
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during the ARA burst portion of the GCR density ID, while 

in PE it involves only the PE density identification track. 
Reference DAC U2 in Fig :i is stepped upwards until com- 
parator Ul indicates the reference voltage is equal to the 
ACC voltage for this track. This procedure is repeated two 
more times and the median voltage for * jved 

as one of the raw gain voltage levels for the current tape. 
The average gain of the raw gain voltage levels is calculated 
and saved as one of the average gain voltage levels for ihe 
current tape. A gain profile calibration is performed using 
the tape drive diagnostics. An all-ones pattern at the appro- 
pnate density is written for 100 feet beginning at BQT The 
pattern is read and the gain voltage for each track Is deter- 
mined 256 times. The average gain voltage for each track 
is the gain profile for the specific density. Each time the 
average gain voltage level is calculated during a tape ID 
routine, the gain profile for the specific; density is nor- 
malized such that its average value across all nine tracks 
matches the average gain voltage level for the current tape. 
Standard Ceil Phase-Locked Loop. When we decided to do 
a follow-on design to the HP 79 78 A Tape Drive^ our experi- 
ence with its design afforded us a great opportunity to 
make technological advancements in the areas of price, 
performance, and reliability. With this in mind, we selected 
areas of the HP 7978 A design where we could get the best 
return using some sort of LSI integration. One area was in 



lock recovery section of ihe HP 7978A. 

Designing a multidensity phase-locked Loop for a nine- 
track tape drive presents a problem when one is trying to 
minimize cost and maximize performance. However, with 
this as a design criterion, it was evident that some sort of 
id was needed. What we decided on was a digital 
analog hybrid phase-locked loop design. This gave us the 
flexibility to set important third -order loop parameters 
using normal analog loop parts, that is, an op amp, filter, 
and VCO, We then wanled to use LSI technology to inte- 
grate ail the digital logic associated with the recovery 
scheme. We chose HP's own CMOS standard ceil process 
whose density allowed us to put three tracks of phase- 
locked loop logii into one chip, and thus use only three 
chips per system. 

The major elements of the standard cell are shown in 
Pig. 6. Two programmable digital delay lines are used in 
the multidensity design, one for qualification of incoming 
read channel flux data, the other for calculating the phase 
error to be used by the phase detector. A pattern detector 
is used to detect the presence of data to control the phase- 
locked loop as well as the automatic gain control circuitry 
of the read channel. A two-state phase detector that can 
switch between frequency-lock and phase-lock modes vir- 
tually eliminates the chances of a harmonic lockup. An 
output encoder is used to synchronize all the channel out- 
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Fig, 5. The layout of the ground planes for the NE592 front ends ts critical for reducing noise 

and maintaining common mode rejection. 
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puts with a common syslern clock. The bus interface section 
contains fnur control registers for full microprocessor con- 
trollability- A center-frequency generator provides a refer- 
ence when there is an absence of data. 

The standard cell is implemented using HP's CMOS pro- 
cess. The implementation and testing of parts took only 
six engineer-months. Our first silicon version was produc- 
linn cjuality. Because the chip was flexible and usable the 
first time, we designed it into the HP 7976B as a cost-reduc- 
limi action before the HP 79H0A was introduced. The cost 
savings achieved in the HP 7978B was enough to pay for 
the chip's development cost. The gains achieved in the HP 
7980 A over the original HP 7978 A phase-locked loop in- 
clude a reduction in printed circuit board area by a factor 
of 10, an improvement in reliability by a Factor of 8, a 
reduction in factory cost by a factor of 5 ; and an increase 
in performance by a factor of 2< 

Digitally Controlled Write. Each of the nine data signals 
from the write formatter is fumed into two write currents, 
each driving one half of the write coil for one track. The 
positive halt (if (he write current is t:c.m I rolled by the signals 
GCR. PE, and CQMP while the negative half is controlled 
by the signals GCR\ PE', and COMP*. These six signals are 
generated for each channel by a PROM which is driven by 
common address signals comp, bp/g*, and dack. The other 
two address lines, data-1 and data-2, arc channel dependent. 
data-2 lags data-1 by one bit w r indow time period, allowing 
the PROM to know the value of the preceding data bit. For 
data equal to one, GCR is asserted for the full bit window 
time, PE is asserted for 100 percent of the bit window* time 
in PE mode but for only 16 percent of the bit window time 
in GCR mode, If the preceding bit was a zero, CQMP is 
asserted lor 80 percent of the bit window time in GCR 
mode or for 90 percent of the bit window time in PE mode. 
For data equal to zero the control signals GCR 11 and PE* are 
asserted instead of the signals GCR and PE. If the preceding 
bit was a one. COMP* will be asserted. 



Each track has two groups of three control signals. These 
signals control a group of three open -collector drivers. The 
open-collector drivers are connected to one end of a WTite 
coil through currant-limiting resistors. The center tap of 
each write coil is connected to 5 V. The half-coil current is 
then the sum of the three currents generated by each open- 
collector driver pulling its resistor to ground. This tech- 
nique allows the write currents to be generated using 5V 
TTL logic. 

To guarantee the write and erase currents will turn on 
only when commanded to, a POWER GOOD signal is hard- 
wired to o protection circuit Tfafl POWER GOOD signal is a 
signal from the pownr supply which will not be valid until 
volt ages are at rated values. The protection circuit turns 
on an enhancement mode n-channel MOSFET which 
supplies five volts to the center taps of the nine write coils. 

The same type of circuit is used in the erase circuit to 
turn on an npn transistor which provides a ground for the 
erase circuit. The dc erase current must ramp to the steady- 
state current in less than 0,2 ms and stay at this value until 
turned off. The erase current is controlled by two circuits. 
The two circuits operate in parallel until the steady-state 
current level is reached. The ramp circuit is then turned 
off, leaving the steady-state current circuit on until the 
erase circuit is turned off. 

Hybrid Servo System 

The HP 7980A servo system loads and unloads the tape 
automatically, controls the tape velocity and tension, and 
perforins diagnostic functions. The drive controller uses 
the servo system in conjunction with the read electronics 
and formatter to perform accurate control of tape position. 
The servo moves the motors independently during the tape 
loading process under the control of a 6809 microprocessor. 
Sensors monitor the progress of the tape and check for 
error conditions. Tension is established and the velocity 
and tension servos are enabled. The tape is under closed- 
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loop control thereafter. The tape velocity is primarily con- 
trolled by the take-up motor, and the tension is controlled 
only by the supply motor. A block diagram of the servo 
system is shown in Fig. 7. 

Closed -Loop Operation. Tape velocity is controlled by an 
_ral plus control scheme, This yields a 

ite velocity error. The long-term tape velocity 
error is well under one percent* The controller is Im- 
plemented as a digital-analog hybrid to take advantage of 
the best features of each type of control. In the HP 798QA. 
a design requirement for ihe servo is that the 6809 micro- 
processor be able to manage all the drive control functions 
and also control the velocity with minimal involvement. 
Fur this reason, the error signal and integral portion of the 
loop are performed by the processor, and the proportional 
feedback and multiplication by coefficients are performed 
b\ analog circuits. An optical position encoder is used by 
the processor to measure linear tape velocity. This is done 
by counting the number of pulses of the encoder every 5 
ins. The measured velocity is compared with a velocity 
command at each sample period and accumulated nu^r 
time. The resulting value is sent to a 12-bit UAC. The pro- 
cessor has complete control over the velocity of the tape 
by changing the value of the velocity command. This 
scheme performs the integral port ion of the control with 
high accuracy' The processor involvement is minimal be- 
cause onlv simple fixed-point additions and subtractions 
are required. The proportional part of the control is done 
in analog form. The optical encoder drives a tachometer 
circuit which gives a bipolar voltage proportional to tape 
velocity. The tachometer circuit does not have to be very 
accurate since it is only used to generate a proportional 
feedback term, The velocity control accuracy is determined 



by the integral term computed by the processor. The inte- 
gral and proportional terms are summed with appropriate 
scaling factors by an operational amplifier and the resulting 
signal is then sent to the motor amplifier. 

The output voltage of the velocity controller drives the 
take-up and supply motors. This reduces buffer arm deflec- 
tions during velocity ramps. If the velocity and tension 
servos were completely decoupled, the velocity servo 
'd start a ramp with the take-up motor, and the ten- 
sion arm would then deflect and start a ramp of the supply 
mi iU>t. This would require large arm deflections. In the HF 
7980A the velocity control voltage is scaled and se 
both motors. The buffer arm then only needs to deflect a 
small amount to correct the tape velocity at the supply 
motor. 

The variable-velocity rewind is an important feature of 
the HP 7 980 A tape drive because it allows a reduction of 
rewind time for a 2400-foot reel from 121) second's for a 
constant-velocity drive to 9D seconds. The maximum re- 
wind velocity attainable by a tape drive primarily depends 
on the maximum voltage available to drive the motors. 
Each motor generates a back emf proportional to its angular 
velocity. This back emf counteracts the supply voltage and 
limits the maxim lj m angular velocity the motor can 
achieve. 

The tape radius of this tape drive can vary over a 2:1 
ralio from beginning of the tape to its end. At each end 
one or the other motor can limit the maximum rewind 
velocity attainable to about 250 ips, In the middle of the 
rewind, the tape radius for both motors is roughly equal, 
and the tape can attain velocities of over 400 ips, tn the 
HP 79ft(JA, the processor maximizes the velocity of the 
drive during rewind by reading the voltages sent to the 
motors with an ADC, comparing them against the supply 
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voltage, and adjusting the servo velocity command several 
times a second to use the available supply voltage most 
effectively, Both servo loops operate at all times; only the 
velocity command is varied during the rewind, 

The tension loop measures the deflection of the holier 
arm and corrects the velocity of the supply motor until the 
arm returns to its nominal position. This effectively 
matches the linear tape velocities at both reels as the tape 
winds from one reel to the other. If both reels were identical, 
no buffer arm deflect ioo would be required. The buffer arm 
applies a relatively constant force of 10 ounces to the tape 
and thus passively controls the tension. The arm force is 
constant as the arm travels ±\ inch. The velocity and ten- 
sion loops are approximately decoupled since each motor 
sees a relatively constant load of 10 ounces. 

The tension loop is implemented completely in analog 
form. The buffer arm uses a linear Hall-effect sensor and 
a magnet to convert position to a voltage. This voltage is 
processed by an analog compensator and combined with 
the velocity control voltage to drive the supply motor* The 
key feature of the compensator is that it has an integrating 
term, This always returns the arm to its null point after a 
ramp, It also reduces arm deflections during ramps and 
eliminates final velocity-dependent arm deflections. If the 
integrating term were not present, the high rewind veloci- 
ties of the HP 7980A would not be attainable. 
Open-Loop Operation and Tape Sensors. The autoloading 
feature of the HP 7 980 A requires individual control of the 
two motors until the tape is loaded. The velocity and ten- 
sion control loops are each followed by a gain/load control 
stage, In the gain mode, this stage can output a voltage 
proportional to an 8-bit number. The processor can perform 
the autoload operations in this mode. When the loops are 
closed, these stages are switched to the gain mode. In this 
mode they operate as programmable gain stages. The pro- 
cessor can adjust the velocity and tension loop gains to 
account for different reel sizes. The gain/load stage is im- 
plemented with an 8-bit multiplying DAC. In the gain 
mode, the voltages generated by the servo loops are fed to 
the reference inputs of the DAC. and the processor adjusts 
the loop gain by writing to the DAC. In the load mode, the 
reference voltage is fixed, and the processor controls the 
motor drive voltage by writing to the DAC, 

Several sensors monitor the progress of the tape during 
the autoload process. A three-line encoder on the supply 
hub has several functions. It determines w T hether the feet 
on the supply hub have properly engaged the reel, it deter- 
mines the presence of the reel, it is used to control the 
velocity of the supply hub during the tape threading oper- 
ation, and it is used by the processor to determine the reel 
inertia. The three-line encoder is also used with the velocity 
encoder during closed-loop operation to determine the sup- 
ply reel radius instantaneously. 

The beginning-of-tape [BOT] and end-of-tape (EOT) sig- 
nals are determined on half-inch tapes by reflective mar- 
kers. The markers are detected by the tape drive using two 
reflective sensor assemblies. Some tapes are almost as re- 
flective as the markers, and this has been a problem for 
marker detection circuits. The HP 7980A employs a circuit 
that responds to changes in reflectivity, rather than absolute 
levels , to detect the markers. This means that the circuit 



is insensitive to component variations and can handle a 
wide variety of tapes without adjustments. The circuit is 
showm in Fig, 8. The sensor consists of an LED and photo- 
transistor pair. The phototransistor receives infrared light 
generated by the LED which bounces off the tape. As a 
marker passes the sensor, the phototransistor emits a cur- 
rent pulse. This current is converted to a voltage VI by 
amplifier Al, V'S is an attenuated version of Vl. V2 is a 
filtered version of VI. A2 is an amplifier used as a com- 
parator. A sufficiently large pulse at VI will cause a pulse 
at V {) since V2 does not rise as rapidly as V3. The output 
does not respond to slow variations in tape reflectivity or 
extraneous light sources, but it does respond to the rapid 
change across a reflective marker. 

A tape-in-path sensor determines when tape has entered 
the path. The sensor is a phototransistor and LED pair that 
shines infrared light across the tape entry point near the 
supply reel. The sensor is designed without adjustments 
by using a pulsed technique for detection. The detection 
circuit is similar to the circuit of Fig, 8. In this case the 
LED is pulsed by the processor for a short time with high 
current. If there is tape in the path, no output is detected 
at V[j, Otherwise, a fast change is detected at V [r 

An eight-port ADC is used to perform diagnostics on the 
analog circuits. It measures the buffer arm position during 
load and signals the processor to close the servo loops. It 
measures the motor voltages during rewind to adapt the 
rewind velocity, It is also used during self-tests to perform 
a motor-drive loopback, a tachometer loophaek, and a 48 V 
supply voltage check. 

Tape Positioning Control. Unlike a stop/start drive, stream- 
ing tape drives do not have the ability to stop and then 
restart the tape within the distance of an interblock gap. 
Each time the tape is stopped, a streaming tape drive must 
reverse the tape until it is well ahead of the next record to 
be read. The tape can then be ramped up to speed and the 
read electronics activated upon reaching the interblock gap 
where stopping was initiated. While stopping and starting 
the tape the read head may traverse an unknown number 
of records. Unfortunately, the 6250 GCR and 1600 PE tape 
formats do not support a standard mechanism for identify- 
ing individual records. Therefore 1 tape repositioning must 
be performed based on physical distance along the tape. 
The HP 7980 A Tape Drive uses an optical shaft encoder 
mounted to the capstan wheel in the tape path. This en- 
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coder provides two quadrature pulse trains. A quadrature 
decoder chip (QDC] converts these puises into a relative 
tape positioning count with the use of an up'down counter. 

Because of the high density and small gap size of the 
GCR recording format, it is important that relative tape 
positions supplied by the QDC be as accurate as possible. 
When a block-to-gap or gap-to-block boundary is sensed 
by the read circuitry* it is essential that the position ob* 
tained by the QDC at that instant be captured and stored 
away for later use. Although the QDC is a bus-compatible 
device, time delays from the occurrence of a block bound- 
ary to the point at which the microprocessor senses that 
event and reads the position can introduce unacceptable 
errors. To avoid the latencies caused by the additional mi- 
croprocessor cycles, position capture circuitry was de- 
signed around the QDC. This circuitry monitors inputs 
from the read formatter which signal when a block bound- 
ary has been crossed. When [he event occurs, external 
hardware extracts the position Information from the QDC 
and stores it in the appropriate location of a 4 x 12-bit regis- 
ter file. The position is captured within 2.4 pis of the event. 
This translates into a maximum position error of only 
0.0003 inch. 

Upon stopping the tape, the last valid block-to-gap posi- 
tion is used to mark the gap that follows the last block 
read. The tape is stopped, then restarted in reverse. Moni- 
toring the up/down counter, the drive controller waits until 
the read head is in the target interblock gap, The read elec- 
tronics is then activated and the tape position at the nexl 
gap-to-block boundary is latched. This position will be the 
same as the block-to-gap position previously latched plus 
any positioning error encountered in stopping and starting 
of tiie tape. The tape is then brought to a stop and restarted 
again in the forward direction upon receiving a new read 
command. Again, using the up/down counter to monitor 
tape position, the drive controller activates the read elec- 
tronics in the interblock gap before the next record to be 
read- By recapturing the bJock-to-gap position in reverse 
as a gap-to-block position, the drive controller prevents the 
accumulation of tape positioning error. Each time a posi- 
tion is recaptured, the amount of positioning error is calcu- 
lated by subtracting the two tape position counts. Ii this 
error becomes excessive it is reported to the host, thus 
indicating a problem with the drlvr. 



The VFD consists of an electron -emitting cathode filament . 
a screened phosphor anode, and a grid (located between 
the anode and cathode) to control electron flow. Therefore. 
the VFD is structurally similar to old radio or television 
vacuum tubes, yet it performs a totally different function. 
Besides the fact that VFDs are optical displays and not 
amplifiers, VFDs are often constructed with multiple triode 
structures included within the same evacuated tube. This 
feature allows the implementation of cost-effective, complex 
displays. Examples of these inexpensive yet complex \ 
can be found in all areas of modern consumer electronics. 

VFD segments are illuminated when electrons emitted 
by the heated cathode (filament) are accelerated by the 
positively biased anode and are allowed by the positively 
biased grid, or mesh, to collide with the screened phosphor. 
This causes photons to be emitted from the phosphor at a 
frequency that depends on the type of phosphor used. This 
feature allows the design of multicolored displays by screen- 
ing different phosphor patterns on the same display. AJ~ 
though this multicolored feature is attractive, it is not totally 
free. Standard available phosphors emit photons at lumin- 
ous efficiencies that can vary by as much as 50 percent. 
This causes a brightness balance problem between different 
phosphors on the same display- Brightness is typically bal- 
anced using optical filters, different power supplies, or 
different duty cycles 

The cathode filament, a thin oxide-coated tungsten wire, 
emits thermal electrons when heated to around B00°C by 
B ti lament power supply. The filament power supply is a 
major part of the power consumed by a VFD, 

The WD display used in the HP 7980A (Fig. 9] uses an 
ac filament supply providing 2 + 8Vrms and 0\47W at 715 
Hz, a imVde, 0.6W power supply which provides power 
far the anode and grid drivers, and a 40-pin IC which pro- 
vides level shifting (0-5V to OV30V] and digital control for 
1H anodes and seven grids. 

Why a Vacuum-Fluorescent Display? There are five reasons 
for choosing a vacuum-fluorescent display— cost, reliabil- 
ity, versatility, market separation, and human factors, The 
VFD is very cost-effective compared to other options avail- 
able. One other display technology, the liquid-crystal dis- 



Vacu urn -Fluorescent Front Panel 

We had two basic objectives in designing the HP 7980 A*s 
front-panel display. The first was to provide an HP-quality 
solution, especially since the display is a focal point for 
the user of the drive. We analyzed a number of factors, 
ranging from the angle of the display to its color and clarity. 
The second objective was to provide a convenient means 
of displaying tape drive status that would be viewable from 
across the length of an average computer room. Because of 
the long-distance viewing requirement, special attention 
was placed upon character size, color, and contrast, and 
upon the physical placement of annunciators. After a thor- 
ough analysis of the options available to us, we converged 
on a design employing a vacuum-fluorescent display. 
Display Operation, In simple terms, the vacuum-fluores- 
cent display (VFD) is a directly heated triode vacuum tube, 
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play, was slightly less expensive, but caused concern re- 
garding low-light viewability and electrostatic discharge 
sensitivity during assembly. 

Reliability is probably the VFD's greatest asset With a 
glass-encased body and an evacuated interior, it is imper- 
vious to electrostatic discharge damage and can endure 
harsh environments. Its low power consumption relative 
to LED displays also allows enhanced reliability. 

Since the VFD's light output is multicolored, it can be 
easily filtered to almost any color desired with standard 
colored lenses, In the HP 798UA Tape Drive, the blue-green 
output is filtered t<» almost white without significantly af- 
fecting the amber light output. 

Many other tape drives have LED displays and annun- 
cirifors. By using a more consumer-oriented display (the 
VFD is standard fare on VCRsJ, we were able to achieve a 
mare contemporary look. 

The contrast of the VFD's bright phosphor with the dark 
filtered background gives excellent viewability from up to 
twenty feet while not being too bright at close range. 
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as a production engineer, 
as a product designer for 
the HP 9876 Thermal Primer, and as a mechanical 
design engineer for the HP 7978A and 7980A Tape 
Drives. His work on the HP 7980A has resulted in 
i m .; f y-. j I e nt a op 1 1 c at ion s for ih e doo r a ssem b I y and 
the integrated tape path John lives m Greeley, Col- 
orado, and enjoys the usual Colorado outdoor ac- 
tivities of hiking, skiing, camping, and satJing He 
is engaged to be married this summer. 

Leslie G. Christie Jr. 

J Leslie Chnsiie was born in 
I Fayette vi lie, Arkansas, and 
attended the University ot 
Arkansas where he studied 
mechanical engmeermg 
(BSME 1970) and en- 
| gmee?mg science (MSES 
' 1979). Before joining HP, 
he taught meehnica 1 en- 
gineering courses at 
Kansas State Universifyand wasa design-engineer 
for a company that made fishing tackle He came 
to HP m 1984. contributed to the HP 7978A Tape 
Drive and worked on the butter assembly tor the 
HP 7980A Tape Drive. His work an the HP 7978 A 
Tape Drive resulted m a patent for an anti-adhesion 
apparatus for magnetic-tape drives Leshe fives in 
Greeley, Colorado, is married, and has two chil- 
dren. His free time is mainly taken up with family 
activities. 




43 _ Tape Drive Electronics I 
Bradfred W. Gulp 




Brad Gulp completed work 

for his BSE E degree from 
Virginia Polytechnic and 
State University in 1978. 
and with an HP fellowship 
earned his MSEE from the 
University of Illinois at 
Champaign-Urbana in 
1984 He joined HP In 1978 
and worked on the control- 



Wayne Thomas Gregory 

I- .. Tom Gregory pined HP's 

( Q^U^A Civil Engineering Division 
If m H m 1979. after receiving a 
^^^^L BSEE degree from Arizona 

^H^K State University, and 

worked as the production 
£ engineer far Ihe HP 3805, 
| 3B0B. and 3810 distance 
and angle meters. After 
transferring to Greeley Divi- 
sion R&D, he worked on the HP 7978 A Tape Drive 
and was responsible tor the power and display 
electronics for the HP 7980A Tape Drive. He was 
bom in Oceanside. California, and now lives with 
his wife in Fort Collins. Colorado. His outside in- 
terests include hiking, skiing, and !able tennis, 
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Gerod C> Melton 

Gerod Melton joined HP in 
1 9 /9 with a background in 
radar and navigation flight 
test instrumentation He 
worked first as a production 
engineer at the Loveland 
Instrument Division, ana 
then became en R&D en- 
gineer at the Greeley Divi- 
sion, where he was respon- 
sible for the design of the read'wnte electronics for 
the HP 79S0ATape Drive He studied electrical en- 
gineering at Texas Tech University (BSEE 1969) 
Gerod l vesin Loveland. Colorado. is married, and 
has one daughter His interests include hiking 
cross-country skiing, photography, and traveling 
with his family 

Peter Way 

Currently a product man- 
ager with HP's Greeley 
Storage Division, Peter 
J Way contributed to the 
A i servo design ana testing of 
1 \ - i-eHP797BA ( 9l44A.and 
„ ^ " 798QA Tape Drives He 

I joined HP m 1979 alter re- 
— ceiving an MS degree in 
aeronautics and astronau- 
tics from the Massachusetts Institute of Technol- 
ogy. He afso holds a BSE degree in aerospace and 
mechanical engineering from Princeton University 
(T977). Peter was corn m Caracas, Venezuela is 
married, and Irves in Fort Collins, Colorado. He «s 
active in local church and civic activities ana enjoys 
Cat racing, kayaking, and motorcycling 





Virgil K + Russon 

I Bom in Lehi, Utah, Virgil 
Russon studied mechani- 
cal engineering at Bngham 
| Young University (BSME 
\ 979 and MSME 1979) He 
pined HP in 1991 and 
I worked on the firmware tor 
I |heHP7978A,797BB.and 
7980A Tape Drives. Virgil 
: ivos in Greeiey. Colorado 
with his wife and five children He is a scoutmaster 
and active in his church. He enjoys commercial ana 
fine art and is currently landscaping his two-acre 
lot, 

David W. Ruska 

Born m Detroit. Michigan 
David Ruska rece;vec a 
BSEE from Michigan 
Technological University in 
1982 and joined HP the 
^jk ^"■b same year He developed 

jfS^S ^^^^ the rirmware} - 0r thG H?-\B 
interface of the HP 79 78 A 
Tape Dnve and developed 
the software architecture 
and firmware lor she dnve and front panel of the HP 
7980A Tape Drive Dav*d currently lives in Greeley 
Colorado H's interests include raising A radian 
horses Bible study, playmg keyboards, and ar- 
ranging and performing Messianic music 




SB 1= Digitizing Oscilloscope 

Joe K. Millard 

Author's biography appears elsewhere in this 
section 
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John J. Corcoran 

When John came to HP in 
■ 969 h ie |Oined an R&D 
team In the Sanf a Oara Di- 
vision doing laser inter- 
ferometer work He moved 
to HP Laboratories im 9 72 r 
O where he is now a project 
: manager He has been in- 
volved in integrated-circuit 
' design lor anak>g-to~digjtal 
converters and communications circuits, using 
bipolar, NMOS, CCD. and GaAs technologies. 
John hcids a BSEE degree from the University of 
lowaandan MSEE degree Irom Stanford He is a 
member of the IEEE and serves on the program 
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committee of the inseroatKjnai Sc^d-Slate Grcuils 
Confe- D 3ttm Valley. California, 

arid enjoys skwng. tennis m *U 

Ken Pouflon 

| 

-^Labofaio- 

H*og-tr> 

. ~1pO- 

d IC process tes$ 

vehsdes He came I 
- 
" MSEE and BS m physics 
^ ■ ' 

iTj? m £ - i h /// s * rg hi'^ :.--■-.'- ^eurs he 
worfceo at Ampex on SAW devices Ken is a 
member of the IEEE and has published a paper de- 
scribing the 1 -GHz ADC system featured 
issue of me HP Journal in the Journal of Sol »d -State 
Circuits and The Proceedings ol the ISSCC, Greuiis 
for analog/ digital conversion and CAD tools and 
methodologies are Ken's foremost professional in- 
terests He was born in Oakland. California and 
now lives with his wife and one-year-old son in Palo 
Alto Hjs favorite leisure activities are sknng and 
rock climbing 

Knud L. Knudsen 

■■ I in the 25 years Since Knud 

Knudsen joined the HP 

Laboratories, he has been 
working on a variety ol 
product developmenis. 

among them the HP T980 
Oscilloscope Measure- 
ment System, the HP 970A 
Digital Multimeter u 
HP &450 UV-Visibte Spec- 
trophotometer More recently. Knud has been re- 
sponsible for the development of ihe silicon digi- 
tizer chip described m This ssuect the HP Journal 
The chip is used in the HP 54201 . HP 54 1 1 1 D, and 
HP 541 12 Digital Oscilloscopes, and in an optional 
boa rd for t he H P 1 65Q0A Log i c Analyzer Knud was 
i Kin, i 5k|Ern Denmark, and his degree in elec- 
IroriiCb engineering (1957) is from Ihe Aarhus Elek- 
troteknikum there. Before moving to the United 
States, he did research work on radio receivers at 
the Technical University m Copenhagen Knud is 
named as an invenior or comventor on several 
patents lor transistor moise-measuring circuitry 
and a particle analyzer He is a member of the IEEE, 
and his technical interests locus on signal-to-noise 
optimization in measurement circuits, optical de- 
tection schemes, and solid-slate circuits rang ny 
from CCD sensors lo high-speed bipolar circuits 
Knud S married and has two grown children His 
avocations include playing the piano, making wine 
tenros r skiing, and soccer. 
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Lewis Dove came to the HP 

Colo tad o Spri ng s T ec h nol - 
ogy Center m 1982, after re- 
ceiving his BSEE degree 
ram The University at 
| Arizona Currently, he 
attends his alma mater lor 
a master's degree in sci- 
ence and electrical en- 
gineering on an HP resi- 
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Mark E. Mathews 

| Before joining HP's Col- 
I orado Springs Technology 

Center in 1964. Mark 
| Mathews had been in- 
solved in research on in- 
rared materials NMOS 
I process development, and 
S est engineering At HP, he 
handles froth product inn 
ji^* and R&D engineering for a 

variety of hybrid products. He earned his BS de- 
gree in physics at the University of California at 
Riverside ( 1 971 ) and his M SEE degree ai California 
State Polytechnic Universsty at Pomona {1976) A 
paper he wrote about PbTe/ZrO WIS structures for 
infrared sensing appeared m the Journal of Ape iec 
Physics. Mark was born in Hamilton, Ontario, he is 
married and has a daughter in his ofMiours. he 
plays the guitar for the Colorado Springs Jazz 
Ensemble 



Thomas K. Sohfey 

i Thomas Boh ley's BSEE 
and MSEE degrees are 
both from the University o- 
Missouri. Since joining HP's 
Colorado Springs Division 
in 1966. he has served as 
j engineer on the HP 1300 
I SenesX-Y Displays and as 
| project leader on HP 1 701 
1707, and 1980 Oscillo- 
scopes and the HP 1 332A Display Tom's primary 
interest is analog circuu design, and his work re- 
sulted in three patents for trigger circuits. He re- 
cently accepted a position with a medical laser 
company. Born m St Lours, Missouri, he presently 
lives with his wife in Colorado Springs. As an avid 
horseman Tom likes foxhunting and competes in 
dressage and combined training, but is also in- 
terested in skiing, hiking, and scuba diving 



Joe K. Millard 

For the past four years, Joe 
Millard has been projeel 
manager to* the HP 
541 1 1D Digital Oscillo- 
i scope at HP's Colorado 
1 Springs Division, which he 
joined in 1972 His earlier 
development projects fo- 
t cused on IC design for the 

HP 1 740 Series of Oscillo- 
scopes He's now a member ol the technical staff! 
at ihe Colorado IC Division. Born m Athens, 






Tennessee, Joe earned his PhD degree meteem- 

: neermg from the University of Ter 
kn 1 970, tes MSEE and BSEE are from the se. 

n Before he came to ^ ne worked on 
i^on at the Oak Ridge N a 
\Dfk has provide:; 
for some ii papers In off-hour: 
and four- wheeling the back roads in the Colorado 
mour -d and has two sons 



David W. Btgelow 

I Davd B geio w s BSEE de~ 
I gree is from the U^ 
I of Minnesota {1967} After 
receiving his PhD 
from Stanford tn *972. he 
I joined HP$ Microwave 
[ Technology Center m Palo 
Alto, California, where h»s 
respor uded 

thm-fiim process develop- 
ment He later moved to the Colorado Springs 
Technology Center, developing integrate:- - 
p roces s e s and serving a s I C process- en g i r ". I 
manager and as IC production manager. Presently, 
he is work-ng as a circuit design engineer A 
member of Sigma XL Dave was born m Rochester. 
Minnesota. He's married, has three children, and 
6 ac; ve in his church in Colorado Springs, He en- 
joys hiking to and from work — "except for the 
-nakes," as he emphatically points out. 

Donald D. Skarke 

I The product design of 
| oscilloscopes and digitiz- 
ing oscilloscopes, includ- 
ing the HP 180, 1726, 1980. 
I and 541 11D. along with 
CRTs and accessories .has 
been Don Skarke s mam 
professional interest Sffice 
| he pined HP in 1965 His 
1 S worfc on the HP ISO re- 
sulted in a patent Tor a plug-in grounding system. 
Don was born In Garwood, Texas, and received his 
BSEE degree from the University of Texas. He 
served as a sergeant in the US Air Force from 
1952 to 1956 Don has three children and three 
grandchildren and lives m Colorado Springs. Col- 
orado In his spare time, he enjoys building live 
steam models, bicycling and building bicycles 
and collecting watches 
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B. Allen Montijo 

I Asien Montijo came to Ihe 

1 HP Colorado Springs Divi- 
sion m 1981, after having 
f worked (here the previous 
summer as a student He 
receded his BSEE degree 
I from Princeton University in 
I 1981 and his MSEE degree 
I In 1985 Irom Stanford Uni- 
I varsity A design engineer 
in ihe R&D laboratory, Alien's professional interests 
include time-domain digital signal processing and 
the interaction ol subassemblies In system design 
He was born m Denver. Colorado, is marned, and 
has a daughter Allen enioys photography, skiing , 
and hiking. 
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A One-Gigasample-per-Second Digitizing 
Oscilloscope 

This instruments high sampling rate makes it particularly 
useful for analyzing high-speed, one -shot occurrences. A 
blend of state-of-the-art designs was required to achieve 
this performance, 

by Joe K. Millard 



OVER THE PAST FOUR YEARS, the demand for 
digitizing oscilloscopes has grown a I a rale of about 
30 percent per year. For CRT-based arming oscillo- 
scopes, this rate has dropped to about six percent, The 
growing preference for digitizing oscilloscopes is driven 
by several key advantages of their architecture. First, con- 
tinuous digitization of analog waveforms allows the user 
to view and analyze pretrigger waveform information. This 
feature is essential in diagnosing faults in digital systems 
that lead to a failure from which a trigger is derived. Second, 
total bus control of the digitizing oscilloscope. Including 
acquisition and transfer of waveform data, meets the grow- 
ing needs of automatic test and measurement systems, 
Third P waveform record lengths are not time-limited by the 
CRT display, but only by the available high-speed memory 
which can be the equivalent of many CRT screen widths. 
Until now, the primary disadvantage of digitizing oscil- 
loscopes has been the lack of large single-shot band widths 
which are essential for digital system diagnostics, high- 
energy physics studies, transient interference analysis^ and 
many other applications. The minimum bandwidth re- 
quired for a large number pi these applications is in the 



range of 200 to 500 MHz, 

The obvious problem in producing a large single-shot 
(frequently called real-time) bandwidth is the difficulty of 
sampling and analog-to-digital (A-to-D) conversion at rates 
oi about four times the desired bandwidth. Although the 
Nyquist criterion only requires a sample rate of twice the 
signal frequency for accurate reconstruction of sine waves, 
a ratio of about 4:1 (sample ram/handwidlhl is required for 
accurate reconstruction of arbitrary waveforms. Attempts 
to raise the signal bandwidth to values closer to the Nyquist 
frequency result in increasing overshoot and ringing of the 
reconstructed system step response as the Nyquist limit is 
approached, 

The HP 541 11 D Oscilloscope (Fig. 1] was developed to 
nd dress the high -bandwidth single-shot applications indi- 
cated above. The one-gigasample-per -second sampling and 
A-to-D conversion technology used in the product was de- 
veloped by HP Laboratories (see article on next page). The 
high-speed MM OS memory required for long waveform rec- 
ords at minimal power was developed at Colorado Springs 
Division and processed at the Loveland Technology Center. 
A general-purpose front-end attenuator and preamplifier 




Fig. 1. The HP 5411 ID Digitizing 
Oscilloscope can sample high- 
speed transient signals at 1 giga- 
samplers It features a 250- MHz 

sing (e -shot bandwidth and a 500- 
MHz repetitive bandwidth with up 
to 8 bits of vert tea I resolution. 
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were also designed at Colorado Springs Division [see article 
on page 67) using custom bipolar ICs and thick-film hybrid 
technologies from the Santa Clam Technology Center and 
the Colorado Springs Technology Center* respectively, De- 
scriptions of several of these developments will be pre- 
sented in the following pages. Digital signal processing 
techniques incorporated for optimizing signal quality will 
also be described [see article on page 70), 



A One-Gigasample-per-Second 
Analog-to-Digital Converter 

by John J. Corcoran, Ken Poulton, and Knud L Knudsen 



EACH CHANNEL of Hewlett-Packard's Model 5411113 
Oscilloscope contains a one-gigasample-per-secomJ 
anaiog-to-digital converter (ADC) which digitizes 
the signal from the preamplifier and writes the result to 
digital memory, The characteristics of this ADC to a large 
extent determine the horizontal and vertical resolution 
of the oscilloscope. The one-gigasample-per-second sample 



rate, for instance, directly determines the achievable 
bandwidth in single-shot operation, an impressive 25G 
MHz. 

Four main requirements shaped the architecture of the 
HP 54111D ADC: 
■ An 8000-sample memory depth was desired to allow ihi- 

capture oi an fight-microsecond time record al the full 
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Fig. 1.. A traditional flash analog-to-digital converter This ar- 
chitecture requires 2 N comparators and 2 N exciusive-QB 
gates 



Fig . 2 . An ana log encoder Encoding four reference con it ) uri 
sons tnto one analog voltage enables this ADC to be im- 
plemented with fewer devices and operate al lower power. 
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Repetitive versus 
Single-Shot Bandwidth 



There is often confusion among users of digging oscillo- 
scopes about the definitions of repetitive bandwidth and single- 
shot bandwidth, and how these two quite different specifications 
relate to the sampling rate of the ADC used in the instrument 

Single- shot bandwidth is important for those applications 
where the signaf occurs infrequently and It is necessary to digi- 
tize the entire event in a single acquisition cycle These applica- 
tions include those for which iradilional storage oscilloscopes 
were used m the past, and a host of others made possible by 
the digital properties of the new instruments 

Repetitive bandwidth applies when the signal occurs often 
enough that the waveform can be built up trom tens or hundreds 
of successive acquisition cycles to make a composite image on 
the screen. These applications include those where a traditional 
(nonstorage) oscilloscope was used in the past The key differ- 
ence is that single^shot bandwidth is proportional to the ADC 
sampling rate, and repetitive bandwidth is quite independent ot it 

Single-shot bandwidth is directly related to the ADC sample 
rate by the Nyqutsi sampling theorem This theorem states that 
to reconstruct a waveform accurately 1rom a limited set of sam- 
ples of that waveform, the sampling rate must be more than twice 
the highest frequency component in the waveform. It this condi- 
tion is violated, frequency components of the mpui l hat are above 
the Nyquist frequency (half the sampling frequency) will be trans- 
lated by the sampling process to different frequencies below the 
Nyquist frequency This frequency shifting (also known as alias- 
ing) causes errors that cannot be corrected m the sampled data. 
To avoid these errors, the signal bandwidth must be limited to 
less than half the sampling rate, either by its intrinsic characteris- 
tics or by an anti-aliasing filter in the signal path Then, to recon- 
struct the entire input signal correctly from the limited set of 
samples taken a digital tow-pass filter is used to interpolate 
between the actual sample points, Both the anti-aliasing filter 
and (his reconstruction filter must have Significant attenuation at 
the Nyquist frequency to avoid errors. As a practical matter to 
avoid excessive overshoot in the step response of these filters, 
their composite bandwidth cannot exceed about one -fourth the 
sample frequency Thus the filter requirements reduce the 
achievable single-shot bandwidth from the theoretical one-half 
of the sampling rate to the more practical one-fourth of the sam- 
pling rate 

Repetitive bandwidth, by contrast is totally independent of the 
A-to-D sampling rate In repetitive sampling, an image of the 
signal is built up by interleaving the samples acquired in multiple 
occurrences of the input signal. The proper placement (on the 
time axis) of each set of samples is guaranteed by a special 
circuit in the instrument that measures the time delay between 
the internal sampling clock and the trigger which occurs at the 
same point on the signal for each acquisition. The composite 
image thus built up has a very high effective sampling rate, and 
there is no need for either anti-aliasing or reconstruction filters, 
The bandwidth in this mode then depends only on analog com- 
ponents like the input attenuator, preamplifier, and sample-and- 
hold circuit, and not on the A-to-D sampling rate. 

John J. Corcoran 

Project Manager 

Hewlett-Packard Laboratories 



single-shot bandwidth. 

■ The power consumption of the analog-to-digital subsys- 
tem was required to be less than 25 watts per channel 
because of instrument power < nn.slr.iints. 

The resolution of the ADC had to be consistent with 
accurate waveform reconstruction. 

■ A 1-GHz analog-to-digital conversion bandwidth was de- 
sired so that the overall goal of a 500-MHz instrun 
bandwidth for repetitive signals could be readily 
achieved [see box at left for the distinction between 
single-shot and repetitive bandwidth]. 

The recj nirement for an 8000-sample memory system 
made it diffh n 1 1 lo use an indirect A-to-D architecture, that 
is. one th;it requires inlerrnediate analog storage. In an in- 
direct A-to-D system, the samples are captured at the full 
sample rate and stored in an analog memory [for example, 
a charge-coupled device). The samples are then read cut 
slowly and digitized by a low-speed ADC This architecture 
offers the attractive advantages of a lower-speed A-ln-D 
t;rter and digital memory. However, errors that in- 
crease with the storage time of the analog storage device 
and the power and complexity of the associated dock driv- 
ers limit the total practical storage capacity- He 
decision was made to build a direct ADC. where the input 
signal is digitized at the full sample rate and the data is 
stored directly in high-speed digital memory. The memory 
depth then can be increased without affecting the analog 
accuracy by simply adding digital memory chips. 

The power consumption constraint immediately made 
it clear that custom integrated circuits would be required 
throughout the A-to-D system, A board-level implementa- 
tion, using commercially available components, would 
have required too much power and probably could not 
have met the speed objectives. Next, each !C would need 
to be fabricated in a state-of-the-art technology with excel- 
lent power-delay product, whose performance characteris- 
tics are well-matched to the needs of the individual ICs. 
This led to the use of three different IC technologies: 

■ GalUum arsenide (GaAs) is used for the front end of the 
system, where the speed requirements are the greatest, 

- Silicon bipolar circuitry is used to digitize the signals, 
because of its excellent device matching and high trans- 
i I inductance. 

Silicon NMOS is used for the memory, since it offers an 
excellent combination of speed and device density * 
The silicon bipolar digitizer IC, which uses a modified 
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Fig. 3. A stmply interleaved A-to-D system Each digitizer is 
driven from a separate sample -and-hofd circuit 
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flash architecture, sets the resolution of the system, In the 
flash architecture there is one comparator for each input 
level to be resolved. Thus the transistor count and p 
are proportional to 2 X \ where X is the number of bits of 
resolution, It was decided early that a six -bit digitizer could 
both meet the power budget of the system and provide 
quite acceptable waveform reconstruction, especially when 
combined with some dithering and digital signal process- 
ing features (see the article on page 70), With these addi- 
tions, the effective resolution of the system is seven or even 
eight bits for slowly varying signals, where the higher res- 
olution is most useful. 

To achieve the 1-GHz A-to-D system bandwidth goal, it 
became quite clear that a sample-and-hold circuit would 
be required. A flash converter intrinsically samples the 
input signal in its clocked comparators, but the usable input 
bandwidth is typically much less than 1 GHz. For the digi- 
tizer designed for this system* the bandwidth for accurate 
conversion without an external sample-and-hold circuit is 
about 100 MHz. Thus, it was necessary to add an external 
sample-and-hold circuit t and investigation showed that 
only a GaAs FET IC could supply the needed bandwidth 
and acquisition time, 

Digitizer 

The silicon bipolar integrated circuit implementing the 
digitizing function for the system uses a flash architecture, 
which generally yields the highest A-to-D sample rate pos- 
sible within a given IC technology, However, even with 
HP's oxide-isolated 5-GHz IC process (the fastest bipolar 
process available at the time), a digitizer could not be built 
that could be clocked at the required one-gigasample/sec- 
ond rate. An even more severe constraint was the require- 
ment to drive and settle the input capacitance of the digi- 
tizer with the step-like input signals supplied by the sam- 
ple-and-hold circuit. Simulations showed that the best 
achievable settling time, with margins appropriate for a 
production instrument, was about 4 ns. Thus it was decided 
to build a digitizer that GOuld he clocked at 250 megasam- 
pl as/second (a 4-ns period), and interleave four of these to 
achieve the required net sample rate of 1 GHz. 

A block diagram of a traditional flash ADC is shown in 
Fig. 1. For six-bit resolution. 64 clocked comparators are 
required. The reference voltage V REF , equal to the full-scale 
input of the converter, is divided by the resistor network 
and distributed to the comparators. When the comparators 
are clocked, their outputs go to a logical one or zero state, 
depending on the relative value of the comparator reference 
and input signal at the time the clock occurred, resulting 
in a "thermometer code' 1 with all zeros above the com- 
parator nearest the input voltage, and all ones below that 
comparator. The thermometer code contains a single point 
where an adjacent one and zero occur, and this is detei ted 
by one of I he exclusive-OR gates. The high exclusive-OR 
output then drives the encoder to produce the correct bi- 
nary value, which is the output of the ADC. 

Although a Hash converter is very fast, it has several 
disadvantages. The device count and power are both high, 
because of the fully parallel nature of the architecture, For 
the same reason, the input and clock capacitances are high, 
and these lead to practical problems in driving the con- 



verter from other circuits. To minimize these problem- 
use a circuit that we call an analog encoder 12 to reduce 
the number of clocked comparators required. 

A block diagram of an analog encoder is shown in Fig. 
2. It consists of four differential pairs, connected with col- 
lectors alternately crossed as shown. The input signal 
drives one side and the re ! iltages drive the other 

side of the differential pairs. The resulting transfer function 
of this circuit, also shown in Fig. 2, changes between a 
high and a low level four times. These changes occur when 
the input voltage is equal to one of the reference voltages 
applied to the differential pairs. This circuit essentially 
implements the functionality 7 of four comparators and four 
exclusive-OR circuits in the traditional flash converter, ex- 
cept that the output is analog. Therefore, a single clocked 
comparator is placed after the analog encoder; The net 
result is that four differential pairs and one clocked com- 
parator replace four clocked comparators and four exclu- 
sive-OR gates. This is a net gain of about a factor of four in 
transistor count and a factor of at least two in power. Fur- 
thermore, since the number of (locked comparators is about 
four times lower, the clock input capacitance is similarly 
reduced. 

The si x-bit digitizer is implemented using seventeen of 
these analog encoders and seventeen clocked comparators. 
To produce the least-significant bits of the output data 
(where more than four transitions must be encoded), the 
outputs of several encoder-comparator groups, with ap- 
propriately interleaved thresholds, are logically combined. 

The digital encoder produces a Gray-coded output. Gray 
code is a method of coding 2 N levels in N bits. Its main 
characteristic is that a transition from any level to a neigh- 
boring level causes only one output bit to change. The more 
common binary coding has points where going from one 
level to the next can change all output bits simultaneoush . 
If a comparator in a binary digitizer should be undecided 
because its reference voltage is equal to the input, then the 
(mt put can produce large errors if it is interpreted differ- 
ently by the encoding paths for the different output bits 
affected. But in a Gray-coded digitizer, each comparator 
affects only one output bit, and hence, if a compamhir 
should be undecided, the output can only be interpreted 
to be one of two adjacent codes t which can be considered 
to be equally correct measures of the input voltage. 

Although analog encoded architecture reduces the de- 




Time t 



Fig. 4. Effect of uneven sample timing in an interleaved sys- 
tem 
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vice count, power, and clock capacitance, it does nothing 
to reduce the high capacitance at the analog input, which 
is typical ol flash converters, Simulations showed it would 
be impossible to drive and settle the expected inpul capaci- 
tance through the inductance of the bond wire used to 
connect the digitizer's input to its package. To eliminate 
this problem, an on-chip input buffer, optimized for step 
settling, was added to ihe digitizer. This buffer reduces the 
effective input capacitance from 25 pF to 4 pF and allows 
settling in the required 4 ns. 

The resulting digitizer chip 2 is quite small, measuring 
only 2.2x3,5 mm. This chip is a complete ADC in its own 
right; it can be clocked at rates up to 400 MHz* and 11 can 
be used without a sample-and-hold circuit on signal band- 
widths of up to 100 MHz. It is currently used in the HP 
54201 and HP 54112 Digital Oscilloscopes and in an op- 
tional board for the HP 16500 A Logic Analyser. Power 
dissipation of the digitizer is 2.5 watts, so the lour chips 
required for the 1-GHz A-to-D system consume a tola! ol 
10 watts. 

Interleaving Digitizers 

To meet the required sampling rate of 1 GHz, four of the 
silicon digitizers are interleaved, each operating at 250 
MHz. A simply interleaved A-to-D system is shown in Fig. 
3 t where each digitizer has its own sample-and-hold circuit. 
The problem with L his arrangement isthatthe four samplers 
sample the input signal separately, so any errors in their 
relative timing will be reflected as inaccurate sampling. As 
shown in Fig, 4, a deviation At in the desired timing of 
sample 3 results in an error AV in the voltage displayed 



for that sample. For high-slew-rate signals, such as a full- 
scale, 500-MHz input signal, maintaining six-bit accuracy 
requires that the relative timing be accurate to 5 ps. With 
the 1C and packaging processes available, timing errors of 
around 50 ps were expected- 

To avoid the 5-ps timing accuracy requirement, a fifth 
sample-and-hold circuit, operating at 1 GHz, was placed 
in front of the four sample-and-hold circuits operating at 
250 MHz (Fig. 5), This first-rank sample-and-hold circuit 
samples for 500 ps, acquiring the input signal on its hold 
capacitor. Then it goes into hold for 500 ps, which creates 
a flat spot on an otherwise quickly moving waveform. The 
four second-rank sample-and-hold circuits each sample for 
1 ns and hold for 3 ns. They will each finish their acquisi- 
tion time during one of the 500-ps flat spots. Since they 
finish their acquisition on a signal that is no longer moving, 
the timing of their clock signals is no longer critical. Pro- 
ducing a sufficiently stable clock lor the single first-rank 
sample-and-hold circuit is relatively easy. 

Sample-and-Hold Circuits 

The heart of each sample-and-hold circuit is the diode 
bridge sampling gate/ A simplified schematic is shown in 
Fig. 6. In the sample mode, Icontrol flows down through 
the four diodes , Dl-D4 t forming the bridge. If the voltage 
Vim is the same as V auTf the current is split evenly between 
diodes Dl and D3 and no charging of the + 3-pF hold 
capacitor takes place. However, if V [N is higher than V OUT . 
the voltages across Dl and D4 will be smaller than the 
voltages across D3 and D2. Hence, the current in D3 will 
be larger than the current in D4, which causes a net current 
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Fig. 5. Sample-and-hoid chip 
block diagram- Ail buffers are unity- 
gain source followers. The bridges 

are diode bridge sampling gates. 
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flowing into the hold capacitor until V FN and Vo&f are the 
same. 

The bridge is turned off (put into hold mode) by reversing 
the current I c untkui ■ In this case thr lows through 

diodes D5 and D6. which places a reverse bias on diodes 
D1-D4, breaking the connection between V tN and V, r 

A diode bridge sampling gate has several key advantages 
over simpler circuits such as pass-FET circuits: 

■ The circuit turns on and off very quickly. 

■ The diode % r oltages need only decrease by approximately 
200 mV (8V T ) to turn off the bridge effectively* 

■ The circuit is intrinsically balanced, so there is (to first 
order) no clock feed thro ugh onto V nuT . 

■ The hold isolation is quite good, Diodes Dl, D5 T and D;j 
effectively form a series-shunt-series isolation circuit in 
the off state, as do diodes D2, D6. and D4. 

A distinction is sometimes made between a sample-and- 
hold and a track-and-hold approach. If the charging of the 
hold capacitor is allowed to occur until V| N and V nirr 
equalize, the circuit is often called a track-and-hold circuit. 
If the on interval is shorter than the settling time, the circuit 
is called a sample-and -hold circuit Using CJaAsICs allowed 
us to design our circuit as a I rack-arid -hold circuit — the 
voltage on the hold capacitor settles completely to V LN dur- 
ing the 500-ps on time. The track-and-hold circuit has a 
higher gain than a sample-and-hold circuit and it is insen- 
sitive to variations in the pulse width of the sampling clock. 
In I he rest of this article, we shall use the terms track-and- 
hold and sample-and-hold interchangeably 

HP's GaAs [C process uses l-;im gate lengths to build 
depletion MESFETs (metal-semiconductor FETs) wilh 
a pinchoff voltage V JP of - 2V and a unity-cur rent-gain fre- 
quency f i of 14 GHz. Since a MESFET contains a Schottky 
diode junction, a fast diode (R ON C aFF * 1 ps) can be made 
from a MESFET by connecting its drain and source together 
as one electrode and using the MESFETs gale as the other. 
This process also includes integrated resistors and 
capacitors, A GaAs FET process is particularly well-suited 
to building sample-and-holrl circuits because: 

■ There are fast Schottky diodes for building diode bridge 




circuits. 

■ The FET speed allows construction of fast switching 
circuits for control of the bridges and very fast source 
follower circuits for buffering the inputs and outputs of 
the bridges, 

■ The semi -insulating GaAs substrate reduces parasitic 
capacitances. 

This process is described in more detail in reference 

The five sample-and-hoJd circuits shown in Fig. 5 are 
integrated on one GaAs chip. A companion clock genera- 
tion chip, also built in GaAs, accepts the input clock (at 1 
GHz or below), amplifies it to GaAs levels (3V swings), and 
provides the clocks for the five sample-and-hoid circuits 
and the four digitizers. Various characteristics of the chips 
used in this ADC are given in Table I. 





Table 1 






Chip: 


Sampler 




Control 


Digitizer 


i ■ Imology: 


(JaAsIl-MESFET 


GaAsD-MESFLT 
l -pLm Gate 


Si Bipolar 
5-GHzf, 


Chip Size: 


14 > 1.7 mm 




M-17 mm 


2.2X3.5 mm 


Device Count: 


460 




JHO 


920 


Number of Pads: 


4lJ 




35 


52 


Power: 


2.9W 




2.5W 


2.5W 



Increasing Resolution 

The major problem with an interleaved ADC is the need 
to provide good matching between the four digitizers, both 
in time and in voltage. The need for matching in lime is 
addressed by the use of two ranks of sample-and-hold cir- 
cuits. The need for matching in voltage is addressed by 
run trolling the voltage references used by the four digitizers 
wilh Individual digital to-analog converters that are set by 
the system processor during a self-calibration cycle. 

However, having multiple digitizers can he an advantage. 
The four 6-bit digitizers' reference voltages can he offset 
from each other in steps of 1/4 LSB 1 1 -ast-significant bit), 
that is, in steps of l/256th of the input voltage range. This 
causes the four digitizers to be inteHea% r ed in voltage, as 
well as in time. This allows 256 levels lo be resolved for 
slowly changing input signals, so the resolution for low-fre- 
quency input signals can beeight bits rather than six bits. 

To display a resolution of 256 levels, the data from all 
four digitizers must be combined. This can be done with 
a nonlinear smoothing filter, which moves each sample 
value towards the average of four adjacent samples, but 
never moves a sample by more than half a six-bit LSB (see 
"Dithering in the HP 54111D/' page 72). Although this 
mechanism can give true 8-bit resolution on sine waves 
only if the signal frequency is below l MHZ, it preserves 
the instrument signal bandwidth and provides a significant 
enhancement in resolution on the flat portions of pulse 
waveforms of any speed. 



Fig. 6. Simplified schematic of the diode bridge sampling 
gate Icohtroi ^verses direction to turn the sampling gate off 



System Performance 

The frequency response of the ADC, sampling at 1 GHz, 

I. .,■•■' ... ! , ■ 



JUNE 19BS HEWLETT PACKARD .mi.iHf ;.■ 63 



)Copr. 1949-1998 Hewlett-Packard Co. 



Digitizer Hybrid 



The three-inch-square digitizer hybrid circuit (see Fig 1 ) in the 
HP 541 1 1 D Digitizing Oscilloscope contains the GaAs track-and- 
hold chip, the GaAs clock and mode control chip, and Ihe four 
six 'bit, Gray code, bipolar anaJog-to-dsgital Hash converters that 
operate at 250 million conversions per second. Two main con- 
tributions tnat the thick-film technology makes to the circuit are 
high-frequency performance and high power dissipanon 




Fig* 1. The one-gigasample-per-second digitizer hybrid cir^ 
cuit used in the HP54U1D Digitizing Oscilloscope The GaAs 
sampie-and -hold and clock generation chips are in the center. 
The four 250- MHz bipolar digitizer chips are located in each 
comer. 

The high-frequency performance is made possible by mount- 
ing the IC chips directly onto the thick- film substrate and then 
wire-bond»ng the chip-to-substrate connections This eliminates 
the need for an intermediate level of packaging, which would 











m 



Fig. 2, Qose-up view of the two GaAs IC chips m the center 
of the digit ne r hybrid circuit of Fig. 1 

add delay and parasitics to the circuit Fig 2 shows a close-up 
view of the two GaAs chips mounted in the center of the substrate. 

The track-and-hoid and clock controller chips are wire-bonded 
directly together (chip-to-chip). This keeps the lead lengths be- 
tween them as short as possibte to minimize delay, the parasitic 
inductance of each lead, and the skew between the lines 

The circuit uses 50 and 75-ohm transmission lines for the 
analog input signals and clock lines, respectively. Thick-film tech- 
nology provides well-controlled impedances for transmission 
lines. It also makes possible the fabrication of very small resistors 
to minimize parasitics. 

The second contribution that thick-fitm technology makes to 
the circuit is in dissipating power — lots of power. The hybrid 
circuit dissipates 17 watts. The two GaAs chips, dissipating 6W, 
are attached directly to a heat spreader that fits through a note 
in the center of the substrate (refer to Fig. 2) This heat spreader 
has a thermal coefficient of expansion that nearly matches that 
of the 96% alumina substrate and has a thermal conduction value 
between copper and tungsten. Attaching the chips dfrectly to 
the heat spreader results in a decrease in junction temperature 




Fig. 3. Data acquisition board for 
the HP 54111 D. The digitizer hy- 
brid circuit is located beneath the 
shield in the center 
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of nearly 2£TC compared to mounting the chips on the aL- r a 
The chips are mounted to the heat spreader u&ng conductive 
•he epoxy process has a higher thermal resis- 
tance than a solder proce: npier 

The "nnectJons between the hybrid Circuit a r 

jft ooard are made using flex circuits Rex c 
consist of flat beryllium copper conductors etched on a polyimide 
~ne four flex circuit s: ase softfered around 

the perimeter of the substrate, and * r~ ds of the strips 

are soldered to the printed circuit board The flex circuits provide 
excellent high-frequency performance and alkDw stress rehef be- 
tween the hybrid circuit and the printed circuit board on which 
the hybrid is mounted 

A low-profile, large-area heat sink mounts on the back of and 
extends through the printed circuit board The heat sink contacts 



the bottom of the hybnd circuit at the heal spreader and directly 

beneath each of the -g-to-digjtai converters each of 

dissipates nearly 3W Five ceramic fids protect the hybrid 

s wire-bonded tC chips A shield covers the hybric 

and presses dc g thai good contact 

ffw heat sink and the Pack of the hybnd 
Rg 3 shows a picture of the data acquisition board Tne a 
hybrid circuit is unaemeath the metal shield m the centi 
the ftex circuit reads are just visible around the penmate' 
shield 

Lew/5 R Dove 
Mark E. Mathews 
Development Engineers 
Colorado Integrated Circuits C 



ted from page 63} 

is shown in Fig. 7. The 3-dB bandwidth is 1,7 GHz and 
the distortion ts below -40 dBc [13 LSB for B-bit resolu- 
tion) for input signals up to 500 MHz [The performance 
measurements quoted here are for the A-tu-D system alone T 
and do not include the HP 54111D*s input attenuator and 
preamplifier circuits,] Any timing misalignments between 
the four digitizers would show up as frequency components 
(sampling products] around 250 MHz; these components 
are suppressed to below -40 dBc by the use of the two-rank 
sample-and-hold architecture, 

Effective bits of resolution is o measure of analog-lo-dig- 
ital conversion accuracy:'' it measures an ADC's perfor- 
mance in digitizing a sine wave, compared against a 
hypothetical ideal ADC of that resolution. Any noise* non- 
linearity, or missing codes will cause an ADC's effective 
bits of resolution to fall below its nominal resolution. Thus, 
a perfect, noiseless, 6-bit ADC will show 6,0 effective bits 
of resolution; a poor nut' might show 3,5 effect iv6 bits. 

The resolution, in effective hits, achieved by this ADC 
while sampling at 1 CH/. is shown in Fig. 8. Here. 5,2 
effective bits are measured all the way up to a 1-GHz input 
frequency. The limiting factor in this case is the mismatch 
betw r een the four digitizers, since a single digitizer achieves 
5*9 effective bits of resolution at 250 megasamples/second. 
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Fig. 7. ADC frequency response and distortion 
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At lower input frequencies, the number of effective bits 
rises to nearly 7 as the higher resolution afforded by the 
nonlinear smoothing filler comes into effect. At low fre- 
quencies, the number of effective bits is limited mainly by 
the QOGlinearities of the digitizers. 

Fig. 9 shows the result of a beat frequency test using the 
voltage interleaving and nonlinear filtering described 
tibove. The larger sine wave is a 250-kH/ sinusoidal signal 
sampled at 1 GHz, The smaller sin<- ivave is the result of 
sampling a 1.000250-CHz sinusoidal signal. Even at a 1^ 
GHz input frequent y, the digitized waveform is clean with 
no visible aberrations. 

Packaging 

The A-to-D system (that is. the GaAs sample-and-huld 
chip, the GaAs clock generation chip, and the four bipolar 
digitizer (hips) is packaged on a single thick -film hybrid 
substrate (see bo\ on page 64). The two GaAs chips are 
side by side in the center and surrounded by monoblcN k 
bypass capacitors and thick-film termination resistors The 
four digitizer chips are placed near the lour corners ot the 
hybrid, a long with their associated capacitors and resistors. 

Acknowledgments 

We gratefully acknowledge I he contributions of Roman 
KagarHtsky, who built the test system for the prototype 




o 

o 
< 



• Single Digitizer 

A. Four Interleaved Digitizers 
with Smoothing Filter 



10 100 

fnput Frequency (MHz) 



1000 



Fig. 8. ADC effective bits of resolution Single ADC data rep- 
resents every fourth point from an interleaved I -gigasample- 
per -second record 



)Copr. 1949-1998 Hewlett-Packard Co. 



JUNE 1988 HEWLETT-PACKARD JOURMAL 65 





o 



u 

D 

< 



60 



50 



40 



30 



20 



10 



Input = 250 kHz 

Input = 1 GKz + 250 kHz 









Sample Frequency 
1 GHz 












12 3 4 

Time (/AS) 

Fig. 9. Reconstructed sine waves with resolution enhanced 
by interleaving the digitizers tn voltage. 

A-to-D system and wrote the software for it, Paul Clark and 
Doe Hi Her. who contributed to the design of the digitizer 
chip, Mike Kreikemeier, who designed the early versions 
of the hybrid packaging, Pauline P ml her and Dorothy Hol- 
lars, for amazing bonding work, the CaAs IC group at HP's 
Microwave Technology Division and the bipolar IC group 



at the Saflta Clara Technology Center, Joe Millard and his 
team at Colorado Springs Division, and Tom Ilornak, for 
his ideas, enthusiasm, and guidance. 

References 

1. T.W. Henry and M-P. Morgenlhaler, "Direct Flask Anulug-to- 
Digital Converter and Method*" US. Patent 4,388.339, Hewlett- 
Packard Company. 

2. J. Corcoran and K. Knudsen, "A 4110-MHz Gb A-to-D." iSSCC 
Digest qfTecfiliical Papers, February l!JH4 h pp. 294-295, 

3. JR. Cray and S.C, Kitscpoulos, "A Precision Sample and Hold 
Circuit with Subnanosecond Switching/ lEltli Transactions on 
Circuit Theory, September 1964, pp. 389-:jyfi. 

4. R.L, Van Tuyl, et al., "A Manufacturing Process for Analog and 
Digital GaAs Integrated Circuits,' IKF.E Transactions on Micro* 
wave Theory and Techniques, Vnl. MTT-30, no. 7, July 1982, pp. 

5. B. PeetZ, "Dynamic Testing of Waveform Recorders," IEEE 
Transactions an Instrumentation and Measurement, Vol. IM-32, 
no. l, March 1983, 

6. J. Corcoran, K Poulton, and T. Homak, "A I -GHz Mb A-lu-D 
System." IBBQC Digest of Technical Papers, February 19B7, pp. 
102-103. 

7. K + Poulton, J, Corcoran, and T. Homak, "A 1-GHz 6-bit A-to-D 
System, 1 ' IEEE journal of Solid-Stata Circuit Vol. SC-22, no. 6, 
December 1987. 



66 HEWLETT-PACKARD JOURNAL JONE IMS 



)Copr. 1949-1998 Hewlett-Packard Co. 



Front-End Signal Conditioning for a 
High-Speed Digitizing Oscilloscope 

by Thomas K. Bohley, Mark E. Mathews, Lewis R. Dove, Joe K. Millard. David W. Bigelow. and Donald 
D Skarke 



^^ RITICAL TO THE ACHIEVEMENT OF high-perfor- 
mance In any digitizing oscilloscope is the design 
^^^ of the front-end signal conditioning stages. The de- 
sign of these stages in the HP 5411 ID Oscilloscope was 
particularly difficult because of the HP 54111D*s 500-MHz 
required signal bandwidth. Key elements in the design are 
the attenuator* the switch assembly, the hybrid circuit upon 
which the attenuator is fabricated* and the preamplifier. 

Attenuator 

The theory of operation of a compensated one-megohm 
attenuator, shown in Fig. 1. can be found in any text on 
linear wave-shaping. Why then have they not been used in 
500-MHz oscilloscopes before? The answer can be found 
in Fig, 2; not in the divide-by-10 attenuator sections that 
can be cascaded, but rather in the straight-through divide- 
by-1 path. It is in that path that the inevitable series para- 
sitic inductance is at a maximum, and it is in that path 
that the series inductance and the input capacitance of the 
following preamplifier form a passive resonant circuit that 
introduces ringing on input pulses and limits the bandwidth, 

An effective solution to this problem is to raise the reson- 
ant frequency high enough so that the oscilloscope's ampli- 
fier no longer responds to that resonant frequency. Gener- 
ally, attempting to damp the resonant circuit by introducing 
a series resistance in the inductive path results in creating 
an RC corner in the passband of the following amplifier, 

To raise the resonant frequency we must make both the 
series inductance and shunt capacitance as small as possi- 
ble. The shunt capacitance, determined mostly by the input 
FET. is already small since a high-frequency FKT has been 
chosen. This leaves the series inductance to be minimized 
and the most effective way of doing that is simply tn liniil 
the physical length of the straight-through path. 

The attenuator is constructed on a thick-film 0.5 x 1.0- 
inch substrate with the input and output separated by about 
0.5 inch, To build an attenuator of this size would not be 
possible without one more technique — moving the ac cou- 
pling capacitor from the attenuator itself, where the capact- 
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Fig* 1 . Compensated attenuator circuit. 



in-*— <J 



xv 



Fig. 2* input attenuator with 
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Preamplifier 



-10 sectfons and straight- 



tor also contributes series inductance, to the preamplifier 
where it is placed in a low-frequency path. The di vide-by-10 
sections require considerable care to ensure that they, be- 
cause of their own parasitic inductance [Fig, 3), also do 
not form resonant circuits. In this attenuator, resistors R3 
and R4 are sized to damp the parasitic inductance associated 
with the capacitors and form a resistive divide-by-10 attenu- 
ator for signals well beyond the amplifier passband. Resis- 
tors Rl and R2 cam only relatively low-speed signals and 
are less important to the high-frequency performance of 
the attenuator than capacitors CI and C2, The divide-by- 10 
sections have greater bandwidth than the straight -through 
path and must be intentionally slowed down to keep the 
response as seen on screen consistent from range to range. 

Switch Assembly 

The challenges in designing the switch assembly were 
to provide small size, repea table characteristics, reliability, 
long life, and remote switching. The small size required 
that contact movement be minimal and repeatable. Other- 
wise, contact stress would he excessive and electrical 
characteristics would vary. Minimal contact movement 
also complicated the remote control requirement. The solu- 
tion to these problems w r as accomplished by the design 
shown in Fig. 4. Contact lift is controlled by the shuttle 
width, contact pin location, and dimensions o\ I in- J firmed 
contact, These parameters are controlled by the tooling 
used in their manufacture and thus are repeatable. Shuttle 
travel is not critical and is easily within the bounds of the 
solenoids and associated switching apparatus. The contact 
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Fig. 3, Resistors R3 and R4 added to circuit of Ftg 1 to damp 
parasitic inductances associated with capacitors CI andC2 
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Fig. 4. Attenuator swttch assem- 
bly for HP 5411 ID. 



pins are held and positioned by Ihe pin carrier and connect 
directly to the attenuator substrate to keep the elect rical 
pat It short. The switch contacts and contact puis are made 
of high-gold-content alloys for reliability and to reduce 
vulnerability to plating variations and wear. 

Attenuator Hybrid Circuit 

The attenuator hybrid circuit is shown in Fig. 5. This 
hybrid circuit has a switchable input impedance of one 
megohm or 50 ohms. There are three decades of attenuation 
al one megohm with 1:1, 1:10, and 1:100 divider ratios. 
Switch contacts are made as described in Fig, 4 using small 
pins that extend through the hybrid substrate, The elec- 
tromechanical actuators move a shuttle to open or close 
the appropriate contact. 

The attenuator hybrid circuit's bandwidth exceeds 500 
MHz. The high bandwidth is achieved by implementing 
the circuit elements in thick-film patterns and by 
miniaturizing the contact assembly, both of which mini- 
mize parasitic inductance and capacitance. 

The contact pins are attached to the ceramic substrate 
with conductive epoxy. A print process using a stencil was 
developed to apply epoxy automatically to the fifteen pins 
and into the pins' holes. To keep the epoxy from running 
completely through the holes to the backside of the sub- 
strate, positive air pressure is maintained on the backside. 
The pressure is turned off after the print squeegee passes 
over the holes so that it does not blow the epoxy back out 
of the holes. 

Most of the thick- film resistors are passively trimmed to 
a predetermined resistance value. However, one resistor is 
actively trimmed to set a circuit parameter. The resistor is 
used in the divider network of a circuit that senses danger- 
ous input overvoltages which could harm the 50-ohm T 1% 
input resistor. The divider resistor is actively trimmed by 
applying 10 volts to the input and trimming until 100 mV 
is obtained at the divider's output, Active trimming is typ- 
ically much less expensive than the alternative of manually 



adjusting a variable component. 

Preamplifier Hybrid Circuit 

The preamplifier hybrid circuit shown in Fig. 6 is mounted 
in the attenuator/preamp assembly. This amplifier boosts 
the signal from the attenuator to the amplitude required 
by the GaAs track-and-hold slage. The hybrid circuit has 
one-megohm input impedance, a bandwidth of 750 MHz, 
and an overall gain of 36 with a 50-ohm output load. It 
also provides trigger pickoff and conditioning. The 1.5 X 
3.0-inch hybrid circuit dissipates 5 watts. This hybrid con- 
tains a custom input MOSKKT, two custom bipolar lGs, 
and two pnp and one opn discrete microwave transistor 
chips. 

The custom MOSFET designed at the Loveland Technol- 
ogy Center includes a limited overload protection circuit. 
The first custom bipolar chip is a preamplifier that has a 




Fig. 5. Top (left) and bottom {right) views of attenuator hybrid 
circuit. 
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Fig, 7. Preamplifier hybrid circuit 
(right) and attenuator (left) mounted 
in the preampi attenuator assembly 



programmable l t 2,4 attenuator and a continuous vernier. 
This allows a user to adjust the input signal amplitude to 
full screen for maximum resolution. In addition, this chip 
generates a logic trigger signal used to reduce trigger jitter 
caused by noise contamination of the analog sync used in 
most instruments. Finally, an auxiliary input on the 
preamplifier chip can be switched on to calibrate the entire 
vertical deflection path from the input of the preamplifier 
chip through the analog-to-digital converters (ADCs), thus 
ensuring excellent gain and offset stability. 

The second custom bipolar chip is a driver chip which 
provides the final gain needed to drive the GaAs track-and- 
hold circuits. Included on the chip is a selectable x5 
amplifier. 

As a final part of the signal conditioning the differential 
signal from the ADC driver chip is level-shifted to near 
ground and converted to a single-ended signal, The circuit 
used is a pnp level-shift circuit, The final output voltage 



of the preamplifier hybrid is ±0,64 volt with a 50-ohm 
load. 

The preamplifier hybrid circuit is connected to the at- 
tenuator hybrid circuit through a compressed leaf spring 
that contacts gold pads on each of the hybrid circuits. Fig. 
7 shows the placement of the preamplifier next to the at- 
tenuator in the assembly (the attenuator is to the left of the 
preamplifier). The trigger output from the preamplifier IC 
drives 50-ohm differential transmission lines on the thick- 
film substrate. These lines lead to two SMC connectors 
epoxied to the back of the substrate. The connector pins 
protrude through holes in the substrate so that they can be 
soldered to the transmission lines on top. The linear signal 
output is connected to a third backside SMC connector 
through a 50-ohm rnicrostrip. All of the bias and ground 
connections to the hybrid circuit are made via two dual 
in-line connectors. These are soldered through holes in the 
substrate. 
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Digital Filtering in a High-Speed Digitizing 
Oscilloscope 



by B. Allen Montijo 



IN THE LAST FEW YEARS the digitizing oscilloscope 
has emerged as a useful tool for waveform acquisition 
and characterization of repetitive events. As a real-time 
(single-shot) tool, misconceptions about its abilities and 
limitations are restricting its use. Confusion over 'effective 
bits'' and ''interpolation" causes the user to question the 
validity of data and measurements. To add to the confusion, 
digital signal processing techniques support concepts such 
as trading excess bandwidth for extra resolution, 

Two events occur when a waveform is digitized; the 
signal is sampled and it is converted to a digital format, 
Sampling should be thought of as an analog process. This 
procedure transforms the signal from the s- plane to the 
z-plane and can create aliasing. Once aliasing occurs, no 
operation can recover the original signal without some 
knowledge about the signal. After the signal is sampled, 
an analog-to-digital converter (ADC) generates a digital 
code that represents the analog sample. This process is 
simply a format conversion, similar to converting a real 
number to an integer on a computer. The conversion adds 
quantization noise to the signal which is dependent upon 
the qualily of the ADC and the number of bits the ADC 
uses to represent the sample, The digitized signal can be 
thought of as the sampled signal with noise added. This 
article deals with the signal processing aspects of digitizing 
and assumes that aliasing is not a concern- With this as- 
sumption, the digitized signal contains all of the informa- 
tion necessary to reconstruct the original continuous signal, 
with some additional noise, Then the digital samples can 
be interpreted as the original analog signal with quantiza- 
tion noise added. This concept forms the basis for believing 
the information displayed by a digitizing instrument. So 
long as you know what the signal processing is doing to 
your signal and that aliasing is not a problem, you can 
believe w r hat is on your screen, 

Effective Bits and Noise 

Effective hits has become a popular term for describing 
the performance of digitizing instruments. The term com- 
pares the noise characteristics of the instrument with those 
of an ideal ADC. To find an expression for effective bits, 
begin with the equation for the mean squared error of an 
ideal ADC. The error can range from -1/2 to +1/2 of a 
quantization level (Q] with equal probability. This error is 
squared, integrated, and divided by the interval: 
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Q = A/2 N 
so that 

side = A 2 /(12 x 4 K ) 
The equation can be solved for N to give: 

N = Jog 4 (A 2 /s^ c ) - 1.7925 
For a full-scale sine wave input, the rms voltage is 

B = AJlVl. 
Then 

N = Log4 (B 2 /s* d J - 0,2925 
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To find effective bits for a system, the noise term s^ is 
replaced with the total system noises^ 

N,fr - log 4 [B 2 /s^ lBm ] - 0.2925 



The variable term (B 2 /s| ys!l > m ) is simply a sign a I -in- noise 
ratio [SNR), Effective bits is an expression of SNR; a factor 
of 4 in SNR (a factor of 2 in voltage) is one effective bit 
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If the ADC covers a voltage range of A and has N bits, then 



Fig, 1. Effective bits versus input noise for ideaf anaiog-to- 
digitai converters. The dB scale is for a fuil-scaie sine wave 
input. 



70 HEWLETT-PACKARD JOURNAL JUNE 19BB 



)Copr. 1949-1998 Hewlett-Packard Co. 



The SNR for an analog system can be plugged into equation 
7 for comparison with a digitizing system. 

For a digitizing system, the SNR is often found using the 
sine- wave curve-fit test, A spectrally pure sine wave is 
digitized, then the best-fit sine wave is calculated for the 
data- The best-fit routine finds the ideal sine wave that 
minimizes the rms error between the data and the ideal 
wave, This procedure cancels any systematic errors 
in offset, gain* frequency (sampling and signal), and phase, 
The residual rms error is substituted into the equation as 
the system noise voltage s 5yslem to obtain effective bits at 
the signal frequency. The test is repeated at various signal 
frequencies to obtain a curve of effective bits versus fre- 
quency. 

The test is not as straightforward as it sounds. Signal 
amplitude is an important parameter because the slew rate 
of the signal and the amount of distortion depend upon 
the amplitude. High slew rates put stress on the analog 
components and the ADC. and magnify the aperture jitter. 
Distortion is important because the harmonics generated 
are considered noise, For most statistical operations, more 
data is better. In this case, more data places stricter require- 
ments on the phase noise of the sampler and the test signal. 
Since the ideal sine wave is spectrally pure, any relative 
phase noise between the sampling clock and the signal 
decreases the number of effective bits- Finally, the purity 
of the test signal must be better than the unit under test, 
preferably by a factor of ten. This often requires the use of 
filters in the signal path, 

Care must also be taken in gathering data. The signal 
frequency should be selected so I hnt each output code ap- 
pears several times in the data; it is not a good test to input 
exactly the Nyquist frequency since only two output codes 
[ideally] will be produced by the ADC. 

The system noise being measured consists of two compo- 
nents: analog noise and quantization noise. Digital process- 
ing also adds a noise component, but this will be ignored 
since if is insignificant if the processing is done with enough 
bits. Analog component blocks, such as signal conditioning 
and sample-and-hold circuitry, create the analog noise 
while the ADC creates the quantization noise. An often 
neglected source of analog noise is the user's signal. With- 
out special processing the best SNR that can be hoped for 
is the SNR at the instrument input. Fig. 1 demonstrates 
how the number of effective bits is determined by both 
sources of noise. The x axis represents the rms noise at the 
input of the ADC relative to the ADC's voltage range (also 
expressed in dB), Feak-to-peak noise is the rms noise multi- 
plied by five. Noise is shown increasing to the right. The 



y axis displays the effective bits for ideal ADCs of four to 
eight bits. On the horizontal portions of the curves, system 
performance is limited by the ADC while on the slanted 
sections, the limitation is analog noise. All the curves over- 
lap when analog noise is dominant. This indicates that 
extra bits of resolution are of no value if they are only 
measuring noise. 

The overall signal-to-noise ratio, rather than the number 
of bits* is the important parameter. The effective resolution 
and measurement accuracy of a system depend on this 
value. More .ADC bits may give a prettier display, but may 
be meaningless as well as costly. 

Improving SNR by Filtering 

The seven and eight-bit filters in the HP 54 11 ID Oscillo- 
scope are low-pass filters that allow the user to trade off 
excess bandwidth for improved noise performance. Al- 
though this Is often done in analog systems, digitization 
adds complications. 

How do spectrum analyzers achieve a high signal-to- 
noise ratio? By limiting the bandwidth during a measure- 
ment, most of the noise is filtered out. allowing the detec- 
tion of small signals. For a wideband system, such as an 
oscilloscope, the major analog noise contributor can be 
assumed to be white in nature, that is, every frequency 
band of a given width contains the same noise power. This 
is reflected in the typical specification for noise t which is 
expressed as volts per square-root hertz, or in terms of 
power as watts per hertz. The noise Is present at all frequen- 
cies within the bandwidth of the system, whether or not a 
signal actually exists. For a reaMime sampled system, the 
maximum bandwidth is the Nyquist frequency f n . If the 
full system bandwidth is not needed for a signal H then the 
noise can be reduced by limiting the bandwidth, thereby 
increasing the SNR. If the noise is white! there is a direct 
correlation between bandwidth reduction and noise reduc- 
tion. SNR improvement by a factor of 4 gives one effective 
bit improvement; that is, filtering out 75 percent of the 
noise improves the performance by one effective bit. 

The linear model for the system is shown in Fig, 2, Analog 
noise is added to the signal by the preamplifier, then quan- 
NzatioM noise is added by the ADC. The quantization noise 
may be a function of the input, but the model is still linear. 
Superposition can be used to create the model shown in 
Fig. 3. Each component of the digitized signal can be 
analyzed easily. Assume thnt the filter ts a perfect low-pass 
filter at one fourth of the Nyquist frequency, If the input 
is within the bandwidth of the filter; it passes through the 



Input 



Input 



Signal 
Conditioning 



^ 



Analog 
Noise 



Processing 



Law-Pass 
Filter 



Quantization 
Noise 



^ Output 



Output 



Fig. 2. Unear model of digitizing system showing the sources 
of noise 
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Fig. 3. Superposition applied to the serial model of Fig 2 
results in this parallel model of an ADC system with noise 
sources. 
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Dithering in the HP 5411 ID 



In general, drthering is the addition of a small signal, possibly 
noise, to a system to improve measurement accuracy. For exam- 
ple, many DVMs average a number of samples to give a more 
accurate measurement Random noise helps the averaging pro- 
cess work, but a controlled signal provides better accuracy with 
fewer samples. 

Most ADCs have reference inputs that allow the user to deter- 
mine the voltage range for conversions. All of the output codes 
are referenced to this voltage range. If the references are raised 
by one volt, then the meanings of ail of the output codes are 
also raised by one volt That is, it takes one more volt on the 
input to get a particular code. Alternatively, one volt can be 
subtracted from the signal before the ADC to achieve the same 
effect. To get a particular code, the input signal must be one 
volt higher, The ADC acts as a buffer with differential inputs. The 
input minus the reference determines the output code Adding 
a signal to the user's signal can be accomplished by either 
summing the two signals before the ADC or by subtracting the 
signal from the ADC references. Therefore, the dither signal can 
be injected by subtracting it from the ADC references. 

The dither signal desired has increments of LSB/4 between 
samples to achieve an equivalent resolution of eight bits. With 
LSB/4 steps, the period of the dither signas should be bur sam- 
ples. The HP 541 1 1 D uses four ADCs sampling 90 degrees out 
of phase to achieve a 4 x improvement in sample rate Therefore, 
if the diirter works correctly, each ADC will "see" the same point 
of the dither signal each time it samples. Since the period of Ihe 
dither signal is 4/f s , all frequencies within the signal are at multi- 
ples of f £ /4. Each ADC >s sampling at fJ4\ so the dither signal 
spectrum appears at multiples of the individual ADC sampling 
frequency. All of these frequencies will alias to dc, To each of 
the ADCs, the dither signal appears to be a dc signal. Then the 
desired result can be obtained by adding the appropriate dc 
signal to the references of each ADC. This gives a stable, accu- 
rate dither signal for optimum performance. 

To prevent increasing the rms quantization noise, the dither 
signal must be subtracted from the ADC output code. The six-oit 
ADC code ts augmented with two bits so that the four dither 



values can be represented. With the final eight-bit code, the 
mean quantization and rms errors remain unchanged 

The order of the dither values is almost as important as the 
step size The pattern should give consistent results whether an 
edge is rising or falling The frequency content should also be 
spread out for the greatest attenuation by the filters (only two 
frequencies, f r , and f r /2. are avai&able with the given period) The 
d (her pattern at the ADC reference that best fits this requiremen* 
is 0, % v., 3& This gives dither bits of 00, 10, 01 , 1 1 Note that 
the farger part of the signal, the MSB, is changing at f n whiie the 
LSB is changing at f n /2. 

Dmenng improves ail modes of operation m, the HP 541 1 1D 
The seven-bit and eight-bit modes are described in the accom- 
panying article. 3n the six-bit (real-time) mode, a nonlinear smooth- 
ing filter improves the accuracy of low-siew-raie signals with only 
a minimal effect on fast edges. The filter averages four samples 
(two before and one after the current sample) If the average is 
wlthm LSS/2 (six-bit) of the current sample, then the current output 
is the average. Otherwise the current output is the current sample 
changed toward the average by LSB/2. Although the smoothing 
filter would work with random noise, the quality of the dither signal 
improves the performance 

In repetitive mode, drthering improves averaging by providing 
a resolution of LSB/4 with repeated acquisitions Persistence is 
the mode with the least improvement Dithering allows all levels 
of the screen to be filled in and increases the accuracy of some 
measurements made by the user. For example, an eye pattern 
is a multivalued function, so averaging mooe does not work 
correctly, When the user measures a portion of the waveform, 
the dithering provides LSB/4 resolution. 

Dithering could be carried further, providing LSB/B or LSB/16 
resolution with a longer signal period, However, since real ADCs 
are not perfect, a practical limit exists for any ADC. The ADCs 
rn the HP 541 1 1 D are tested to LSB/4 accuracy, or one LSB at 
eight bits. With this specification, dithering to LSB/4 resolution 
is an improvement over LSB/2. but LSB/8 would not provide a 
significant advantage (particularly since the HP 541 llD's four- 
ADC scheme accomplishes LSB/4 Oither^ng so easrly). 



filter untouched. The analog noise is white t so it is reduced 
by a factor of 4. 

For many signals, quantization noise will appear to be 
white noise, uncorrected with the input- For the moment, 
assume that this is the case. Then the filtering process will 
reduce the quantization noise by a factor of 4. Since both 
noise components have been reduced by 4 f the SNR has 
been improved by a factor of 4, or one effective bit. Simi- 
larly, the eight-bit filter has a bandwidth of one sixteenth 
of the Nyquist frequency and improves the SNR by a factor 
of 16, or two effective bits. 

Statistically, this can be explained as follows. The filter 
is an FIR [finite impulse response) filter. The filter performs 
a weighted moving average on the input to find the out- 
put, Let the filter have N coefficients, where N is odd, and 
let M = (N - l)/2. The filter coefficients are h if where 
i = -M t , . ., M, To achieve unity gain at dc: 



n: 

x. 



(8) 



The filtering equation for a symmetrical filter is; 



1= — M 



(9) 



A simple example of this would be the filter defined by 
the coefficients 0,25, 0.50, 0.25. With input x and output 
y, the value for y k is 

y k - PlJX^ 4 (t/2)x k + (l/4]x k ^ [10) 

If the filter input x is the noise, the noise at the output is: 



s* = var(y k ) = var % hpc^ 
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where var is variance and gov is covariance. The noise 
characteristic for all samples is identical t so Sj = s in . If 
the noise is white, all noise samples are independent, so 
covfs,, Sj) = 0. The equation reduces to: 



4- §„*?«?■ 



- v hf 

- -M 



(12) 



The output noise power is the input noise power times 
the sum of the filter coefficients squared. To gain an effec- 
tive fait, this sum must be less than 0.25 + This is not difficult 
to achieve since all of the coefficients have a magnitude 
less than one and most are positive. The filter given above 
reduces the noise to (1 4) 2 + (12)" - [1 4) 2 = 37.5 percent 
of its original value, However, the filter has poor passband 
characteristics so that the signal is also reduced at all fre- 
quencies but dc; the SNR is improved by 8/3 only at dc. 

When is the overall system noise white? This occurs 
when either white analog noise or white quantization noise 
dominates. With either of these cases, the improvement in 
SNR is straightforward. Quantization noise is generally 
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Fig. 4. interpolation by a factor of 4, (a) Example signal spec- 
trum before sampling, (b) Spectrum after sampling at fre- 
quency f s . (c) Sampled data pants, (d) Adding three zeros 
between pants shown in (c). (e) New spectrum (f) Desired 

low-pass filter (g) Final result. 



considered to be white whenever the signal changes by at 
least a few quantization levels between samples. In this 
case, the error from one sample to the next will take on a 
random pattern and the error signal will appear to be un- 
corrected. In general, the more bits in the ADC. the higher 
the signal frequency, and the larger the signal amplitude, 
the more random the quantization noise. 

Subharmonics can cause problems creating white noise, 
depending upon the order of the subharmonic and the filter 
length. For example, an input at one eighth the sampling 
frequency fV will give at most eight different codes (in the 
absence of analog noise). The quantization error for every 
eighth sample is correlated. If the filter has a length longer 
than eight, then the covariance term in the noise equation 
is nonzero. The worst-case signal for generating white 
quantization noise is dc with no analog noise. For this 
signal, the ADC outputs the same code for each sample 
and the quantization error is exactly the same (this is the 
degenerative case of a sampling frequency subharmonic]. 
The quantization noise signal is a dc signal — certainly a 
poor approximation of white noise. For a dc signal, the 
covariance term in equation 11 exactly cancels the increase 
in the SNR. 

In the HP 5411 ID, this problem is solved by adding a 
dithering signal to the analog signal before the ADC [see 
the box on page 721, The signal is precisely designed to be 
three fourths of an LSB peak to peak with its entire fre- 
quency content in the stop band of the low-pass filters, 
Sampling takes place so that four successive conversions 
sample the dither signal at LSB/4 intervals, The six-bit code 
from the ADC is augmented with two bits to cancel the 
error from the dither signal. Far a dc signal, four successive 
samples give a resolution of LSB/4 for eight bits of resolu- 
tion. For example, assume that an ADC has a Q level of 
one volt and that the code is expressed in volts such that 
code N represents N volts and will be present whenever 
the input voltage is ^(N-0.5) volts and <[N + 0.5J volts. 
If the input voltage is 1.5V, then the ADC outputs code 2 
for an error of 0,5V. If dithering is added, then analog 
voltages of (», -0.5, -0.25, and -0,75 volt are added to 
successive samples. The voltages that the ADC sees are 1.5 t 
1, 1.25, and 0.75 volts so that it outputs codes 2, 1, 1, 1. 
Note that the average value of the samples is 1,25V, The 
average value of the dither signal is -0.3 75V, so the average 
value of the signal is 0.375 volts higher than indicated, 
That is, the measured voltage is 1,625 volts for an error of 
0.125V. 

This example gives the maximum error for both cases, 
This is a factor of 4 in voltage, 16 in power; or two effective 
bits! If we augment the output codes with two bits to cancel 
the dither, then the outputs are 2. 1,5. 1*25, and 1.75 for 
an average value of 1,625 volts. When the outputs are av- 
eraged (filtered), the extra two bits neither help nor hurt 
the resolution; the extra resolution comes from providing 
enough extra code thresholds over four samples to match 
the number of thresholds for an eight-bit ADC, 

In the frequency domain, the dithering signal pushes 
most of the noise at dc (15/16 of it) to the dithering frequen- 
cies. Since these frequencies are in the stop band of the 
filters, the noise is eliminated from the output. Examination 
<it the dithmng characteristics in equation 1 I shows that 
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the co variance term is negative so that the overall improve- 
ment is better than with white noise- 
Under conditions where white noise dominates the 
input, the dithering signal has no effect; it simply appears 
as a high-frequency signal that is filtered out by the low- 
pass filter. Hence, the dithering has taken the worst-case 
signal and turned it into the best-case signal without hurt- 
ing the original best case. A signal at any frequency is 
improved in relation to the quantity of noise that is corre- 
lated. 

The filters used in the HP 5411 ID do not come close to 
ideal low-pass filters (see the box on page 75), The goal 
was to optimize SNR improvement while maintaining a 
good step response. The desired SNR improvement deter- 
mines the minimum stop-band width: 3f n /4 for the seven- 
bit filter and 15f fl /lb for the eight-bit filter. The stop-band 
attenuation is not critical; above 20 dB it is easier to improve 
the SNR by widening I he stop band than by increasing 
the attenuation, (Infinite attenuation is not much improve- 
ment when the major source of noise is in the passband,] 
The maximum amount of overshoot and ringing on the 
step response is known. These factors determine the width 
of the transition band and, in turn, the passbaud. To provide 
the required SNR improvement over as wide a frequency 
band as possible, the passband is as flat as possible. The 
bandwidths of the seven-bit and eight-bit filters are 100 
MH/, and 25 MHz at full sample rate and they provide one 
and two extra effective bits to frequencies of 83 MHz and 
18 MHz, respectively. 

The seven-bit and eight-bit filters in the HP 541 11 D, 
along with dithering, provide increases of one and two 
effective bits at the expense of bandwidth. In Fig. 1, the 
filters move the operating point at least one or two ADC 
lines up as a result of reducing quantization noise and one 
or two tics [6 or 12 dB] to the left as a result of reducing 
analog noise. Note that any noise in the user's signal is 
included whenever measurements are made and that the 
filters reduce this noise too. The HP 5411 ID has better 
effective -bits performance in the seven-bit filter mode than 
it would have with an eight-bit ADC without filters, al- 
though a penalty is exacted in lower bandwidth. 

Interpolation Process 

One of the most questioned functions of digitizing oscil- 
loscopes is interpolation. Users feel uncomfortable with 
the "fictitious'" data that is placed between samples to fill 
in the waveform. Most interpolators are linear functions 
that can be described in the same sequence of steps. This 
includes common filters such as sin(irx]/TTX [i.e., sincfx}), 
Gaussian, and Nth-order polynomials that fit the sample 
points exactly (LaGrange interpolation), A specific example 
of interpolation by a factor of four will demonstrate the 
Steps in the frequency domain. 

For simplicity, assume that the analog signal has B trian- 
gular spectrum as shown in Fig, 4a, After sampling the 
signal, the spectrum is repeated at intervals of f s . the sam- 
pling frequency, as in Fig. 4b. Mathematically, any combi- 
nation of the spectrums could account for the sampled 
waveform. The purpose of interpolation is to increase the 
sample density as though the signal were sampled at a 
higher rate. The first step in interpolation is to increase 



the sample density, Zero- valued data points are inserted 
between the actual data points (Fig. 4c) to get the desired 
density. For the example, three zeros are added to increase 
the density by a factor of four [Fig, 4d). in the frequency 
domain, this simply places the new sample rate. T h , at 4f s 
and the new Nyquist frequency f r ', at 2f 4 [Fig, 4e)> There 
are now four images of the original spectrum below the 
Nyquist frequency. This means that a combination of spec- 
trums is required to produce the zero-filled signal. The 
SNR is 1/3 h since there are three undesired spectrums to 
the one desired spectrum. The final step is to filter the 
result to leave the original spectrum at a higher effective 
sample rate (Figs. 4f and 4g). This low-pass filter is often 
a sine function or a Gaussian function. 

The most difficult concept of the process to understand 
is the spectrum of the zero-filled signal. A simple example 
may clear up the confusion. The signal 

J /:, + cos(2irf s t)/2 + cos(47rf s t)/4 

can be used to modulate the original signal before it is 
sain pi t-d (Fig, 5a), If the signal is sampled at frequency f^ 
in phase with the modulation signal, the result is identical 
to Fig. 4d. The spectrum of the analog signal before the 
digitizer is the convolution of the original signal spectrum 
(Fig. 5b) with the spectrum of the modulation signal (Fig. 
5c) n giving the 4 'zero-filled 11 spectrum (Fig. 5d], After sam- 
pling, the energy from 2f & to 2>5f s aliases back to I he range 
2f s to 1.5f b . doubling the height of the short triangle. Sam- 
pling also causes this new r spectrum below 2f s to be copied 
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Fig. 5, (aj Multiplying input signal by modulation signal. Orig- 
inal spectrum (bj convolved with modulation signal spectrum 
(c) gives a "zero-filled" spectrum (d) The result after sampling 
a! £ is identical to Fig. 4e. 
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Digital Filters 



The digital filters in the HP 5411 1D are FIR (finite impulse 
response) filters The FIR Kite*" has no unsampled analog counter- 
pan since a has no potes — only zeros. One of the advantages 
of FIB filters ts that they can easily be designed with perfectly 
linear phase, providing a constant defay through the filter regard- 
less of frequency. Constant delay is important for making accu- 
rate measurements Without constant delay, the appare 
between two waveform edges could be a function of the differ- 
ence m edge speed For example, in a digitaf circuit this might 
be a function of logic family or signal loading 

One characteristic of linear phase is that the step response is 
ant i symmetrica I, the preshoci is a mirror image of the overshoot, 
in an oscilloscope, this appears as noncausal oehavjor since 
the filtering routine takes the constant delay into account when 
placing the filtered points on screen. Although this is unsettling 
when first encountered, it is an advantage over an analog filler 
since the user can often separate the matching overshoot/pre- 
shoot response from the actual signal overshoot 

Another characteristic of the digital filters is that the bandwidth 
is related to the sample rate When the sample rate changes, 
the bandwidth changes, Some users consider this an advantage 
and some a disadvantage. Generally speaking, the sample rate 
is not lowered until the screen shows only one or two pertods of 
a waveform. In ihrs situation, the lower bandwidth does not ap- 
preciably alter the signal, but still provides the increase in effec- 
tive bits. For example, at 100 ns.'div. the sample rate drops to 
500 megasamples per second Now the seven-bit filter has a 
bandwidth of 50 MHz and a rise time of approximately 7 ns. The 
largest change in response occurs for a perfect step, where one 
sample is an the bottom of the step and the next is on the top, 
In this case, the filtered rise time will be 7 ns, or less than a tenth 
of a division. The rise time for the eight-bit filter will be about 30 
ns, or one third of a division 

The FIR filters achieve high throughput by using lookup tables. 
A table exists for each filter with each coefficient value multiplied 
by all possible data values At run time, the microprocessor need 
only perform lookup and add operations. The sizes of the tables 
were minimized by taking into account properties of the filters 
and data The linear phase characteristics are achieved by mak- 
ing the filter symmetrical about the center coefficient. This allows 
a reduction in size by a factor of two. The eight- bit data values 
are shifted to a two's complement format to give 128 nonzero 
magnitudes instead of 255, Altogether a data reduction of nearly 
four to one was accomplished Further reduction is possible, but 
would tower the throughput of the instrument. Throughput could 
have been increased by using a math coprocessor or a digital 
signaJ processing chip, but this was not done because of cost 
and space considerations. 

Many of the limitations to digital filters parallel limitations to 
analog filters . The characteristic of most concern ts that the steep- 
ness of the transition region between the passband edge and 
the stop-band edge determines the amount of overshoot and 
ringing m the step response, The filters in the HP 5411 1D were 
designed with this in mind, trading off bandwidth for pulse fidelity. 



at intervals of 4f s (fj), giving a sampled spectrum identical 
to Fig. 4e. 

The requirements for the HP54111D's interpolation filter 
are more stringent than for the seven-bit and eight-bit fil- 
ters. The filter can interpolate up to a factor of 100, For 
this case, the SN r R before low-pass filtering is 1/99, or - 20 
dB, The filter slopband must, attenuate the 99 spectral rep- 
licas to achieve an SNR of 42 dB. This minimizes the reduc- 
tion in effective bits and the measurement uncertainty in 
both time and voltage, Since there is no anti-aliasing I 
in the instrument (to allow a 500-MHz repetitive band- 
width), the stopband must begin near the Nyquist fre- 
quency of 500 MHz, To maintain good pulse fidelity, the 
maximum amount of overs boot and ringing is determined. 
As with the other filters, these quantities determine the 
width of the transition band and the passband. The filter 
is down 2 dB at 250 MHz and 4H dB at 515 MHz with only 
6% overshoot. 

The sinc(x) and Gaussian filters often used are extreme 
types of filters. The sinc(x) filter gives the maximum 
bandwidth (f ri ), but has a poor step response. The Gaussian 
filter adds no overshoot to the signal, but requires a very 
wide transition hand, giving an unnecessarily narrow pass- 
band. The HP 54111D f s interpolation filter is somewhere 
between these two extremes. One advantage to having a 
bandwidth of f 5 /4 rather than ij2 is that the transition region 
provides some anti-aliasing protection. Any signals in the 
500-to-750-MHz range will alias to the 25(Mn-50U-MHz re- 
gion, The attenuation of the transition band will perform 
an anti-aliasing function for this frequency range. Beyond 
750 MHz the attenuation of the analog signal path provides 
aliasing protection. 

A side effect of having a bandwidth lower than the 
Nyquist frequency is that filtered points do not always 
overlay the actual samp hit J points. Many users believe thai 
this is a shortcoming of the interpolation, although it is a 
necessary phenomenon. Any movement in samples results 
from reduced bandwidth after interpolation; no one would 
expect the bandwidth-limited trace on an analog oscillo- 
scope to toot like the full-bandwidth trau 1 To h'riw I In: 1 
actual samples in place, the filter's impulse response must 
have the exact zero crossings of the sinc(x] filter and its 
characteristics will approach those of the sinc(x) filter. It 
is very limiting on the filter design and fidelity to enforce 
such a requirement. 

Interpolation is not a difficult process, but must be done 
with care. The result is consistent with the response of the 
interpolation filter so long as aliasing is not a concern. 

Summary 

Digitizing a waveform includes two separate effects; sam- 
pling ami format conversion, Sampling is an analog process 
that transports the signal from the s-plaue to the z-phme. 
resulting in the possibility ot aliasing. Thf* format conver- 
sion simply changes the signal to a form that can be handled 
by a computer, adding noise in the process. If aliasing can 
be ignored, the digitized data contains all of the information 
necessary to reconstruct the signal, with a little added 
noise. For any measurement system, the signal-to-noise 
ratio is an important parameter. With digitizing instru- 
ments, this is more important than the actual number of 
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bits or the resolution of the display, Since digitized data 
is merely a different representation of a waveform, the same 
techniques that have been used to improve the SNR in 
analog systems can be used with digital systems. The seven- 
bit and eight-bit filters are digital equivalents of the 
bandwidth-limit button on traditional analog oscillo- 
scopes, 

A function that has often been questioned by digitizing 
oscilloscope users is interpolation. Waveform fidelity is 
highly dependent upon the quality of the interpolation 
filter. There is a direct trade-off between the bandwidth-to- 
sampling-rate ratio and the quality of interpolation, The 
HP 541 11D does not present the best filter for all types of 
signals, but it does provide a high-bandwidth, good-quality 
filter that works very well for the types of signals viewed 
most often on oscilloscopes — sine waves, pulse trains, and 
logic signals. 

The digital filters in the HP 541 11D provide benefits that 
are not easily available in analog counterparts. First, the 
filters have no delay at any frequency and timing measure- 
ments are not altered because of nonlinear phase charac- 
teristics. Second, a stable system response is achieved, not 
only from unit to unit, but with time and temperature. 
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